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By exploiting in large-A^ YM the change of variables from the gauge connection to the ASD part 
of its curvature by a non-SUSY version of the Nicolai map, we show that certain twistor Wilson 
loops supported on a Lagrangian submanifold of twistor space are localized on lattices of sur- 
face operators of Z^ holonomy that form translational invariant sectors labelled by the magnetic 
<~i . charge k = 1,2, ...,A^— 1 at a point. The localization is obtained reducing the loop equation in 

the ASD variables in the holomorphic gauge, regularized by analytic continuation to Minkowski 
space-time, to a critical equation, by exploiting the invariance of the v.e.v. of twistor Wilson loops 



> 

(•^ \ by deformations for the addition of backtracking arcs ending with cusps on the singular divisor 

^N ■ of surface operators. Alternatively the locaUzation is obtained contracting the YM measure in the 

tH- ■ ASD variables on the fixed points of a semigroup that acts on the fiber of the Lagrangian twistor 

r~^ ■ fibration which twistor Wilson loops are supported on and leaves invariant their v.e.v.. The renor- 

^^ . malized effective action induced by the localized YM measure in the ASD variables scales accord- 

ing to a large-A^ beta function of NSVZ type that reproduces the first two universal perturbative 
coefficients. Because of a non-trivial Jacobian due to the lack of supersymmetry a multiplicative 
is^ _ renormalization by a Z factor of the ASD field occurs. The masses squared of the fluctuations of 

^ ■ surface operators in the sectors labelled by k, supported on the Lagrangian submanifold analyti- 



cally continued to Minkowski space-time, form a trajectory linear in k that does not include any 
massless state. The glueball propagators in the holomorphic/antiholomorphic sector defined by 
correlators of a complex combination of the ASD curvature and its adjoint saturate at short dis- 
tances the logarithms of perturbation theory by a sum of pure poles. The anomalous dimensions 
of long gauge invariant operators belonging to the holomorphic/antiholomorphic sector that are 
implied by the Z factor coincide with the anomalous dimensions of the scalar operators that occur 
as the antiferromagnetic ground state of the Hamiltonian spin chain in the thermodynamic limit, 
that it is known to provide the anomalous dimensions in the ASD one-loop integrable sector of 
large-A^ YM. In this framework Regge trajectories of higher spins are related to fluctuations of 
surface operators with pole singularities of any order. 
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1. Introduction 

A technical presentation of the main ideas and results of this paper is in the synopsis and in a 
very sketchy way in the list of contents. 

This section is more of introductory nature, in order to convey in a simpler way part of the 
meaning of the technical ideas. 

The problem of the Yang-Mills (YM) mass gap as reported in [[T]] has an infrared and an ultra- 
violet nature at the same time. 

Indeed the renormalization group {RG) requires that every mass scale of the YM theory must 
depend on the canonical coupling constant, gYM, only through the RG invariant scale, Aym'- 

AyM = Aexp(-— ^)(/3og2^) W(i + ...) (1.1) 

whose dependence on the coupling constant is equivalent to the knowledge of the exact beta func- 
tion of the theory in some scheme. The dots refer to the higher loop contributions irrelevant in the 
ultraviolet, while the non-analytic two-loop result is explicitly displayed. Eq.(l.l) in turn implies 
that an amazing asymptotic accuracy, as gyM vanishes ^ when the cutoff. A, diverges, is needed to 
solve the mass gap problem and that the mass gap is zero to every order of perturbation theory. 

One possibility is that such finest asymptotic accuracy may be achieved only by an exact 
solution both on the ultraviolet side, for the beta function, and on the infrared side, for the mass 

gap- 
While an exact solution of the YM theory, to use just an euphemism, seems completely outside 
the reach of the present techniques, in this paper we propose an exact solution in the large-A'^ limit 
of the SU{N) YM theory for the beta function and for the glueball spectrum restricted to a special 
sector of the theory. 

1.1 A purely tield theoretical presentation: summing planar diagrams combining the Nicolai 
map with the holomorphic gauge 

It has been known for a long time that the large-A^ limit of the pure YM theory can be defined 
by the celebrated Makeenko-Migdal loop equation |§^, Eq.(10.8): 



/ 



, dxa<^Tr{^^^>{x,x;A))> 
L,, 2g^ OAa{x) 



i f dxa f dya5'''^\x-y) < Tr'¥{x,y;A) >< Tr'¥{y,x;A) > (1.2) 

'^ ^ V.V '^ AvA- 



Therefore the mass gap problem is not a strong coupling problem. 
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where: 

^>{x,y-A)=P&xpij Aadxa (1.3) 

is the Wilson loop, i.e. the holonomy of the gauge connection, Aa- The left hand side of the 
Makeenko-Migdal loop equation contains the critical equation for the YM action and the right 
hand side contains the contribution of the "change to loop variables" ^. It has been known for long 
that the Makeenko-Migdal loop equation is a compact way of "summing the planar diagrams" |Q]. 

However, the solution of the Makeenko-Migdal loop equation lives in a von Neumann algebra 
that is not hyperfinite, i.e. that is not the weak limit of a sequence of matrix algebras. Such non- 
hyperfinite algebra is too large a mathematical object to handle, sect.(1.3). 

Nevertheless, the algebra of local single trace operators at next to leading ^ order is, in a sense, 
the most simple as possible. Correlators of such operators, 0{x), are conjectured to be an infinite 
sum of free fields, saturating the logarithms of perturbation theory, sect.(12. 1): 



Goip") = je'P' < 0{x)OiO) >conn d 



\ 



= I^-2o(/)/"-^log(^) (1.4) 

The basic philosophy of this paper is to try to disentangle the more limited, but very interesting 
information on the spectral side, contained in the free propagators at next to leading order, from the 
overwhelming information, but less relevant from a spectral point of view, contained in the algebra 
of non-local Wilson loops at the leading ^ order. 

This disentanglement is obtained constructing new kinds of Wilson loops, called twistor Wil- 
son loops for geometrical reasons, built by a connection B^, defined in Eq.(4.28): 

^>{B;,;L„„)=Pexpi f {A, + XDu)dz + {A-, + X-^Da)dz (1.5) 

with the property that their v.e.v. is trivially 1 in the large-A'^ limit, but non-trivial at next to leading 
order. 

The twistor Wilson loops are defined initially in L'^(A^) YM theory on Euclidean space-time, 
R^ X R^Q, where the first factor is the ordinary Euclidean plane and the second factor is a non- 
commutative plane with non-commutative parameter, Q. The limit — ;• oo is the most relevant 
for us, since this limit is equivalent to the large-A'^ limit of the SU{N) YM theory on commutative 
space-time, sect.(4). 

The aim of this paper is to solve exactly a new loop equation in \wg&-N YM, that involves the 
new twistor Wilson loops. The basic strategy is to solve the new loop equation for twistor Wilson 
loops by a change of variables [^]. We set for simplicity B = Bi in the following. 

In fact despite the twistor Wilson loops satisfy a Makeenko-Migdal loop equation of standard 
type, Eq.(10.13): 



There is a way of deriving the Makeenko-Migdal loop equation as a change of variables from the gauge connection 
to the trace of its holonomy, i.e. the trace of a Wilson loop [0]. 



Glueballs in large-N YM by localization on critical points Marco Bochicchio 

if dz I dzS'^^^z-w) < Tr^>{w,z;B) >< Tr^>{z,w;B) > 

(1.6) 

the Makeenko-Migdal loop equation is as difficult to solve for twistor Wilson loops, ^{x,x;B), as 
it is for ordinary Wilson loops, *P(x,x;A). 

However, twistor Wilson loops satisfy a new loop equation, called holomorphic loop equation, 
obtained by a change of variables, actually two changes of variable. It is this holomorphic loop 
equation that is possible to solve exactly [^]. 

The first change of variables, worked out in sect.(3.3), is a non-supersymmetric {SUSY) ver- 
sion of the Nicolai map [H, Hi, previously known in .yV = 1 SUSY YM theory [P, 110] sect.(3.2), 
from the gauge connection, A„, to the anti-self dual (ASD) part of the curvature, n^g. This change 
of variables is not one-to-one everywhere in function space even after gauge fixing. We discuss 
later in this introduction how this technical difficulty is solved in this paper. 

There is a further change of variables to a holomorphic gauge, from n = n^y + //Iq3 to /i' |^], 
that has the scope of avoiding that the right hand side of the loop equation, i.e. the contribution 
of the "change to loop variables" be field dependent. Indeed this would make the reduction of the 
loop equation to a critical equation impossible. 

The holomorphic gauge, sect.(6), is the gauge in which B^ =Ai + Da = 0. This gauge can be 
reached by a gauge transformation in the complexification of the unitary gauge group. Therefore 
strictly speaking the choice of a holomorphic gauge defines a change of variables and not a proper 
gauge transformation of the YM theory. 

The result is the new holomorphic loop equation reported in sect.(6). It is obtained following 
the Makeenko-Migdal technique, as an identity that expresses the fact that the functional integral 
of a functional derivative vanishes: 

The new holomorphic loop equation for twistor loops follows: 

< Tr{^^^^{B';L,,)) >=- j -^ < Tr^>{B';L,„) >< Tr^>{B';L„,) > (1.8) 
oijL'{z,z) nJL--z-w 

where ^{B';L^^) is the holonomy of B in the gauge B'^ = 0. The Cauchy kernel arises as the kernel 
of the operator d^^ that occurs by functionally differentiating *P(B';L,,). 

Also the holomorphic loop equation contains in the right hand side the contribution of the "change 
to loop variables" as the Makeenko-Migdal loop equation does. In fact the right hand side contains 
the contribution of a certain subclass of planar diagrams that contribute to the twistor Wilson loops. 
The other planar diagrams being already included in the logarithm of the Jacobian of the change of 
variables that occurs in the "effective action", F, in the left hand side. 

However, a main point of this paper and of [^ is that the right hand side of the holomorphic 
loop equation vanishes provided the twistor Wilson loop has a backtracking cusp, i.e. cusp with 
cusp angle, D., at the point where the loop equation is evaluated. We follow the convention that a 
straight line has cusp angle n while a backtracking cusp has cusp angle 0. 
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The vanishing of the right hand side impUes the localization of the holomorphic loop equation, 
since for the cusped twistor Wilson loop the holomorphic loop equation reduces exactly to a critical 
equation for an effective action: 

We can choose the cusped loop with any number of backtracking cusps. In particular we can choose 
the standard Makeenko-Migdal loop with the shape of the symbol oo, but with a double backtracking 
cusp at the non-trivial self-intersection point and the equation of motion for the effective action 
inserted precisely at the cusp, sect.(6) and sect.(lO). 

We explain in this introduction why the localization of the holomorphic loop equation occurs. 

The twistor Wilson loops have the following fundamental properties at large N. 

They are non-trivial as operators. Indeed the curvature of the connection, B, that occurs in the 
twistor Wilson loops for A = 1 is: 

F{B) = ^ (1.10) 

where /i is a field of ASD type, Eq.(5.4). If the twistor Wilson loops are chosen in the adjoint 
representation, their v.e.v. factorizes in the large-N limit in the product of the fundamental and 
conjugate representation. The sector generated by the local operators that are polynomial in (/x,/i) 
in the fundamental representation plays a special role in this paper and it is referred to as the 
holomorphic/antiholomorphic sector. 

The twistor Wilson loops have trivial expectation value at large A'^, i.e. their v.e.v. is exactly 1 
at large 6, Eq(4.29): 

lim <-Tr^'{Bx;K„)> = 1 



< —Tr^^{B^;L„„) > = < —Tr^^>{Bi;L„„) > 



(1.11) 



and their v.e.v. is A-independent for any 6. This is shown in sect.(4), in particular in sect.(4.5). 

Therefore their shape can be deformed at will without changing their v.e.v.. In particular they 
can be deformed by adding an arc that backtracks, i.e. a path that is oriented in one direction and 
then it comes back along the same path in the opposite direction. This invariance property is known 
as zig-zag symmetry ^ and it is well known classically for any Wilson loop, but for twistor Wilson 
loops it holds also at quantum level. The reason is that in a regularized version the backtracking 
arc becomes the boundary of a tiny strip ending with a cusp. 

Now, the v.e.v. of an ordinary cusped Wilson loops has an extra divergent logarithmic con- 



tribution, with respect to smooth loops, known as cusp anomaly [12]. The cusp anomaly, i.e. the 



coefficient of the logarithm, actually diverges when the cusp backtracks, i.e. when the cusp angle 



D. tends to (see for example 1 13, 14]) 



< TrW{x,x;A) >~ exp(--log(ZA)) (1.12) 



The zig-zag symmetry plays an important role in a string ansatz for the solution of the Makeenko-Migdal loop 
equation due to Alexander Polyakov |fll|]. 
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where / is the length of the cusp and 5 a numerical factor. In the Makeenko-Migdal loop equation 
of standard type the right hand side develops extra divergences if a Wilson loop develops a back- 
tracking cusp. Indeed, performing the two contour integrations along the loop in the right hand 
side of the Makeenko-Migdal equation: 

f , N ^ , oSym ,^,, ,,, 

Jl,, 2g^ 5Aa[x) 

= i I dxa j dya5^'^\x-y) < Tr^'{x,x;A) >< Tr\ > (1.13) 



it is obtained [13]: 



, dxa < :^Tri^^Wix,x;A)) > 
~ KPA'+Y^ — -^(7r-a™,p)A2) < Tr^'ix,x;A) >< Trl > (1.14) 

cusp ^^^^'■cusp 

where / is the perimeter of the loop and ^cusp the cusp angle at a cusp. In the limit in which the 
cusp angle ^cusp reaches the cusp backtracks and the cusp contribution to the contact term of the 
Makeenko-Migdal loop equation is divergent. 

On the contrary, the right hand side of the holomorphic loop equation, after the gauge invariant 
regularization by analytical continuation to Minkowski described in sect.(6), vanishes at the back- 
tracking cusp. Indeed the right hand side of Eq.(1.8) is a contour integral of a d^^' first taken with 
one orientation and then taken with the opposite orientation: 

1 w+{s+^,p) 1 w+is-^,p) 



dw+{s)5{z+{scusp)-w+{s)) = - ^ '"'^ +- 

J 2\w+(Scusp)\ 2 



2\w+ {Scusp ) I 2 I IV+ {Scusp ) I 

= ^4=0 (U5, 

The delta function arises by analytical continuation to Minkowski, sect.(6). The light cone Minkowski 
coordinate is chosen to be parallel to the tangent to the cusp. Thus the holomorphic loop equation 
for a cusped loop reduces exactly to a critical equation for an effective action, provided the equa- 
tion of motion for the effective action, F, is restricted to the subalgebra generated by twistor Wilson 
loops: 

< Tr{nB;Ll,,)j^^^^nB;l^U)) >=0 (1.16) 

where L' and L" are the two petals in which the Makeenko-Migdal loop with the shape of oo is 
decomposed. All the information about the original "change to loop variables", that is contained in 
the planar diagrams that contribute to a cusped twistor Wilson loop, is encoded in the logarithm of 
the Jacobians that enter the effective action: the Jacobian of the Nicolai map from A to /x and the 
Jacobian to the holomorphic gauge from /i to }x' . 

Thus an effective action, F, that represents the "sum of planar diagrams" for cusped Wilson 
loops has been obtained. F contains the interesting information of the localization of the loop equa- 
tion, since we can compute in principle the glueball spectrum, in the holomorphic/antiholomorphic 
sector generated by (/x,/i), by fluctuations around the critical points of F. 

8 
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However, there are some technical difficulties. To interpret the preceding equation strongly ^: 

(1.17) 



5ju(z+,z- 

we must realize explicitly the restriction to the subalgebra of twistor Wilson loops and in addition 
we must actually find the critical points. 

Afterward we must be able to check that the fluctuations around the critical points are actually 
suppressed by ^. 

Finally, we must be able to understand the loci in function space where the Nicolai map is not 
one-to-one. 

Surprisingly, the three problems are actually deeply linked, in such a way that the solution of 
the third one furnishes the solution of the other two. 

Indeed the first change of variables, the one from the connection A^, to the ASD field, /i^g, 
may not be one-to-one even after gauge fixing. 

Firstly, we work out the supersymmetric case in sect.(3.2), just as a simpler exercise. It turns 
out that the Nicolai map is not one-to-one and at places where multiple solutions of A^ occur for a 
given n^n, there are moduli of Aq. and correspondingly zero modes in the Jacobian of the change 
of variables. However, in the SUSY case and for a special observable, the gluino condensate, the 
only relevant zero modes arise by instantons. But these zero modes are essential to reproduce the 



exact supersymmetric beta function of Novikov-Shifman-Vainstein-Zacharov (NSVZ) [15]. 

Secondly, we work out the non-supersymmetric case, in sect.(3.3). It turns out that, as in the 
SUSY case, even the first coefficient of the pure YM beta function cannot be reproduced from F 
unless zero modes occur, sect.(3.4). But unlike the SUSY case localization on instantons does not 
reproduce the second coefficient of the YM beta function, sect.(3.4). Thus in the pure YM case 
we need to understand the moduli of the Nicolai map generically in function space. This is more 
difficult, since it is known that there is no Hausdorff, i.e. separable, moduli space of all bundles 
even in two dimensions. Therefore we have to introduce a dense set, in the sense of distributions, in 
function space, on which the moduli problem has a complete answer, by the standards of differential 
geometry of fiber bundles, sect.(7). On this dense set also the second change of variables, from /i 
to n' is well understood, sect.(7). 

In fact this leads to understanding zero modes also in the pure \aige-N YM case and to the 
large-A^ exact beta function of sect.(9) and sect.(l 1), via the effective action obtained by the cusped 
holomorphic loop equation. 

The aforementioned dense set involves interpreting the non-SUSY Nicolai map as defining a 
hyper-Kahler reduction on singular instantons, with singularities of magnetic type: 



ap-L^apiP)S'''K^-^p) (1-18) 



P 



The restriction of the YM measure to the lattice hyper-Kahler reduction, that is dictated by the need 
of reproducing the beta function at the critical points of F, has two fundamental consequences, 
sect.(7). It realizes in an explicit and mathematically well defined way the algebra of twistor Wilson 



If all the matrix elements of an operator vanish, the operator actually vanishes. 
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loops, via the equivalence of the hyper-Kahler quotient to the holomorphic quotient. This solves 
the first problem too. 

In addition the lattice hyper-Kahler reduction of the YM theory in the Nicolai variables is 
locally abelian, up to zero modes, because the triple, i^^Ap), is commutative at any given lattice 
point, p, because of fundamental properties of the Hitchin equations, sect.(7). Thus there is a 
gauge, referred to as the singular gauge in this paper, in which fluctuations of ix~q [p) reduce to 
fluctuations of eigenvalues, that are suppressed in the large-A^ limit, sect.(12). This solves the 
second problem. All the remaining non-abelian degrees of freedom of the theory on the lattice 
hyper-Kahler quotient are the zero modes associated to the moduli of the local systems determined 
by the hyper-Kahler reduction. Therefore the solution of the third problem furnishes the solution 
to the other two problems as anticipated. 

Finally, the critical points of the effective action can be explicitly found, using the A -independence 
of the v.e.v. of the twistor Wilson loops, Eq(l.l 1), as fixed points for the action on the Nicolai vari- 
ables of the semigroup that rescales A, sect.(8). They are a "lattice of surface operators" with 
Zn holonomy labelled by a magnetic quantum number, k, sect.(7). The surface operators with 
Zi\i holonomy condense because of asymptotic freedom (AF), sect.(12). In turn this leads to the 
physical interpretation of the magnetic condensate of surface operators in terms of 't Hooft elec- 



tric/magnetic duality [|1^, |17[ |18p that requires that, if the YM theory has a mass gap, then either the 
electric charge condenses (Higgs phase) or the magnetic charge condenses (confinement phase). 

The Jacobian of the non-SUSY Nicolai map gives rise at the leading large-A^ order to an anoma- 
lous dimension, sect.(3.4), that contributes to the one-loop exact Wilsonian beta function, sect.(9), 
and to a canonical beta function of NSVZ type, once the effective action is restricted to surface 
operators, sect.(ll). At next to leading order this Jacobian gives rise to the glueball kinetic term, 
sect.(12). 

Instead the logarithm of the Jacobian of the holomorphic gauge is precisely the glueball poten- 
tial at leading order, whose second derivative is the glueball mass matrix at next to leading order. 
Both are computed in sect.(12). A detailed study shows that a mass gap, proportional to the den- 
sity, pic, of surface operators of magnetic quantum number, k, arises, sect.(12). The density scales 
as {kN)^'i times the square of the /?G-invariant scale, but the laige-kN vanishing of the density is 
compensated by large-^A^ degeneracy factors in the mass matrix and by large-A^ degeneracy factors 
in the kinetic term, in such a way that the mass gap survives the laige-N limit. 

Actually the spectrum is of pure poles, because the non-local terms in the effective actions are 
suppressed by higher powers of the density of surface operators in the expansion of the effective 
action in powers of p^, that occurs as the usual expansion of a functional determinant in terms of 
one-loop graphs with multiple insertions of the background field. 

In fact the spectrum of fluctuations of composite surface operators forms a trajectory linear 
in k that does not include any massless state. While the spectrum does not depend on the choice 
of a particular observable of the YM theory, in our case surface operators, in [ [T9| ] we describe 
the correspondence between the two-point correlators of long composite surface operators and the 
glueball propagators. One ingredient of this correspondence is the coincidence of the anomalous 
dimensions of long surface operators with the anomalous dimensions of the ground state of the 



Hamiltonian spin chain in the thermodynamic limit [20, 21], that provides anomalous dimensions 



of local operators in the one-loop ASD integrable sector of large-A/^ YM, sect.(12). 

10 
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Introducing a lattice of surface operators allows us to write a lattice holomorphic loop equation. 
There is indeed a lattice version of holomorphic loop equation, Eq.(lO.l), in which the lattice 
points are the locations of singularities of surface operators. This lattice version arises through the 
discussion presented in sect.(7): 

< '^''^airrrTT^'^^^''^^)) >= ^ / -^<Tr^'iL...)><Tr^"{K,^)> (1.19) 
o/i [Zp ,Zp) n Jl,,,,„ Zp-w 



'cpcp 



The lattice holomorphic loop equation equation allows to build more on the geometrical side of the 
localization, discussed more extensively in sect.(1.3). 

By homological localization of the loop equation we mean a deformation of the loop that is 
trivial in homology and for which the term that arises as the "change to loop variables" vanishes, in 
such a way that the loop equation is reduced to a critical equation for an effective action [^]. Hence 
the needed homological deformation has to satisfy the following properties. It has to be trivial in 
homology. It has to leave the expectation value of the loop invariant. It has to imply the vanishing 
of the quantum term in the loop equation, i.e. of the term that contains the contour integral along 
the loop, sect.(lO). 

Now we deform the twistor loops by adding a backtracking arc ending with a cusp at each 
lattice point. The lattice points associated to the divisor of surface operators become the cusps that 
are the end points, p, of the backtracking strings, bp, that perform the deformation of the loop, 
L. Adding the backtracking strings implies that the contribution of the "change to loop variables" 
vanishes for the modified loop: 



Si^'{zp,Zp) 



< M s .u[ J, 'i''{LU [bp])) >= (1.20) 



This phenomenon is called homological localization of the holomorphic loop equation because, in 
geometrical language, is dual to the cohomological localization that is described in sect. (1.3) and 
in sect.(3.2) for the gluino condensate of the SUSY YM theory. 

In the cohomological localization of sect.( 1.3) functional integrals are reduced to critical points 
by adding "a total differential in function space". 

In the homological localization of sect.(lO) loop equations are reduced to critical points adding 
to the loop "vanishing boundaries", i.e. the aforementioned backtracking arcs ending with cusps. 

The duality is no mystery. It is just the Stokes theorem in the form: 



d(0= (0 (1.21) 

V JdV 

that relates differentials, i.e. the cohomological side, to boundaries, i.e. the homological side. 

In the holomorphic loop equation of the YM theory the boundaries are the backtracking arcs 
that are added to the loop to get localization of the holomorphic loop equation. In the SUSY YM 
theories the coboundary is the differential furnished by the (twisted) supersymmetric charge. This 
is not a physical duality. It is purely a mathematical duality of the two different localizations. 

1.2 Summary of results 

From a purely computational point of view we may present our results in terms of the glueball 
propagators in a certain sector of the large-A'^ YM theory. The glueball propagators that we refer to 
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are initially defined in Euclidean signature, constructed by means of fluctuations, 5/i ^, of certain 
surface operators ^ supported on a Lagrangian submanifold of four-dimensional space-time ^: 

A%J < -^Try{^ll){z,z,z,z)-^Tr^{m{0,0,0,0) >conne*^''^''^d'z 

(1.22) 

where (z = XQ + ix\,z = xq — ix\ ,u = X2 + ix3,u = X2 — 1x3 ) are complex coordinates in Euclidean 
R^ X R^, (f-z = dx^dxy and 

i" = 2('"oi-«M(B) (1-23) 

with ji the Hermitian conjugate. The following identification holds up to (infinite) factors ^: 

I ^i{z,z,z,z)e'^P^'+P^'^d^z ~ ^{F,- - iF^3){p„p-„p„Pz) (1-24) 



where 



(^01 - iFo3)iPz,P-z,Pu,Pu) = J{F,, - iF^,)iz,z,u,a)e'^P-^~'+P-^'+P"'+P^'"Uhd\ (1.25) 
is the Fourier transform with 

Fap=Fap-Fap 

1 

Fap = 2^aPysFaji (1-26) 

the anti-selfdual (ASD) part of the curvature of the gauge connection. The identification holds 
because the fluctuating field of surface operators, 5/i, has non- vanishing momenta dual in the 
Fourier sense to the Lagrangian support of surface operators and it occurs in the effective action 
with zero momenta dual in the Fourier sense to the manifold normal to the support. ^ = NN is 
the total rank of the gauge group. The first factor of N is the rank of an SU (N) gauge bundle that is 
embedded by non-commutative Morita equivalence into U(N x N). The construction is explained 
in sect.(7) and sect.(12). After analytically continuing to Minkowski space-time in its simplest 
form of our result reads ^: 

A^J < —Tryr{nll){x+,x-,x+,x-) — Tryr{nll){0,0,0,0) >conne'^P^'-+P-'+Ux+dx- 



^The precise definition is in Eq.(12.39). 
Tlie idea of integrating in YM theory on local systems associated to an arbitrary parabolic divisor appeared for 

26l El in the YM 



the first time long ago in some papers by us pj, E3], by embedding the Hitchin fibration p4| , pq, 

functional integral, and physically corresponds to integrating over surface operators ante litteram \ki ] . Explaining how 

surface operators arise in large-A' YM is in fact the subject of this paper. 

^The plane (z, z,z,z) is Lagrangian for the symplectic form dzAdz — du/\ du in four dimensions. 

**The identification extends to composite operators in a certain asymptotic sense, see below and sect. (12). The 
infinite factors are actually defined and regularized in sect.(12). 

The symbol --^ stays for "equal up to constant irrelevant numerical factors" or "equal up to irrelevant additive terms" 
depending on the framework. 
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1 - kp^^Kl 



-^^ lt\ -OC'p+P- + kAl, 
1 Y^ k A(^ 



-^^ i=i -a'p+P- + kKlf 

(1.27) 

where, (x+ = X4 + xi ,x_ = X4 — X1) are light-cone coordinates, (;?+ = pi\ + p\^p^ = p4 — pi) are 
light-cone momenta, and A^ is the renormalization group invariant scale in the Wilsonian scheme. 
Pfe is the density, in units of A|,, 

p=£5P)(z-z^) (1.28) 

p 

of surface operators carrying at each lattice point, p, magnetic charge k and holonomy valued in 
the center, Z^, of the gauge group, i.e. such that: 

e^'^p=e'— (1.29) 

with: 

\{F,;,-iF,;,)=Y^^p5^^\z-Zp{u,u)) (1.30) 



"^ 2 
and Zp{u, u) = Zp- Pk scales with k as: 



Pk^l (1-31) 



and the dimensionless inverse "string tension" in units of A;^^ is: 

a' = —n (1.32) 

371 

where n > 2 is a finite positive integer that depends on the choice of the renormalization scheme. In 
fact n can be reabsorbed into a redefinition of Ai^ in any of the countably many possible inductive 
sequences that define the SU (00) group as a limit of a sequence of finite dimensional SU {N) , sect.(7) 
and sect.(12). 

The peculiar support, (;c+ , ;c_ , ;c+ , x_ ) , of the correlator arises as the projection with Minkowski 
signature on the base of a Lagrangian submanifold of the twistor space of (complexified) Euclidean 
space-time that occurs in our approach. The field /x is dimensionless and normalized in such a 
way that the correlator in Eq.(1.7) be renormalization group invariant ^^. The mass spectrum, 
—a'p+p- +kA^ = 0, and the multiplicity, k, that occurs in the numerator of Eq.(1.7) in the second 



'in the canonical normalization also an extra factor of the anomalous dimension occurs. This is reported in the 
formulae below and discussed in sect.(12). 
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line, are in fact exact ^Mn a certain asymptotic expansion in powers ofN^'^ of the gauge connection 
of the surface operators that occur in our approach ^^. 

At large N there is a Wilsonian scheme in which the Wilsonian beta function is one-loop exact 
[§] and a canonical scheme in which the beta function has a form NSVZ that reproduces the first 
two universal perturbative coefficients [§] '^: 

^'"^ -Pogi (1.33) 



<9logA 
and 

,3 J 1 „3dlogZ 



dg _ -Pog + (4iF^ ^toJA ,, , .. 



with: 

1 11 



5 log A 



(471)2 3 

(1.35) 

where g = gyM^ i^ ^^^ '^ Hooft canonical coupling constant and ^j°f^ is computed to all orders in 
the 't Hooft Wilsonian coupling constant, gw, by: 

-,, ^ 1 10 „2 

5logA l+cgl, 

with c a scheme dependent arbitrary constant. Indeed since ^^"^^ to the lowest order in the canon- 
ical coupling is: 

aiogZ _ 1 10 2 



5 log A (471)2 3 



+ ... (1.37) 



the correct value of the first and second perturbative coefficients of the beta function [ plj [32i p3| , |34| ] 
arise: 

dg o 3 / 1 1 10 o 4 , . 



a log A ^"* '(471)2(471)2 3 '^"(471)2- 



1 11 , 1 10 44, c 



(471)2 3 " (471)4 ' 3 3 



1 11 , 1 34 



+ ... (1.38) 



(471)2 3 * (4;j;)4 3 

' ' Large-A' exact linearity of the spectrum, rather than only asymptotic linearity, may look surprising. However, the 
ratio between the masses of the two lowest scalar states in pure S(/(8) YM has been found numerically to be compatible 
with the value Vl M]. We would like to thank Michael Teper for a clarifying discussion about this point at the Galileo 
Galilei Institute workshop on "Large-A' Gauge Theories" (201 1), hereafter referred to as the GGI workshop. Indeed the 
best fits in [M] for the continuum limit ratios of the masses (in units of the lattice i?G-invariant scale. Table (7.14)) of 
the 7^'' glueball, nio++« = 4.71(29) and mQ-+ = 4.72(32), to the mass of the lowest scalar, m()++ = 3.32(15) agree with 
very good accuracy with \/2. Yet, there are larger statistical errors than this agreement may suggest. There is also a 
m2++ = 4.65(19) glueball, essentially degenerate in mass with the two aforementioned scalar states, that may in principle 
couple to our operators. Nevertheless, for reasons explained in the paper, we suggest that the spectrum in Eq.(1.7) is all 
made by scalars. 

'^The precise meaning of this statement is clarified in sect.(12). 

'^It has been known for some time that the Wilsonian, gw, and the canonical, g, coupling constant have different 
beta functions in general yS(M. 
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In fact a whole family of correlators of the Fourier transform of composite operators of naive 
dimension 4L, 0^{p+,p-), constructed by surface operators supported on the aforementioned La- 
grangian submanifold, are computed in sect.(12). In the Wilsonian scheme the result reads: 



~ A 2 /,2(2L-1) a4(2L-1) 



w 



{-p+p^f^-Hog. 



A2 



(1.39) 



where the dots stand for contact terms, i.e. distributions whose inverse Fourier transform is sup- 
ported at coinciding points. 

In the canonical scheme for the same objects anomalous dimensions arise: 

< Tr^rO^{p+,p^)n^rO^{-p+,-p-) >'£Ln 

= g\-p+P-)Z-'-^ {-P+P-) < Tr,j.O^{p+,p^)Tr,_^0^{-p+,-p^) >S„ 

~/(-p+P-)Z"'^(-p+P-)(-p+P-)'^"'log^^ (1.40) 

where g{—p+P- ) and Z{—p^p^ ) are the {RG improved) momentum dependent canonical coupling 

in Eq.(1.14) and renormalization factor in Eq.(1.16). A deeper and more complete discussion of 

the relation between correlators of surface operators and glueball propagators can be found in [jl^]. 

For large L the anomalous dimensions agree ^^ with the anomalous dimensions of the ground 



state pOj , |21[ ] of the Hamiltonian spin chain in the thermodynamic limit, furnishing an identifica- 



tion, that is at least asymptotic for large L, between composite surface operators ^^ and composite 



local operators of YM in some regularization scheme QIQQ. Indeed the Hamiltonian spin chain is 
an integrable model by which the anomalous dimensions of composite operators of large-A'^ YM in 
the ASD and SD sector can be computed exactly at one-loop. 

The exactness of our formula for the spectrum is not affected by the possibly only large-L 
asymptotic identification of the operators, since the spectrum depends only on the occurrence of 
poles in any correlator of gauge invariant operators. However, the aforementioned asymptotic 
identification suggests that all the glueballs in our spectrum are in fact scalar, since this is so for the 
operators that correspond to the ground state of the Hamiltonian spin chain in the thermodynamic 



limit [21]. 



'^The agreement is at one loop since anomalous dimensions are universal, i.e. scheme independent, only at one loop. 
Actually they agree also for L = 1 , since in this case the anomalous dimension is determined by the beta function via the 
factor oi g . 
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We actually mean that we identify the Fourier transform of our composite surface operators with the Fourier 

m + 



transform of composite local operators in the same fashion as in Eq.(1.4). For example, 0^{p+,p-) --^ \{Ff 



'Foj,)\^^{p+,P-,P+,P-)- 
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Yet, there is an extension of our approach to fluctuations of surface operators defined by con- 
nections with wild singularities (i.e. pole singularities of any order), that are naturally associated 
to Regge trajectories of higher spin (sect.(13)). Their contribution is not computed in this paper. 

Now we re-explain the basic ideas underlying our computations in geometric terms as follows. 

1.3 A geometric point of view: localization by cohomology and by homology 

In the last thirty years we witnessed the geometrization of theoretical high energy physics. 
This geometrization has several faces but the one that we refer to consists in computing exactly 
functional integrals by geometrical methods. The key idea in solving the analytical problem of 



performing an integral by geometrical methods lies in the work of Duistermaat and Heckman [ |35| ] 
on exact localization of the integral of the exponential of the Hamiltonian of a torus action on a 
compact symplectic manifold on the fixed points of the torus action '^: 

1 /■ „ . -., ^.271 



f (o"exp{-eH)=Y,{ — rPr\o)p'd^Hp) (1.41) 

J p c 



where Pf is the Pfaffian of the skew matrix, (Op ^d^Hp, at the fixed points, P, and d^Hp the Hessian 
of the Hamiltonian at P. Such localization has the following cohomological nature according to 



Atiyah and Bott [37] and Bismut [38, 39 1. In finite dimension the integral of the exponential of a 
closed form ft), dco = 0, on a compact manifold without boundary, M, defines a cohomology class 
invariant for the addition to the closed form of an exact differential, da, (i.e. of a coboundary), 
since da is trivially closed because d^ = 0. Rescaling the exact differential by a large factor, under 
suitable positivity assumptions, the integral of the cohomology class gets localized on the critical 
points of the exact differential and the saddle-point approximation turns out to be exact: 

exp{-(o-tda)= / exp(-ft)) (1-42) 

M Jm 

Indeed the first f -derivative vanishes because it is the integral of a coboundary and Jm^^ — ^■ 

— exp{-(o-tda) = - d{aexp{-(o-tda)) = (1-43) 

dt Jm Jm 

Therefore the integral is f-independent and in the limit ? — ;■ oo can be evaluated by the saddle-point 



method. Bismut [^, ^ was the first one to extend rigorously this kind of argument to infinite 
dimensions, actually to a functional integral in one dimension, i.e. to quantum mechanics. In the 
quantum mechanical setting there are (essentially) no existence problems for functional integrals 
and the localization argument is in fact a mathematical proof. 

However, the subject blossomed in quantum field theory only after Witten paper on localiza- 
tion in two-dimensional Yang-Mills theory [ p^ and Witten work on Donaldson invariants [^] that 
introduced localization in four-dimensional supersymmetric gauge theories, by identifying the dif- 
ferential needed to define the cohomology with a twisted super-charge, Q, satisfying Q^ = 0. In 
turn Witten twist of supersymmetry requires to start with at least an ^ = 2 SUSY YM theory. 

Thus the infinite-dimensional field theoretical analog of Eq.(1.42) is: 

Joexp{-SsusY-tQa) (1.44) 



16 
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with: 

QO = 
QSsusY=0 (1.45) 

There have been a number of apphcations of the locahzation idea in four-dimensional SUSY gauge 
theories, among which we mention the Nekrasov computation [ p^ ] of the prepotential in ^ = 2 



SUSY gauge theories, that reproduces by localization methods the Seiberg-Witten solution [ ]43P for 
the same object, and Pestun [ p] ] computation of certain twist-S't/^'F invariant Wilson loops in ^ 
= 4,2,2* SUSY gauge theories. 

From a purely mathematical point of view these exact results state the equality between a 
mathematically not well defined object, the original functional integral, and a mathematically well 
defined and explicit answer, the result of the localization. 

However, from the point of view of theoretical physics, these results are in fact satisfactory 
since, waiting for a realization of the constructive program of quantum field theory in four dimen- 



sions [[45|], the explicit answer that is found by localization defines the functional integral by the 
rules by which it is computed and contributes to fix the properties that the yet-to-come mathemati- 
cal construction of the functional integral has to satisfy: the localization property indeed. 

The aim of this paper is to add, rather surprisingly, a non-supersymmetric chapter to the afore- 
mentioned exact results. 

The simplest way to present our basic result is to compare it with Nekrasov computation of 
the prepotential. In the first part of Nekrasov computation the functional integral that evaluates the 
cohomology of 1 (i.e. the partition function) is reduced by cohomological localization to a sum of 
finite dimensional integrals over the instantons moduli spaces: 



: lim f\exp{-Ssus¥-tQa) = Yexpi- ^^^ ^^ )A^kN f ^^ ^^ ^^^ 



This depends on the supersymmetry and has no analog in the pure YM case. On the contrary, in 
the second part of Nekrasov computation, the finite dimensional integrals over instantons mod- 
uli are reduced to a sum over the fixed points for the action of the torus L'^(l)^^' x U{1) x U{1) 
in SU{N) X 0(4) by applying the Duistermaat-Heckman formula, after a suitable ultraviolet and 
infrared regularization of the moduli space, by means of a non-commutative deformation parame- 
terized by e: 

Here SU{N) is the (global) gauge group at infinity and 0(4) the group of Euclidean rotations. 
These groups are symmetry groups also of the pure YM theory. The result of the localization can 
be resummed into an exact formula for the prepotential [p3|], ^, that is a function of the quantum 
moduli of the theory, that are related to the v.e.v. of the eigenvalues of the complex scalar field in 
the adjoint representation of the ^ = 2 YM theory: 

Z = exp(^^) (1.48) 
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Despite the prepotential is obtained by the locaUzation of a trivial observable, the cohomology of 
1, from a physical point of view it contains the interesting information of the localization. Indeed, 



according to Seiberg-Witten ||43|], the prepotential contains exact highly non-trivial quantum infor- 
mation. It determines an exact beta function and the low energy effective action in the Coulomb 
branch of the theory as a function of the translational invariant condensate of the eigenvalues of the 



scalar field '^. Thus the prepotential is used to reach conclusions about the physical theory [43] 
that by far exceed the very limited framework of its derivation by localization. 

A general feature of cohomological localization is that the saddle-point computation can be 
employed only for the specific observables that satisfy Eq.(1.45). A fortiori in pure YM, that has 
no SUSY, there is no hope that localization may hold, if any, but for very special observables. 

To say it in a nutshell, our basic idea for pure YM is to construct special trivial observables, 
called twistor Wilson loops for geometrical reasons, since they are supported on a Lagrangian 
submanifold of twistor space of complexified Euclidean space-time. In a technical sense the trivial 
twistor Wilson loops are in the homology of 1, rather than in the cohomology of 1, since in pure 
YM there is no SUSY and thus no interesting cohomology '^. The loop equation for twistor Wilson 
loops can be solved, since they are trivial, in the sense that it can be reduced to a critical equation 
for an effective action, i.e. it can be localized. Despite the effective action is obtained by trivial 
observables, it carries highly non-trivial quantum information, that exceeds by far the framework 
of localization of 1 . 

The effective action determines an exact large-A^ beta function and turns out to be a function 
of the density, pk, in units of A^, of the condensate of surface operators of magnetic charge k that 
occur in the localization of the twistor Wilson loops. 

In addition the effective action restricted to fluctuations of surface operators supported on a 
Lagrangian submanifold with Minkowski signature, obtained by a certain Wick rotation from the 
Lagrangian submanifold which the twistor Wilson loops are supported on, determines the glueball 
spectrum. 

The analytic continuation to Minkowski space-time is the only way to regularize gauge invari- 
antly the holomorphic loop equation for the twistor Wilson loops, that in turn leads to localization 
on the critical points of the effective action. 

We describe now in more detail what the twistor Wilson loops are. 

They compute the holonomies along loops of a modified non-Hermitian YM connection, the 
twistor connection. Its curvature is a non-Hermitian linear combination of the ASD part of the cur- 
vature of the ordinary gauge connection. These loops are supported on a Lagrangian submanifold 
in twistor space of complexified Euclidean space-time, locally the product of a two-dimensional 
surface immersed in complexified space-time with local (complex) coordinates (z,z) and of a one- 
dimensional curve immersed in the fiber of the twistor fibration with (not necessarily real) coordi- 
nate A. 

The twistor Wilson loops are chosen in the adjoint representation. The operator definition of 
twistor Wilson loops involves the parameter A, but their vacuum expectation value (v.e.v.) is A- 



'^In the Coulomb branch the eigenvalues are generically all different in such a way that the unbroken gauge group is 

'**There is in fact the Becchi-Rouet-Stora (BRS) cohomology associated to gauge-fixing, that leads to localization on 
gauge-fixed slices of gauge orbits, but it is not relevant for our purposes. 
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independent, in fact trivially 1 at large-A'^. Hence there is a non-compact real version U{1)r, of the 
complexification U{l)c, of one of the aforementioned U{1) '^ that acts by rescaling A in such a 
way that the v.e.v. of twistor Wilson loops is invariant under the aforementioned action. 

As a consequence we show that twistor Wilson loops in pure YM are localized on the sheaves, 
defined by the change of variables from the gauge connection to the ASD part of its curvature in 
the functional integral, fixed by the action ofU{l)R. In addition we show that there is a dense set 
in function space in a neighborhood of the fixed sheaves ^", that at \aige-N is classified by local 
systems on a sphere with a very large number of punctures and with fixed conjugacy class of the 
holonomy of the twistor connection around the lattice of punctures, with values in the complexifi- 
cation, SU{N)c, of the gauge group, modulo the global action oi SU{N)c- 

In the physics terminology the local systems are lattices of surface operators satisfying the 
self-duality (SD) equations with singularities: 

Fap=L^apiP)S^'^(^-^p) (1-49) 

P 

In addition the v.e.v. of the aforementioned twistor Wilson loops in the adjoint representation 
factorizes in the large-A'^ limit in the product of the v.e.v. in the fundamental and conjugate repre- 
sentation. Then to each factor the following argument applies. 

On the dense set described by local systems, by translational invariance we can assume that 
all the conjugacy classes of the holonomies are a copy of the same adjoint orbit, and that the orbit 
for a holonomy around one arbitrarily chosen point can be put by the global action of SU{N)c in 
canonical form, i.e. either in diagonal or in Jordan form. Now the global compact SU{N) gauge 
group acts on such diagonal or Jordan holonomy by conjugation. 

If the global gauge group is unbroken, as it is believed to be the case for pure YM, only the 
holonomies that are fixed by the entire SU (N) may occur at large-A^. Thus these holonomies at a 
preferred point are in fact valued in the center of the gauge group and their orbits reduce to points. 
But then by translational invariance all the orbits reduce to the center ^' . 

Besides we show that there is a homological explanation for this localization on fixed points 
based on a new localization theory of the loop equation for twistor Wilson loops, such that the actual 
fixed points that contribute to the twistor Wilson loops are the critical points of a certain effective 
action determined by the loop equation [^. The localization by homology of the loop equation, i.e. 
its reduction to a critical equation, is obtained deforming the loop by adding vanishing boundaries 
that are backtracking arcs ending with the cusps of the local system ^, an operation allowed by 
the large-A^ triviality of twistor Wilson loops, by dualizing the idea of deforming a closed form by 
a coboundary in the cohomological interpretation of the Duistermaat-Heckman localization. 



It is a U{1) in the Cartan subgroup of 0(4). 

'We refer to the support of the fixed measure as fixed sheaves to imply not any manifold structure for such a locus. 
However, for a dense set in function space the fixed sheaves at large-A' are in fact a manifold that is parameterized by 
the disjoint union of moduli of local systems with fixed conjugacy class of the holonomy of the twistor connection. We 
refer to fixed points instead when the fixed locus has no moduli and it is a set of disconnected points. 

The same conclusion is reached by an inductive argument on the holonomies around each point, without assuming 
translational invariance, since once the holonomy around a point is shown to be in the center by the assumption of 
unbroken gauge group, the global SU{N)c still acts on the holonomies around each of the remaining points. 
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In order to get localization, the first main technical innovation of our approach is a reformu- 
lation of the YM theory in terms of a change of variables that in the ^ = 1 SUSY YM theory 
has been known as the Nicolai map [^ ||]. The Nicolai map in the ^ = 1 SUSY YM theory was 



worked out by De Alfaro, Fubini, Furlan and Veneziano Vq, [lOp as a change of variables from the 
gauge connection to the anti-selfdual (ASD) part of the gauge curvature, that needs a gauge fixing 
to be locally invertible, with the property that the Jacobian of the map cancels precisely the fermion 
determinant in the light-cone gauge. 

As a preparatory exercise, the Nicolai map allows us to introduce localization also in the pure 
^ = 1 SUSY YM theory ^^ by means of the tautological Parisi-Sourlas supersymmetry associ- 



ated to the cancellation of the Jacobian with the fermion determinant [ |48| , [49| ]. While it has been 
known for some time that the Nicolai map can be associated to cohomological localization ^^, the 
localization by the Nicolai map has never been worked out in asymptotically free gauge theories 
because of the following difficulty. Naively the Nicolai map maps ^ = 1,2 SUSY YM into a free 
theory, that cannot hold true hterally. The question arises for example how to reproduce the NSVZ 
beta function [15, ^ by means of the Nicolai map. Our simple but key observation is that the 
cancellation of determinants occurs only up to zero modes. Therefore the divergences associated 
to the Pauli-Villars regulator of the zero modes occur. 

In fact understanding how the NSVZ beta function occurs by cohomological localization via 
the Nicolai map in this paper is only an exercise for understanding localization of the aforemen- 
tioned twistor Wilson loops in large-A'^ pure YM in the ASD variables. The crucial point is that the 
localization on the fixed points of the U{1)r action can be obtained only in the ASD variables. 

The second main technical innovation consists in interpreting our non-^'L'^^'F version of the 



Nicolai map in the pure YM theory [|6|, ^, 23] as hyper- Kahler reduction ||5^, ^ on a dense set in 
function space, that corresponds to a lattice of surface operators in the physics terminology. This 
is an analytical and differential geometric construction that does not need any supersymmetry. It 
reduces the YM functional integral to a finite dimensional integral with respect to a product measure 
on a lattice and it is the analog of the first part of Nekrasov computation in the supersymmetric case. 
The physics interpretation is that the localization of the twistor Wilson loops in the laige-N YM 
theory is described in terms of variables that are of purely magnetic type, realizing, in the technical 
sense of localization of twistor Wilson loops, a new version of 't Hooft long-standing ideas ^^ on 
the YM vacuum as a dual superconductor [16, I^, 18]. 



In particular 't Hooft duality in YM theories with fields in the adjoint representation requires 
that, if the theory has a mass gap, then either the Zf^ magnetic charges condense (confining phase) or 
the Z]\i electric charges condense (Higgs phase). Localization by homology of twistor Wilson loops 



^^Witten already observed in his paper WU] on Donaldson invariants that although most naturally formulated in the 
^ = 2 SUSY theory localization could be extended to certain theories with only ,yK = 1 SUSY, called ^ = 2* theories 
that involve anyway the occurrence of a scalar field. In pure -.A^ = 1 SUSY gauge theory Witten localization does not 
apply directly since there are no scalars. It is always possible to give the scalars of the ^ = 2* SUSY theories large 
masses in order to obtain at low energy pure ._/f = 1 SUSY YM theory. This leads to the modern "weak coupling" 
approach to the computation of the gluino condensate |Hq,kj]. 

In the lectures [pQ] it is worked out the zero dimensional case of the Nicolai map and it is shown indeed that 
coincides with localization. 

In addition to 't Hooft original papers [Ma [O, [la] see also ||54|] for a very neat account of 't Hooft duality. 
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in pure YM realizes the first alternative, in which the electric charge is unbroken ^^ and the magnetic 
charge is broken in superselection sectors labelled by k, the magnetic charge at a (lattice) point, that 
are degenerate ^^ for the laige-N renormalized effective action that occurs in the holomorphic loop 
equation. From the localized renormalized effective action restricted to fluctuations supported on 
the aforementioned Lagrangian submanifold it follows also that in each sector there is a mass gap 
proportional to Vk in units of the common /?G-invariant scale. 

Yet, we should stress that locaUzation is by no means a universal concept, but it applies only 
to special observables. Therefore, it would be completely wrong to employ surface operators oiZjM 
holonomy to compute general observables of the YM theory, as it would be completely wrong to 
employ instantons to compute anything but the gluino condensate in ^/K = 1 SUSY YM. 

The mathematics interpretation is that we are in fact representing the YM functional integral 
as an adeUc integral over (the moduli space of) local systems. 

The two aforementioned technical innovations are crucial for our twofold approach to local- 
ization in pure YM. 

Firstly, as we just explained, following the spirit of the Duistermaat-Heckman idea our new 
kind of localization in the \aige.-N pure YM theory involves the action of a semigroup fixing the 
v.e.v. of twistor Wilson loops and contracting the support of the functional YM measure, resolved 
into ASD orbits by our non-supersymmetric version of the Nicolai map. At technical level the 
hyper-Kahler reduction to surface operators furnishes a structure theory of the locus of the fixed- 
points. 

Secondly, as well as the Duistermaat-Heckman localization on fixed points has a cohomo- 
logical explanation, so the new localization on fixed points in pure YM theory has a homological 
explanation. Indeed there exits a new holomorphic loop equation for twistor Wilson loops that 
can be localized, i.e. reduced to a critical equation, by deformations of the loop that are vanish- 
ing boundaries (backtracking arcs) in homology, in the dual sense to which a cohomology class 
represented by an integral of the exponential of a closed form can be localized by deformations 
that are coboundaries in cohomology. At technical level the localization of the holomorphic loop 
equation for twistor Wilson loops requires that the backtracking arcs end with cusps supported on 
the singular divisor of the surface operators. 

Our new holomorphic loop equation for twistor Wilson loops is derived using the standard 
technique of the celebrated loop equation of Makeenko and Migdal |Q ^ invented long ago, but 
the crucial difference is that the integration variable that gives origin to the loop equation in our 
case is not the gauge connection but instead the ASD field of our non-SUSY version of the Nicolai 
map, in a holomorphic gauge defined by a further change of variables. The resulting loop equation 



resembles for the cognoscenti the holomorphic loop equation of Dijkgraaf and Vafa [ ]55| , pq , p7[ ] 



for the holomorphic chiral ring of ^ = 1 SUSY gauge theories []58[]. 

The homological localization of the holomorphic loop equation completes the analogy with 
Nekrasov computation. As well as the prepotential, i.e. the effective action in the low energy 
sector as a function of the condensates of the ^ = 2 SUSY theory, is computed by cohomological 
localization of 1 , so the \aige-N effective action of the YM theory in the twistor sector, as a function 



^^The adjoint action of the global gauge group leaves invariant the center, Z/y. 
The classical action scales as k times the square density of surface operators, p , and the renormalized square 
density scales as j. 

21 



Glueballs in large-N YM by localization on critical points Marco Bochicchio 



of the condensates of surface operators of Za^ holonomy, is computed by homological localization 
of 1 , the trivial twistor Wilson loops. 

The twistor sector is defined by correlation functions obtained by holomorphic/antiholomorphic 



fusion a la Cecotti-Vafa \ p% |60| ] of the holomoiphic/antiholomoiphic ASD curvature of surface op- 
erators restricted to the aforementioned Lagrangian submanifold of twistor space of complexified 
space-time with Minkowski signature. 

On the ultraviolet side, a striking result that follows from the localization on fixed points is 
that the large-A^ beta function for the Wilsonian coupling constant in the ASD variables is one- 
loop exact, because for twistor Wilson loops, precisely because of the localization, a certain kind 
of saddle-point approximation turns out to be exact. Thus the quantum corrections for these ob- 
servables are completely accounted by functional determinants whose diagrammatic expansions 
contains only one-loop Feynman graphs and possibly the logarithm of the powers of the Pauli- 
Villars regulator of zero modes. Thus in the large-A^ pure YM theory Eq.(l.l) can be replaced by 
the much simpler: 

Aiv=Aexp(-— 1^) (1.50) 

where Ai^ is the RG invariant scale in the Wilsonian scheme in the ASD variables. 

At the same time we show that, in the regularization scheme of the homological localization 
of the holomorphic loop equation, the one-loop exactness for the Wilsonian beta function implies a 
large-A'^ exact beta function for the canonical coupling of NSVZ type, that reproduces the one- and 
two-loop perturbative universal coefficients. In this scheme the \aige.-N canonical beta function of 
the pure YM theory in the ASD variables is given by Eq.(1.14). 

We should stress that the computation of the canonical beta function depends crucially on 
exploiting the gluing rules for functional integrals in the specific case of the localization on local 
systems. 

In particular we show that these gluing rules coincide with the ones of topological strings via 
the gluing of an associated arc complex ^l, 52, 63,^, 65]. Indeed the homological localization is 



based on "the most local part" of the homology of the essential arc complex of a punctured sphere. 

The requirement that the homology be essential, i.e. the exclusion of arcs that can be deformed 
to a puncture, rules out the local relative homology of compact support around a puncture. In fact 
"the most local part" of the essential homology that is relevant for us is the essential homology of 
the arc complex with no polygons, that involves only links ending with two different cusps, one 
in the divisor at the ultraviolet and one in the divisor at the infrared. Naively the Wilsonian beta 
function in our scheme is a purely ultraviolet concept, and therefore does not distinguish between 
the local homology and the essential homology of links. Yet, the canonical beta function, that in 
our scheme involves infrared physics too, does. 

In addition, while fixed point arguments suffice to display localization directly in large-A^^ 
YM theory by the change to the ASD variables in the functional integral and the dense hyper- 
Kahler resolution on local systems, the corresponding effective action, i.e. the logarithm of the 
density of the localized measure, has intrinsic finite holomorphic ambiguities due to the freedom 
of making holomorphic changes of variables and possibly holomorphic anomalies at loci where the 
holomorphic change of variables may develop singularities. We are able to fix these holomorphic 
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ambiguities only using the localization by homology of the loop equation, via the choice of the 
holomorphic gauge that is necessary to write down the holomorphic loop equation. 

On the infrared side, the holomorphic gauge in the loop equation is essential, because the mass 
gap and the glueball spectrum occur precisely because of a non-trivial Jacobian from the unitary to 
the holomorphic gauge in the effective action. 

We end this introduction with some loose heuristic considerations as to why the line of thought 
of this paper may be able to overcome the main difficulties of the ultraviolet and infrared problem 
of YM in the restricted sense specified above. 

While the mass gap problem as formulated in full generality for every correlation function 
and for every compact gauge group in [jT]] appears presently almost hopeless in our opinion, the 
program of solving the SU (N) YM theory in the large-A^ limit has attracted considerable attention 
and efforts. 

A promising avenue is to find an equivalent string theory [ [TT| , |6^ by effectively resumming 
't Hooft perturbative double expansion in powers of g and N^^ |^]. In this string theory the v.e.v. 
of any Wilson loop of the YM theory in the large-A'^ limit would be computed by a string diagram 
that is a disk with the loop as boundary. No other interesting observables, but the /?G-invariant 
condensates, exist at the leading large-A'^ order because of the factorization of the v.e.v. of local 
normalized gauge invariant operators. 

Now, on the field theory side, the knowledge of all the Wilson loops that would be implied by 
the string solution contains a vast information in a mathematical sense. Indeed it has been known 



for some time that the ambient algebra of the master field [ ]67[ ] that solves the large-A'^ Makeenko- 
Migdal loop equation for ordinary Wilson loops [^, ^ is the Cuntz algebra with four generators 
[ l68| , |69| , [70| , [tTI , 72, ^, 74] whose Fock space representation is known to be of type //i but not 



hyperfinite [l75[], i.e. not the weak limit of matrix algebras. Indeed such Fock representation is 
isomorphic to a free group factor with the same number of generators, which is the main explicit 



example of the "elusive" type IIi non-hyperfinite factors [75]. It is clear |23] that obtaining the 
relevant non-hyperfinite information would be extremely difficult in case the von Neumann algebra 
generated by the actual solution shares with the ambient algebra the non-hyperfinite character. 

On the contrary, to the next to leading ^ order, the connected two-points correlation functions 
of local gauge invariant operators are conjectured to be the most simple as possible: a sum of 
an infinite number of propagators of free fields ^Jdj\, saturating the logarithms of short distance 
perturbation theory [[77|]: 

f<^LTrF^p(^)L^TrF^p{(^)>,„^^ ,'>Vx = £-^ -/(;.)/ log(^) (1-51) 

Now any string solution, as it is usually meant, cannot avoid to solve the leading order problem 
for the Wilson loops, in order to solve the much simpler looking subleading problem for the free 
glueball spectrum. This makes such a general, \aige-N exact, string solution very difficult in our 
opinion. We may wonder as to whether we can solve the easy looking subleading problem for the 
free glueball spectrum avoiding and thus loosing the information about the hard looking problem 
for the Wilson loops. Our answer is positive to a certain extent: we construct trivial Wilson loops, 
the twistor Wilson loops indeed, whose v.e.v. is 1 in the leading large-A^ limit. However, they admit 
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non-trivial ^ corrections and thus morally they couple to a certain non-trivial sector of the Xwge-N 
theory. 

By the way, in this restricted sense we believe that there is also an explicitly solvable string 
theory, that captures the sector of YM accessible to the twistor Wilson loops defined in this paper. 
The outlook for this twistor string is described in the conclusions. 

Coming back to the field theoretical framework, we use our localization theory to localize 
the twistor Wilson loops, that are in the "homology of 1", precisely in the dual sense to which 
Nekrasov localized the "cohomology of 1" to get the Seiberg-Witten prepotential. Indeed, although 
the prepotential is obtained by a "trivial" cohomology, it allows one to reconstruct the low energy 
effective theory. Precisely in the same sense, since twistor Wilson loops live in a "trivial" homology, 
they can be localized by suitable deformations. Yet, the interesting information is contained in the 
effective action, i.e. in the localized measure. The quadratic small fluctuations of the effective 
action around the localized loci of the measure furnish the glueball spectrum in the twistor sector. 

2. Synopsis 

We summarize here the main technical arguments in a logic order and the main results. 

In the prologue we describe in some detail the computation of the beta function in the following 
cases. The one-loop beta function in YM by the usual background field method. The NSVZ beta 
function in ^ = 1 SUSY YM by cohomological localization in the Nicolai variables. The one-loop 
YM beta function for the Wilsonian coupling in the ASD variables by the usual background field 
method. 

These concrete examples are used to furnish a comparison with the computation of the YM 
beta function by our new localization. This section contains many definitions and computational 
technicalities that are referred to throughout the whole paper. 

We define also the change of variables from the gauge connection to the ASD curvature in the 
pure YM case. In particular we show that, since in pure YM the Jacobian to the ASD variables is 
not cancelled, as opposed to the ^ = 1 SUSY YM case in the light-cone gauge, a multiplicative 
Z renormalization of the ASD field occurs. In sect.(12) this Z factor is related to the anomalous 
dimensions of a large class of composite operators that occur as scalar polynomials in the ASD 
curvature in the one-loop integrable sector of \aige-N YM. 

In sect.(4) we define twistor Wilson loops in non-commutative gauge theories. The twistor 
Wilson loops are defined on a non-commutative deformation of space-time, that is used as a tool 
to define the large-A'^ limit much in the way Nekrasov used a non-commutative deformation as a 
tool to regularize the instantons moduli space. We recall some features of non-commutative gauge 
theories that we employ in the following sections. 

We display the following properties of twistor Wilson loops. The v.e.v. of twistor Wilson 
loops is fiber independent and trivially 1 in the large-A'^ limit. In addition twistor Wilson loops are 
supported on Lagrangian submanifolds of twistor space. 

In sect. (5) we show that the curvature of the twistor connection is of purely ASD type and we 
describe the localization in large-A'^ pure YM theory of twistor Wilson loops on the fixed sheaves 
of a semigroup acting on the fiber of the Lagrangian fibration which twistor loops are supported on 
and contracting the support of the functional measure in the ASD variables. 
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We refer to this kind of localization as the quasi-localization lemma, since the resulting lo- 
calized measure is still represented by a residual functional integration on a certain complex path, 
supported on distribution valued sheaves in fact, rather than by a sum over fixed points. 

The quasi-localization lemma is a purely formal computation that involves a quite disputable 
formal exchange of the order of limit and integration. However, the exchange of order of limit and 
integration is justified in sect.(8) in the \aige-N limit and for a lattice version of the Nicolai map 
sect. (7), that allows interpreting the Nicolai map as the hyper-Kahler reduction on a dense set in 
function space in the sense of distributions. 

We write down the corresponding effective action, i.e. the logarithm of the density of the 
localized YM measure. Because of the residual complex integration the effective action has an 
ambiguity by holomorphic change of variables that we can solve only through the loop equation of 
sect.(6). The fixed sheaves in the quasi-localization lemma are characterized by the vanishing of 
two of the three ASD fields of the non-SUSY Nicolai map. 

In sect.(6) we write the holomorphic loop equation for twistor Wilson loops. The holomorphic 
ambiguity of the effective action of sect.(5) is fixed by a change of variables to the holomorphic 
gauge, necessary to write down the holomorphic loop equation. It is precisely the Jacobian to the 
holomorphic gauge that generates the glueball potential. 

This implies that the glueball potential in the holomorphic/antiholomorphic sector defined by 
twistor loops in the fundamental and conjugate representation must be singular at the fixed points, 
as it is indeed, for the theory to have a mass gap, since the contribution of the Jacobian to the 
effective action is formally the logarithm of the square of a holomorphic function. 

In sect. (7) we introduce a regularization of the large-A^ functional integral by integrating on 
"infinite-dimensional local systems" ^' on non-commutative space-time. 

The idea of integrating on local systems associated to an arbitrary parabolic divisor appeared 



for the first time long ago in a paper by us [g2|, g3|], by embedding the Hitchin fibration [g4|, g7|] in 



the YM functional integral, and physically corresponds to integrating over surface operators ante 



Utteram [28|. We employ Morita equivalence [78| to reduce to the case of ordinary space-time for 



finite rank bundles. Thereafter we reconstruct the large-A'^ limit of YM as an inductive limit on the 
finite rank local systems. 



Following the mathematical literature [ ]52| , |24] , [79| , |80| , |82| , |83[ , [27| , |84| , |85|, |87|] we discuss 



the topological, holomorphic and differential geometric features of the finite rank local systems 
[ ^ , ^ ]. As topological objects local systems are representations of the fundamental group of 
a punctured Riemann surface. As holomorphic objects they are holomorphic connections with 
regular singularities. As differential geometric objects they are parabolic harmonic bundles, i.e. 



parabolic Hitchin bundles [52] equipped with a harmonic metric by a Hitchin-Kobayashi corre- 
spondence [^4]]. Remarkably in our setting the harmonic bundles arise as the hyper-Kahler reduc- 
tion [B^] induced by our version of the non-SUSY Nicolai map. 



Physically the hyper-Kahler reduction [52] is a resolution dense in function space [22], in a 



neighborhood of the fixed sheaves, of the ASD field [ 23 ] as a linear combination of two-dimensional 



delta distributions supported on a lattice of surface operators [^]. These are local systems that oc- 



We write it in quotes because these infinite dimensional objects admit unstable finite dimensional subbundles, thus 
violating a fundamental property of finite dimensional local systems. 
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cur in the mathematics and physics Uterature for completely different reasons, among which we 



mention the Hitchin-Kobayashi correspondence []84|], non-abelian Hodge theory \ \J9l |80| ], twistor 



D-modules [|87|], and last but not least the physics version [ |28| ] of the geometric Langlands corre- 
spondence [|8H8|, ^, mH, ^, |4[ 0]- 



In sect.(8) we get our localization on fixed points. Indeed we combine the quasi-localization 
lemma of sect.(5) with the idea of sect.(7) of integrating on local systems to get localization on 
fixed points. Reducing to finite dimensional local systems needs Morita duality and is allowed 
implicitly by the triviality of twistor Wilson loops. 

In fact the quasi-localization lemma depends on the aforementioned disputable formal ex- 
change of the order of limit and integration. This is justified by showing that, on the set described 
by the hyper-Kahler reduction, a gauge exists in which the lattice theory is locally abelian, all the 
remaining non-abelian degrees of freedom being zero modes associated to the moduli of the local 
system. Therefore the fluctuation of the eigenvalues of the theory are suppressed in the laige-N 
limit. Of course the integral of the limit is the localized measure in the ASD variables, while the 
Hmit of the integral is the original YM measure on the gauge connections. 

In particular the choice of the approximating sequence by finite dimensional local systems 
(i.e. stable bundles) is essential for localization on fixed points. We show that the fixed manifold 
restricted to the dense hyper-Kahler locus of local systems is a Lagrangian submanifold of the 
moduli space of surface operators. In addition we show that, assuming that the gauge group is 
unbroken, the fixed manifold is in fact a collection of fixed points represented by surface operators 
with Zjv holonomy. This localization is our analog of the Duistermaat-Heckman locahzation. 

In sect.(9) we use it to compute the Wilsonian beta function of the large-A^ YM theory. 

In sect.(lO) we get homological localization of the holomorphic loop equation for the twistor 
Wilson loops by means of the lattice version of the holomorphic loop equation obtained integrating 
over the local systems. 

The triviality of twistor Wilson loops plays a key role here, since it allows arbitrary defor- 
mations of the twistor Wilson loops without changing their expectation value. In this section the 
localization on fixed points in not a consequence of the assumption that the gauge group is unbro- 
ken, but a consequence of the reduction of the loop equation to a critical equation. 

The localization is obtained deforming the twistor Wilson loops by backtracking arcs ending 
with the cusps of the local system. This is the localization by homology that combines the holo- 
morphic loop equation of sect.(6) with the idea of integrating over local systems of sect.(7). This 
is our analog of cohomological localization by deforming by a coboundary. 

In particular we derive the glueball potential by computing the Jacobian to the holomorphic 
gauge of sect.(6) for the local systems of sect.(7). 

At mathematical level homological localization involves the essential arc complex of a punc- 
tured sphere and a combinatorial model of the gluing of the arcs developed in the mathematical 



literature []61|, |62|, p3| ], in turn inspired by (topological) strings [|64j, p5| ]. 

The version of the localization via the loop equation is the one that has, in our opinion, more 
chances to hold in a strictly mathematical sense. The loop equation for twistor Wilson loops occurs 
as a formal Schwinger-Dyson or Ward identity derived imposing that the integral of a functional 
derivative vanishes in function space in the ASD variables in the holomorphic gauge. 
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The left hand side of the loop equation contains the effective action of the theory that implies 
a one-loop exact Wilsonian beta function. The right hand side is still divergent, but it can be 
regularized in a gauge invariant way by analytic continuation to Minkowski space-time and, by 
deforming the loop, it can be made to vanish. Thus, despite the loop equation is obtained only 
as a formal identity (as the Makeenko-Migdal equation is), its solution is defined via its would-be 
properties, essentially the fact that it allows analytic continuation to Minkowski space-time. 

The glueball spectrum occurs only after analytic continuation to Minkowski space-time of the 
effective action renormalized in Euclidean space for fluctuations of surface operators restricted to 
the Lagrangian submanifold. 

The basic idea is that once the laige-N localization is obtained for twistor Wilson loops, that are 
non-local extended objects, the localized effective action is used to compute physical fluctuations 
of local operators restricted to certain channels. 

In sect.(l 1) we compute the canonical beta function of l&rge-N YM by means of our holomor- 
phic loop equation restricted to the local systems of sect. (7) and we check agreement with the first 
two universal perturbative coefficients. The result depends crucially on the gluing rules for local 
systems. 

In sect.(12) we display our main result about the mass gap and the glueball spectrum using the 
effective action of sect. (6) together with its extension to the hyper-Kahler locus in a neighborhood 
of the fixed points of sect.(7). 

To do computations we employ the local model of the singular part of the connections with 
regular singularities around surface operators of sect.(7). We show that at the (renormalized) critical 
points the local model is in fact asymptotic for large N. We find a trajectory with mass squared 
exactly linear in k and residues at the poles determined by the multiplicity of the eigenvalues of 
the ASD curvature at the fixed points and by certain finite counterterms as a function of /c in a 
neighborhood of the fixed points. 

We display in some detail how the glueball spectrum for the trajectory in the twistor sector 
follows from the effective action. We also use our explicit solution to check the long-standing 
conjecture that the sum of pure poles in the large-A'^ limit on our trajectory saturates the logarithms 
that occur in the glueball propagators in perturbation theory. 

In the Wilsonian scheme we reproduce a factor of a logarithm that occurs in perturbation 
theory, that arises in our scheme by the /?G-invariant spectral sum over the glueball. 

In the perturbative canonical scheme the sum of a logarithm (with a coefficient that can be 
normalized to 1) and of its square occurs. The coefficient of the square of the logarithm depends 
on the operator and is related to the anomalous dimension. 

We observe that in the canonical scheme the multiplicative Z renormalization, that occurs be- 
cause of the Jacobian to the ASD variables mentioned in sect.(3.4), implies through the localization 
of composite surface operators the same anomalous dimensions as for the operators associated to 
the ground state of the Hamiltonian spin chain in the thermodynamic limit, which is known to fur- 
nish the one-loop anomalous dimensions of long local gauge invariant operators in the one-loop 
integrable sector of \aige.-N YM. This sector is made by SD or ASD fields. 

This cannot be the whole story, since we get from surface operators just one trajectory. In 



sect.(13) using by now standard results in mathematics [85, 99], we extend the hyper-Kahler re 



duction induced by the non-SUSY Nicolai map to twistor connections with wild singularities, i.e. 
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poles of any order. We suggest that such an extension corresponds physically to the more realistic 
case of an infinite family of Regge trajectories of increasing spins. We write the basic definition of 
the functional integral on wild surface operators but explicit computations are left for the future. 

In sect.(14) we summarize our conclusions and we outline some features of the twistor string 
conjectured to be dual to the YM theory restricted to the sector defined by the twistor Wilson loops 
of this paper. 

3. Prologue 

3.1 One-loop beta function of YM by the background field method 

This computation is now completely standard, but since it is not easily found in textbooks in 
the form that we will need in the rest of the paper we display it here in some detail ^^. The basic 



philosophy is as in [[Tq]. The partition function of pure SU{N) YM is: 

Z=[5Ae-^ (3.1) 

where: 



^^ Jd'x trf{F^^f = ^Jd'x (F^p)' (3.2) 



2g 



The sum over repeated indices is understood. The action has been rescaled by a factor of N in 
such a way that the theory admits a non-trivial large-A'^ limit. The coupling constant, g, is the 't 
Hooft coupling related to the YM coupling, gYM, by g^ = gyM^- The normalization of the action is 
appropriate for a gauge connection in the fundamental representation of the Lie algebra of SU (N) : 

A„=A«r' (3.3) 

with the Hermitian generators in the fundamental representation normalized as: 



The curvature of the YM connection is 





trf{T^T^) = - 5"'' 




;tion 


is: 








[Da,Dp]=iFap 




Fap 


= daAp 


-dpAa + i[Aa,Ap] 




= daA"p 


-dpA"„-f'^A',A 



(3.4) 



(3.5) 
(3.6) 
(3.7) 



where Da = da + iAa is the covariant derivative. To perform the one-loop computation of the 
effective action it is convenient to split the gauge connection into a classical background field and a 
fluctuating quantum field Aa=Aa + 5Aa- The Fourier transform of the quantum field is supposed 
to be supported on momenta much larger than the momenta of the classical background field. The 



We would like to thank Luca Lopez for working out a detailed version of this computation during our course at 

SNS. 
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gauge-fixing is performed by the Faddeev-Popov procedure. It is convenient to choose the Feynman 
gauge with respect to the background gauge field A^: 

D„5Aa-C = (3.8) 

where we denote by the dot the adjoint action in the Lie algebra: 

Da{A)-8Ap = da8Afi+i[Aa,5Ap] (3.9) 

{Da{A)-5Ai,Y = daSAj, -f'^'^A'^SA'p = D''^{A)5A'p (3.10) 



where 



D''^{A) = da5'''-f^'A^ 



a 
dad'^' + iA't 



'a 



with A"^ = if^'-'A^ and [T'',T''] = if^'^T'^. C is an auxiliary gaussian field whose covariance is 
chosen in such a way to cancel a longitudinal term in the YM action quadratic in the fluctuating 
field, by adding ^ J d'^xtrf{D-^8AaY to the action. Quantities such as Da are evaluated at the 
background field Aa- The gauge fixed partition function reads: 

Z= j 8A8Cex^{-SYM)Det{-h-^)exg{-— f d'^xtrf{C^))5{D-„5Aa-C) (3.11) 

where we have inserted the Faddeev-Popov determinant of (minus) the Laplacian in the background 
field: 

-A^ = -A2 - idaA„ - 2iA„da +A„A„ (3.12) 

As a consequence the gauge-fixed action is: 



/ 



^ trfF^piA + 5A) + ^ trfiD„5Aafd'x (3.13) 



and the one-loop partition function reads: 

Zi-io„p = e-r'-'-"(^) = e-^™(^)Der i/2(-A^5„/3 - 2iF^p)Det{-A-^) (3.14) 

where ri_;oop(A) is the effective action for the background connection, Aa, to one-loop order and 
F^pi...) = [Fap,..]. 

It is very instructive to understand the origin of the spin term, —liF^n, in the first functional 
determinant of Eq.(3.14). By the splitting of the connection into Aa = Aa + 5Aa the curvature 
decomposes as follows: 

FapiA + 5A) = FapiA)+DaiAy5Ap-DpiAy5Aa + i[5Aa,5Ap] (3.15) 

Performing the square we keep only up to the quadratic terms in 5Aa in the action, since we are 
doing a one-loop computation, we understand integration on space-time in the following and we 
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freely integrate by parts. We use the equation of motion, D«(A)F„jg(A) = 0, to eliminate the Unear 
term in the action. Therefore we get: 

Flf^{A + 8A)=Flp{A) + {Da{Ay8A^-D^{AydAaf + liF^fi{A)[5AaMii\ 

= Flfi{A) +2{Da{Ay SAnf -IDaiA)- 5AiiDii{Ay 5Aa + 2iFaii{A)[5Aa,5Ap\ (3.16) 

Using 

DaDp=DpDa + iF„p (3.17) 

trf{5Ap[F„p,5Aa]) = trf{F„p[5Aa,5Ap]) = -trf{5Aa[F„p,5Ap]) (3.18) 

the quadratic form in Eq.(3.16) becomes: 

trf{{D-„5Ap-Dp5A„f + 2iFap[5A„,5Ap]) 

= trf{-2dAaA)^dapdAp+25ApD-o.D-p5Aa+2iFap[5Aa,dAp]) 

= trf{-2dAaAASapSAii+25ApDpD-^5Aa + 2idAiiF^pdAa + 2iF„p[5Aa,dAii]) 

= trf{-25AaAA5ap5Ap-2iD„5Aaf-4i5Aa[Fap,5Aa]) (3.19) 

where we skip the label of the background field since no confusion can arise. In the Feynman gauge 
the second term in the last line is cancelled by the gauge-fixing. Finally the quadratic form written 
in components becomes: 

trf{5Aa{-2A-^5aii - 4/ adFap)5Ap) 

= trf{5A''j"{-2AA5ai}-4iadFo,p)dApft'') 

= ^Kii-^ASap-^i adFf^p)5Apy 

= 5A"„{-Aa5„p - 2/ adFo.prSA'^p 

= 5A«„ ( - {AAr5„p + 2r'''F^p)5A'^p (3.20) 

where adF^^pi...) = [i^„/3,-], {adF^^pr = if^'F^p and 

(A^r = iDa)'^{Da)'''=A'5'"'- + idaA'},' + 2iAl'da-Al'Ai' (3.21) 

QED. 

The following identity holds: 

Det-'/^i-AA5ap -2i adFap) = Der'/^i-AA6ap)Der'^\\ -2i{-AA)-' adF„p) (3.22) 
The first factor gives: 

Der'^^{-AA5ap)=Der^{-AA) (3.23) 

Therefore the one-loop effective action reads: 

g-r,_,„„,(A) ^ g-s™(A)^g^-i/2(^j _ 2/(_Aa)"^ adFap)Der\-AA) (3.24) 
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The first determinant is the spin contribution while the second determinant is the orbital contribu- 
tion. 
We can factorize away a trivial infinite constant from the orbital contribution: 

Der^ (-Aa) = Der\-^^ - idaAa - liAada +AaAa) 

= Det-\-A^)Det-^ (1 + {-A^y^ {-idaAa - HAada +AaAa)) (3.25) 

where the operators occurring in Eq.(3.25) are now defined by Eq.(3.21). Using 

Det{\ +M) = e^'-l°g(l+«) = JrM-Tr(Mfl2+... ^^ ^6) 

at the lowest non-trivial order we get: 

Det-\-AA)=Det-\-A)txp{-Tr{{-A)-\-idaAa-2iAada+AaAa))) 

exp (rr((-A)-i {-idaAa - liA^da + A„A„)(-A)-' (-/5„A„ - HAada + A„A„))/2) 

(3.27) 

where the trace is over the space-time, the Lie algebra and the vector indices. The term Det^^ {—^) 
is an irrelevant constant while the Lie algebra trace of the term linear in A„ vanishes. The term 
Tr{{—A)^^AaAa) is a quadratically divergent tadpole that cancels in any gauge invariant regular- 
ization scheme, since it would give rise to a mass counterterm for the gauge connection. Therefore 
it can be ignored. There remains an interesting divergence: 

Det-\-AA)^txp{Tr{{-A)-\idaAa+2iAada){-A)-\idaAa+2iAada))/2) (3.28) 
that evaluated in momentum space leads to: 

A r d'k f d^p {2pa-ka){2pp-kp) 

where the trace tr on the Lie algebra indices refers to the matrices defined in Eq.(3.21). The 
logarithmically divergent part of the integral over d'^p has to be transverse in such a way that 

j- d^p {2pa-ka){2pp-kp) . 2^.,2c , , n, ,o onx 

J (W I^IF^' ^ ^ ^^ 5ap-kakp) + ... (3.30) 

where the dots stand for the quadratically divergent part that can be ignored because of the afore- 
mentioned reasons. Taking the trace over the vector indices one gets: 

d4p 4p^ + k^_4pk _ , , 
(27r)4 p^{p-kr -^'^^HM + - (3.31) 

We are interested in extracting the logarithmic divergencies by expanding the denominator in pow- 
ers of k/p up to the appropriate order: 

d^p 4p^ + k^-4pk f d^p 4p^ + k^-4pk 



(271)4 p^{p-kY J {2nYp\\ + {k^-2kp)/p^) 

d'^p -4k^ + k^ + S{pk)^/p^ f d^p k^ 2k^ A ^ ^ 



f d'^p k^ 2k^ A . , . 

.} 271 4 p4 4^2 &„ V 



(27r)4 p'^ J {2nY p^ {4n) 

(3.32) 
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where we have replaced 

1 2 
PaP^ -^ -^P Saji (3.33) 

into the integral and similarly {pk)'^ — )■ p^k^/4 and we have regularized 

-^=1k' -^^In" log - (3.34) 

J p^ Jo p At 

Therefore: 

Hence the orbital contribution to the beta function is: 

Det-\-AA) ~ exp {-j^^log- J j^triAcc{k^5ap-kcckp)Ap)) 

1 A /■ d'^k 

= ^'^P(-3(W^°'Mi(W^«^'''"'^-'"''^^^'^'^ 

= '^P("3(^^°^^^^'^/'^'^^^"''^'^ ^^'^^^ 

where in the last step we used: 

and 

1 /" 1 /" d^k 

-J d'xF^af^F^^p r^-J ^—^AU-k)A^k){k'5,p-kakp) (3.38) 

at quadratic order. 

Now we have to compute the spin contribution to the effective action. Since trF„p = up to 
quadratic order in F^p we get: 

Der 1/2(1 -2/(-Aa)-i adFap)=exp{-Tr{i-AA)-'adF„pi-AA)-'adFpa)) (3.39) 
At lowest order (— A^) '^ (—A), therefore 

Tr{{-Ay^adFap{-Ar'^adFpa) = / d'^x j d'^ytr{G{x-y)adFap{y)G{y -x)adFpa{x)) 

= -N j d\ j d'yG{x-yfF"^p{y)F-^p{x) (3.40) 

where in coordinate space 

and 

tr{adF„padFpa) = (adF^^pr (adFp^r = if''Pafiif'"Ppa = -^KpFa^ (3-42) 
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Assuming that the background field carries momenta much smaller than the fluctuating field we 
can expand F^p {y) = F„^ (x) + . . . by Taylor series and keep the first term since we are interested 
only in the divergent terms. Thus defining z = x — y 

AN A 



'^^°^^^1/''^(^«/^)' ^'-''^ 



Therefore at this order the divergent part reads: 



Det-'/\l-2i{-A^)-'adF„p)^cxp{j—^log{-f-ld'x{F^pf) (3.44) 

Finally the local part of the one-loop effective action reads: 

Therefore the bare coupling constant, g(A), renormalizes as: 

1 1 11 1 A 

+ ^zr^log- (3.46) 



2g2(A) 2g2(At) 3 (47r)2 ^^ 
or 



^(A) = ^At^ (3.47) 



that is the solution at one loop of the equation that defines the j8 function: 



''« = 51^ = -ft«' + ' 



^ = W^ O"*) 



Eq. (3.46) can be also written as: 



A e 2/Sos"(A) =i_ie ^fts'C' (3.49) 

Thus the combination: 



AyM = Ae 2fc?(I) (-350) 

is independent on the cutoff A and it is a renormalization group invariant at one loop. 
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3.2 A SUSY interlude: cohomological localization by the Nicolai map in ^ = 1 SUSY YM 

Shortly after Nicolai discovered [Q, |8|] that the vanishing of the vacuum energy in an unbroken 
supersymmetric theory implies the existence of a change of variables whose Jacobian formally 
sets the functional integral in ultralocal form, De Alfaro, Fubini, Furlan and Veneziano ^^ [^, [l0| ] 
worked out explicitly the Nicolai map in the case of ^ = 1 SUSY YM. 

They found that in this case the Nicolai map is actually the change of variables from the 
gauge connection to the ASD part of its curvature in the light-cone gauge, with the property that its 
Jacobian cancels the gluino determinant. 

In this section we reconsider the Nicolai map of c/K = 1 SUSY YM paying particular attention 
to the fact that, while generically in function space the aforementioned cancellation occurs exactly, 
in a renormalizable but not finite supersymmetric quantum field theory such as ^ = 1 SUSY YM 
there should exist loci in function space where the cancellation occurs in fact only up to zero modes. 

Indeed if it were not so the theory would be in fact mapped into a theory of free fields with 
zero beta function. 

It is quite clear that the Jacobian of the Nicolai map develops zero modes precisely at loci 
in function space where the Nicolai map fails to be one-to-one. If these loci are characterized by 
moduli then there is a continuous family of zero modes and the Pauli-Villars regularization of these 
zero modes in the functional integral furnishes in general some contribution to the beta function of 
the theory, thus resolving the puzzle that the Nicolai map maps formally the theory into a theory of 
Gaussian fields with vanishing beta function. 

Understanding the distribution of these zero modes as a function of the Gaussian random 
field which the theory is generically mapped on is in fact a non-perturbative problem seemingly as 
difficult as performing the functional integral in the original variables. 

However, we point out in this section that, thanks to the tautological nilpotent Parisi-Sourlas 
[ p8| , p^ ] BRS symmetry ^° associated to the cancellation of the Jacobian with the gluino deter- 
minant, the partition function with the insertion of certain BRS invariant operators necessary to 
saturate the zero modes of the gluino determinant is in fact localized by cohomological localization 
(sect.(l)) on those (Euclidean) instantons ^^ that can be analytically continued to ultrahyperbolic 
signature (this constraint arises because the cancellations due to the Nicolai map actually occur 
only in the light-cone gauge). 

Thus because of the localization the occurrence of the zero modes for these special BRS- 
invariant observables can be understood semiclassicaly as they coincide with the moduli of the 
instantons. 

The immediate consequence of this localization is an exact formula for the beta function of ^ 



= 1 SUSY YM, that quite obviously turns out to be the NSVZ beta function [ |15| , |51| ], by an almost 
verbatim reproduction of their original computation. 



^^We would like to thank Gabriele Veneziano for several discussions about the Nicolai map over the years and at the 
GGI. 

^"We may consider this BRS symmetry as a remnant of the supersymmetry after gauge-fixing. 

^^We refer here to instantons as configurations satisfying the SD equations, without any implication about being 
defined on 5 . In fact in the present framework the instantons are naturally defined on 5 x S , because of the analytic 
continuation from Euclidean (4,0) to ultra-hyperbolic (2,2) signature, that can be handled by twistor techniques | |98| | 
(see also the Appendix in ip9l|). 
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We can now start to work out the details. It turns out that to get localization we need only 
the information that the gluino determinant cancels the Jacobian of the Nicolai map. Firstly let us 
suppose that there are no zero modes. Using the identity: 

Tr{F^I^) = Tr{F^^f/2 + Tr{FafiFaf,) (3.51) 

the partition function reads: 



[j exp{- 






2^2 



^2 j Tr{F,f^F-''^)d\)Det{^)8A]^^^^ (3.52) 



where fields live in the adjoint representation of the Lie algebra of SU{N). In Eq.(3.52) we have 

just expressed the existence of the Nicolai map by inserting its inverse Jacobian, Det{-^), in 
place of the gluino determinant. Thus the theory is mapped into a theory of free fields: 

Z = /exp(-^^^-^/r.(F-,F-«'^)A)5F-, (3.53) 

that holds generically in function space where zero modes do not occur. In fact, taking into account 
the zero modes, we get: 









e^P( 2^ - s^2 J Tr{F,pF-''')d'^)5F^p (3.54) 



[/exp(-^^^-^/j4xrr(F„^F-«/^+p„/3^r7,))5A5p5r7],^^„ (3.55) 



2^2 

where the extra factor is the contribution of the Pfaffians of the bosonic and fermionic zero modes 
^% and of the associated Pauli-Villars regulator, whose origin is explained below. Going back to 
Eq.(3.52) we can write it in a more suggestive form introducing anticommuting fields, (Paj3, TJ/): 

2g2 
The non-topological term can be rewritten as: 

/I r U F 

i-^Jd'xTriE^pE-l^ + iE^pF-P-ip,p^7i,))dE-dAdpd7i]^^^, (3-56) 

The functional integral here has to be interpreted as either in Minkowski space-time (3,1) or in 
ultrahyperbolic signature (2,2), since otherwise the light-cone gauge does not exist. In Minkowski 
signature an overall factor of / in front of the action is understood but not explicitly displayed. In 
ultrahyperbolic signature the Gaussian integral is defined by analytic continuation. In this form the 
partition function enjoys the following tautological Parisi-Sourlas BRS symmetry ^^: 

Qbrs^y = ^Y 



^^To the best of our knowledge the existence of this symmetry in jV = 1 SUSY YM has never been related to 
cohomological localization, presumably because it leads to results inconsistent with the NSVZ beta function if zero 
modes are not taken into account properly. However, the zero dimensional version of the Nicolai map, for which of 
course no zero modes occur, has been related to cohomological localization in [pO|. 
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QbrsTIy = 
QbrsPuP = E^p 

QbrsE^P = (3.57) 

with 2|/j5 = 0. The consequence of the existence of this symmetry is that the term E^^E^'^P is 
a coboundary, since E^oE^"^ = Q{papE^"'P). Thus it can be cancelled without changing the 
cohomology class of the integrand. The resulting functional integral reduces to: 



/I r oF o 

{- — j d'xTr{iE-^pF-^P -ip^f,^r],))5E-5A5p5r]]^^^^ (3.58) 



5F^ 
h 
Thus the complete partition function reads: 

Z = [/exp(-^^^)5(F-^)D.K^)5A],^^„ (3.59) 

that expresses the fact that the partition function is localized on those instantons that can be an- 
alytically continued to Minkoswki ^^ or ultrahyperbolic signature, thus violating the assumption 
that there are no zero modes. In case we assume the existence of zero modes from the start we 
have to insert some fermionic contribution to take into account the fermionic zero modes of the 
gluino determinant. This can be done by identifying lexicographically the gluino zero modes, say 
A, where / a spinor index, with T], where / is a vector index: 



[^A ( {TreijTf^iT^j) 2 exp{ 



Nl6n^Q. 



2g^ ' 

1 r oF 

^M-^ j d'xTr{E-^pE-^^ + iE-^pF-^P -iPap-^^,))8E-5A5p5^]^^^^ (3-60) 

Afterwards everything goes through as before and the localized partition function is: 

Z = exp(-^^15£^)A".-¥ / flliyi^ (3.61, 

where we have explicitly displayed the Pfaffians (i.e. square root of determinants) that occur evalu- 
ating the residual integral on the instantons moduli associated to gauge field and gluino zero modes. 
This is precisely the NSVZ result originally found by evaluating the gluino condensate, but for the 
fact that there is a constraint of analytic continuation for the instantons, say to ultrahyperbolic sig- 
nature, that can be handled by twistor techniques [p8|]. In particular the natural framework for the 
analytic continuation is to start with twistors on Euclidean S^ x S^ and analytically continue to 
ultrahyperbolic signature. Doing so the conformal compactification of ultrahyperbolic space-time 
occurs, that is 5^ x S^/Z2, where Z2 acts by antipodal involution (y{x,y) = {—x, —y) ||98|]. Thus the 
analytic continuation defines in fact the double cover, S^ x S^, of the conformal compactification, 
S^ X S^/Zj, of ultrahyperbolic space-time. As a consequence the possible values of the second 



In this case the gauge group must be complexified since SD equations exist only in Euclidean or ultrahyperbolic 
signature for Hermitian connections [PSl^^. 
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Chern class, Q, (the topological charge) can only be even for ultrahyperbolic instantons |98]. Thus 
the gluino condensate cannot be saturated by single instantons. This is perhaps related to the old 
controversy about the strong versus weak coupling evaluation of the condensate (for reviews see 



16|, |47|]), but it is a matter too far away from our main subject to discuss further. 

In any case the NSVZ beta function is tautologically reproduced. We show the computation 
because it is very useful to understand the pure YM case. 

A striking consequence of the localization on instantons is that the beta function for the Wilso- 
nian coupling constant, gw, is one-loop exact since the only sources of divergences are the zero 
modes via the Pauli-Villars regulator: 

167r^6 l6n^Q , rif A 

TT777T = TT77T-("''-t)1oS77 (3.62) 

Now since ni, = 4NQ and nf = 2NQ the result for the Wilsonian beta function follows: 

' ' ' log^ (3.63) 



2g|,(M) Ig^iA) (47r)2 ^^ 
or differentiating with respect to log A: 

dgw 



(9 log A 
with 



Pog'w (3-64) 



Po = j^, (3-65) 

From the one-loop exactness of the Wilsonian beta function it follows the NSVZ formula for the 
canonical beta function. Indeed the renormalization of the canonical coupling is obtained rescaling 
the fields in canonical form, i.e. in such a way that the quadratic part of the action is normalized in 
order to be g independent: 



2gjy J.^Q Pf<gric,gTlc> 



2gl '" ° J.£q Pf<ri„ric> 

^^ 2g2 ^ /^^ P/<T7„T7,> 



where we have defined: 

NlGK^QigA,) Nl6K^QigA,) 



+ {4NQigA,)-2NQigA,))logg (3.67) 



2g2 2^^ 

or 

Differentiating with respect to log A the NSVZ beta function follows: 



2 
{47C) 

QED 



^logA 1-7^! 



(3.69) 
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3.3 'Non-SUSY Nicolai map in pure YM 

Now we have set the stage for the non-SUSY Nicolai map in pure YM. 

In the pure YM case we define our change of variables in a more general way than in the SUSY 
case [^ |I^, by allowing arbitrary gauge-fixing. Indeed, while in the SUSY theory we need the 
light-cone gauge to get the cancellations of the determinants, in the pure YM theory cancellations 
do not occur at all. Therefore there is no point in choosing a non-covariant gauge. 

The YM partition function is (for definitions see sect.(3.1)): 

Z = /exp ( - -^ - ^/^r,(F„-,)Vx)5A (3.70) 

where we used the identity Tr{F^n) = Tr{F^„Y /l + Tr{FapFap) as in the SUSY case. We change 
variables from the connection to the ASD curvature by introducing in the functional integral the 
appropriate resolution of the identity: 

^ = JS(Fap-^ap)S^ap (3-71) 

Thus 

Z = |exp ( - 1^ - ^ jtrf{li-^^fd'x)5{F^p - ^^p)5^,p5A (3.72) 

Exchanging the order of integration we can now perform the integral on the gauge connection 
because of the delta function. The easiest way to do this is defining the delta function as: 

^i^aH -^afi) = li"i N-\e)e'^^^"'f^''"i^'^"P^' (3.73) 

with A'^(£) an irrelevant normalization factor. Thus the Jacobian of the map to the. ASD variables can 
be evaluated as a Gaussian integral for the quadratic form obtained from the expansion to quadratic 
order of: 



^'■/(^a/3-i"«/3) =^^/(^"(^a/3-Ma/3))'~?ry(p-(Di„5A^]))^ (3.74) 



where: 



and 



1 1 



Fali = P-Fap + P^F^p = -F^p + -F- (3.75) 



Therefore, integrating by parts freely and using the same identities as in sect.(3.1), we get: 

trfP'{Dla^Ap]f = trf{D[jAp] - -SapysD[^5As]f 

= trf{2{D[jAp])^-D^jAji]eaprsD[^SAs]) 

= trf{-45AaAA5Aa-4{D„5Aaf+4iFap[5Aa,dAp]-4D-„dApeapysDr5As) 
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= 4trf{-5AaAASAa - (D-^SAa)^ + iFo,p[5Aa,5Ap] +D-^D-^5Apeapys5As) 

= 4trf{-5AaAASAa-{Da5Aaf + iFo,p[5Aa,5Ap] + -{[Dy,Da]y5Apeyaps5As) 

= 4trf{-5AaAA5Aa - {DaSAaf + iFap [5Aa,5Ap] + i[Fps,5Ap]5As) 
= 4trf{-dAaAA5Aa-{D-^5Aaf + iFap[5Aa,5Ap]+iFap[5Aa,5Ap]) 
= 4trf{5Aai-AA5ap+DaDp-2iadiP+F)ap)5Ap) (3.77) 

The partition function becomes: 

Z = |exp ( - 1^^^ - ^ j trf{^i^^fd'x)Der"\-AA5ap +DaDp - iF+p)5^-{3.n) 

The determinant in Eq.(3.78) does not exist unless the gauge is fixed. This is most conveniently 
done in a background Feynman gauge: 



Z = lim / 5A5C5ix e\p{-SYM)e 

£>ef(-A;i)e-^/rf'^"-/(c')5(D„5Aa -C) (3.79) 






As a consequence the gauge-fixed partition function in the ASD variables is: 

r , \6n^NQ N f . _ ,2^4 X 

= y^^P( p ^ J t^fi^ap)^^) 

Det-'l^{-AA5ap - iF+p)Det{-AA)5^- (3.80) 

that because of the argument displayed below can be rewritten as: 

f , \6n^NQ N r , ,2 ,4 N 
2 =\^M g^-^J^ni^^apfd'^) 

Det-^l^{-AA8ap-ipi'^p)Det{-AA)A">'^^''^(o'^8^- (3.81) 

where g) is a Kahler form on the moduli induced by a Kahler form on the connections. A possible 
choice for (O is: 

1 /■ 4 5A^ 5Ai 5Au 5Au 



(Oi = — d^xtrf {^^-^^)5mi^5mk (3.82) 

2ti J orrii onik ontj onik 

As in the SUSY case zero modes have to occur, but in the pure YM case there are also other 
contributions to the beta function due to the lack of cancellation of determinants (sect.(3.4)). 

Let us introduce the matrices Oa = ( 1 , it) with T the three Hermitian Pauli matrices and the 
self-dual and anti-self dual matrices (we use the same notation as in ||4^): 



(^ap = -^{(^aOp - OpOa) 
(^ap = -^[(yaOp - apGa) 

(3.83) 
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The three variations [H^] : 



can be rewritten as: 



and the forth variation: 



as: 



p-{D[^5Ap]) (3.84) 

iTi)ip-'"'dAaa = Tr{Ti p5 A) (3-85) 

Da8Aa (3.86) 



p""5Aaa = Tr{p-5 A) (3.87) 

in such a way that the four variations can be written together as: 

P^5A^y (3.88) 

where we have defined: 

P = <yccDa 

P = OaDa (3.89) 

Therefore the Euclidean invariant positive semidefinite quadratic form: 

£ ?r/(p-(Di„5A/3]))2 + ^r/(D„5A„)2 (3.90) 

a>i3 

can be written as: 

Y^trj\Tr{x, p5 A)\' + trf\Tr{p-5 A)P 

= l^trf\P'-5Aj^ 

ah 

= -trfTr{5 A P P15 A) (3.91) 

where in the second hne we have used the completeness of the Oa over the 2x2 matrices. The 
trace, Tr, refers to the spin indices and we have freely integrated by parts. Therefore the partition 
function in the Feynman gauge: 

Z =lim [ 5A5C5H- exp{-SYM)e'^'^"''''^"P'^"f^''^''' 

Det{-A-A)e-^^'''"''f^^">d{D,,dAa-C) (3.92) 

becomes: 

e^P(-7 I ftrf{p-{D^jAp^)fd'x)eM-^ f d'xtrf{D„dAaf) 

= J 5^-exp{-^-^^-^ J trf{^^pfd'x)Deti-A-^)DerH- P pd,h) (3.93) 
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where: 



Hence: 



In addition: 



p- p■ = -^J,l + a^fiF„^ (3.94) 



Der^- p p-8,y)=Det-"\-Kj,5afi-iF+^) (3.95) 



- p p = -^■^\ + Oa|iF„p (3.96) 

and symmetrically: 

Det-'^{- p pd,i,)=Det-'/\-A-Ado,p-iF-) (3.97) 

Now, since p p and p p have the same spectrum of non-zero modes, it follows that —A-^Sap — 
iF^A has the same non-zero modes as —A-^Sap —iF^g, but it has no zero modes precisely when 
— A;^5„jg — iF^A has. This is the case when the gauge connection that solves the equation of AS D 
type, F^g(A(m;)) = jx^n, has moduli, m,. Indeed taking the derivative of this equation with respect 
to the moduli one gets: 



S^ap 5A 



r 



dAy dnti 



(3.98) 



that implies that the operator: 



^5A,= ^^^^5A, = P-(D[„5A^]) (3.99) 

has zero modes, and therefore p and {—A'^Sap —2/ ad{P^F)ap) have zero modes too. QED 

3.4 One-loop beta function of pure YM in the ASD variables 

We can now use Eq.(3.81) as the definition of the partition function of YM in the ASD vari- 
ables. We can apply the standard background field method of sect.(3.1) for the computation of 
the beta function in the ASD variables. The field iJ.„p = }Xap + SfXap can be decomposed in a 
background, jX^p , and a fluctuating field, SlXap ■ The correlations of the fluctuating field can con- 
tribute only starting from order of g^. Therefore the only 0(^") contributions, relevant for the 
one-loop beta function, arise from the functional determinants. To evaluate the effective action in 
the ASD variables it is most convenient to compare it with the standard one-loop effective action of 
sect.(3.1). In the standard background field method the quadratic form: 

^ j d^xtrf{F^pf (3.100) 

is expanded around a solution of the equation of motion, leading to the one-loop effective action: 

Zi-i„op =e-ri-w(A) 

= e-''"'^^^Det-'/\-A-^5ap-2iF^p)Det{-A-^) (3.101) 
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In the ASD variables because of the delta function that defines the resolution of identity, the "action" 

is expanded around the background F^n{A + 8 A) = jX^n + 8}X^„, leading to the one-loop effective 
action in the ASD variables: 

2ASD _ g-rn„,(A) 



i^-4/"-/(^)'^'- 



/ 



A'"'^^''^co'^Det-^l^{-^■^5afi-2iad^)Det{-^■£) (3.103) 



Thus the orbital contribution is the same as in the standard background field method. The difference 
is the spin term, 2iad{P^F)fj,p, as opposed to liadFap, and the possible contribution of the zero 
modes, whose occurrence is not generic but depends on the background. Let us suppose at first that 
zero modes do not occur. In sect. (3.1) we have seen that the combination of determinants 

Der^^^{-A-^5ap - 2i adFap)Det{-Aj^) 

(3.104) 

leads to the beta function: 

VM = V(Ay^^3(4^"(4^^^°^M ^^'^^^^ 

where the first term in the brackets is the orbital contribution and the second one is the spin contri- 
bution. 

In the ASD case the orbital part is the same one, while the spin part differs because of the 

^ap -^ -f ^ap — -^ 
ASD variables the orbital contribution is: 



substitution F^r —5- P F^^a = -^. Thus both in the standard one-loop effective action and in the 



On the contrary, while the spin contribution for the standard one-loop effective action is: 



-A^ 



in the ASD variables, because of the substitution F„^ — )■ P F^p = -|^, is: 



N I trf{p-F,pf— -log- (3.108) 



(47r)2 ^n 
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Hence 






3(4;r)2 ^^u {An) 
2> +^j^niP-Pa,f (3.109) 



where: 



2 g" A 

3 (4;r)2 ^^ ^ 






v-l 1 10 / , A ,aiim 

Thus generically the YM beta function is not reproduced in absence of zero modes as in the SUSY 
case. Moreover generically the renormalizations of Q and of lJ.(^p are different. However, if the 
background field satisfies the equation of motion at leading order, ix^p = 0, the corresponding YM 
connection is SD and therefore instantons occur. In this case the zero modes have to be included 
and the one-loop beta function is reproduced in the ASD variables: 



"A4^e^2ive (3.ni) 



(4;t)2weZgl 



Yet, we may wonder as to whether the YM theory can be exactly localized on instantons as in the 
SUSY case. It is very instructive to check that it cannot be so, otherwise the two-loop beta function 
is not reproduced. Rescaling fields in canonical form, as in the SUSY case, we get: 



1 

(4)t)2A'e{szjAc)Zgl 



zf-foop =e ^r^A^^e«(4A,)2^e 



1 

(4)l)2jve(gz|A£.)Z-l 

:2 



1 



A4^e(^zi)4we^(yi^)2ive 



(3.112) 



Now we can define as in Eq.(3.67) 



m6n^Q(gZU,)Z^^ Nl6n^Q(gZU,)Z^^ i i 
V = 2g^ +4JV6(gZ^.)log(gZ|) (3.113) 



that implies: 



4Zn ' 

+ — ^log(gZ^) (3.114) 



2g^ 2g2 - (4nr ^'^ Q 
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and with two-loop accuracy 



ii^-v^F^"'^**"^' <'■"'' 



Taking the derivative with respect to log A and assuming by the localization hypothesis that gw is 
one-loop exact, we get with two-loop accuracy: 

' ^« - -ft + ^i^ + ^£i5M^ 0.116) 



g^^logA {4ny g dlogA (47r)2 5logA 



with j3o = T^- Therefore 



3 I 2g3 d log Zq 



dg _ h>Og -I- (4^ c)logA 



5 log A 



(3.117) 



(4;r)2 



Since 



it follows that: 



aiogZg 2 gl, 2 g 



a log A 3(471)2 3(471) 



11 ^3 4 ^5 4 „ 5 11 g^ 48 g^ 



(3.118) 



aiogA ' 3(471)2 3(471)4 (471)2^°^ 3(471)2 3(471)4 ^^-^^^^ 

Therefore the second coefficient of the beta function, j8i, differs from the perturbative result: 

48 1 34 1 

Thus it is not possible to localize the YM partition function on instantons. QED 

On the contrary, we will see in the following sections that twistor Wilson loops can be localized 
on surface operators with Zn holonomy. 

4. Twistor loops and non-commutative YM 
4.1 Non-commutative Eguchi-Kawai reduction 



We recall some fundamental facts about the non-commutative YM theory [J78|, |100| ] that will 
be used throughout the whole paper. These results will allow us to construct the twistor Wilson 
loops which our approach is entirely based on. 

The non-commutative 7?'^ is defined by: 

[x'',xP] = ie"i^i (4.1) 

Let A(x) be: 

r rl'^k 
A{x)=J ^^e"^e-''^ (4.2) 
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and 

/= f d''xf{x)A{x) (4.3) 



for complex functions of rapid decrease in both coordinates and momenta (Schwartz space). This 
defines an operator/function correspondence such that: 

/! = /*§ (4.4) 

with: 

if*g)ix)=f{x)eM^dy''PdP)giy)\y=, (4.5) 

that can be extended to multiple • products: 

Mxi)*...*Uxn)=llexp{^d^^e''Pd^)Mxi)...Uxn) (4.6) 

needed in the evaluation of Wilson loops of the non-commutative theory in the function represen- 
tation: P^^expi J^ Ao,{x)dxa- 

By the operator/function correspondence translations are represented by unitary operators: 

e''^A{x)e-''^ =A{x + a) (4.7) 

where: 

d'{x^) = 5'h (4.8) 

Thus non-commutative derivations can be represented via Eq.(4.1) and satisfy: 

[da,dp]=^^ap^ (4-9) 

In addition the integration on functions coincides with the operator trace up to a factor: 

{2n)ipf{d)trf = jd'xf{x) 

{2%)ipf{e)fr{k{x)K{y)) =5\x-y) 

Jd'x{f*g){x) = J d'xf{x)g{x) (4.10) 

The YM action of the U{N) non-commutative gauge theory has the function/: 



2g 
/operator representation: 

N ,^ .d 



YjJd''xtrN{Fc,p*Fc,p){x) (4.11) 



-^{2n)-2Pf{e)trNTr{-i[da + iAa,dp+iAp] + d^^\y (4.12) 
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where the non-commutative gauge connection is valued in the tensor product of the Lie algebra, 
u{N), of U{N) in the fundamental representation and of the field ^-algebra. This leads to the 



non-commutative [ ]101| , |102| ] ^"^ Eguchi-Kawai reduction [ ]103| , |104 |105[ ] 

JLAr(^)rffr^rr^(_/[5„ + iA^^dp + iAp] + 0„^ 1)2 (4.13) 

where the trace Trj^ is taken now over a subspace of dimension A^, with 

N{^Y = {27i)ipf{d) (4.14) 

in the large N,d,A limit. The simplest way to understand the occurrence of the inverse power of 
the cutoff ^^ in the reduced non-commutative action ^^ is to study the Makeenko-Migdal [Q, ^ 
loop equation after having reabsorbed the two factors oi N,N into a unique factor, ^ = NN, that 
computes the rank of the tensor product. For this we need to write the Wilson loop of the non- 
commutative theory in the operator notation. In this version the theory is a matrix model of infinite 
matrices. Thus the Wilson loop must involve a connection constant in space-time: 



trNTrf/'¥{A;L,,^,) = —trNTrf/Pexp / {da + iAa)dxa (4.15) 



1 _ .,.,-> - . 1 



Indeed this prescription leads to the correct definition of the Wilson loops of the non-commutative 
theory in the function representation. The proof is as follows. In the operator representation we 
can gauge away the non-commutative derivative that occurs in the definition of the Wilson loop 
by performing a local gauge transformation with values in the infinite-dimensional unitary group 
acting on the Fock representation of the non-commutative theory: 

(/(x) = e^«'^« (4.16) 

where x^ is a commutative space-time coordinate. The operator-valued gauge connection trans- 
forms under this gauge transformation in the usual way: 

A^ = (/(x)A„t/(x)-' +/5„f/Wf/(x)-i (4.17) 

where the partial derivative is the usual partial derivation with respect to the commutative param- 
eter x«. The operator a a must instead transform as a Higgs field in order for —id a +Aa to be a 
connection: 

dl = l]{x)dal]{x)-^ (4.18) 

Correspondingly the Wilson line: 

^'{A;Ly,)=PexpiJ (-/5„+A„)Jx„ (4.19) 

We would like to thank Antonio Gonzalez- Arroyo and Chris Korthals-Altes for discussions on non-commutative 
Eguchi-Kawai reduction and Antonio Gonzalez- Arroyo for a detailed exam of our work at the GGI. 
■'^We would like to thank Yuri Makeenko for discussing this point with us at the GGI. 
See | p7| ] for a modem treatment. Another interesting way to understand the same factor in the quenched version 



of the Eguchi-Kawai reduction is in [ 106 ] 
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transforms as ^': 

JLy- 

= Pexpi f {U{x)AalJ{xY^ -iU{x)daU{xy^+idaU{x)U{xy^)dXa 

JLy- 

= Pexpi U{x)AaU{xy^dxa 

JLy- 

= P^expi Aa{x)dxa (4.20) 

JLy- 

that is the function version of the non-commutative Wilson loop by the operator/function corre- 
spondence. Now the Makeenko-Migdal loop equation [^, |3|] of the large-A'^ commutative theory 
is: 

A^ dAa{x) 
= i [ dya5^'\x-y) < ^trM'i'{A;L^-) >< ^frw*P(A;L,,) > (4.21) 

where the normalized commutative YM action is: 

Sym = -^ j trN{Fajifd''x (4.22) 

and the v.e.v. is defined with respect of the unnormalized action: 

< ... >=Z-' J ...cxp{-^ J trN{F„pfd''x)5A (4.23) 



The loop equation of the non-commutative matrix model is instead M, 107]: 



= ij dya < -^Trj.^'iA-L^.) X -LTr^yi^{A;Ly,) > (4.24) 

where the normalized action, Snc, of the non-commutative theory is: 

1 27r 

SNC = Yl(^yTr_^i-ii^« + '^«^h+^M + ^ap^f (4.25) 

and the v.e.v. of the non-commutative theory is defined with respect to the unnormalized action: 

< - >=Z-' J ...cxp{-^{^y'Tr^{-i[da + iAa,dp+iAp] + d^^\)^)5A (4.26) 

At this point we notice that the factor of (^)'^ in the normalized non-commutative action is essen- 
tial to reproduce the loop equation of the commutative gauge theory, since its effect is equivalent 



We ignore central terms that vanish for large 6. 
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to the insertion of the missing 5^^'' (0) in the right hand side of the non-commutative loop equation. 
As a consequence the d^'^\x — y) of the commutative loop equation is reproduced provided the 
trace of Wilson lines vanishes for xj^y [^ , 107 ] : 



<—Trj.W{A;L^)>=0 (4.27) 

QED 

The occurrence of the inverse power of the cutoff in the matrix model version of the non- 
commutative theory opens the way to saddle-point computations of new kind in which power- 
like divergences cancel against the {^Y factor In particular if the theory is defined on R^ x R^ 
quadratic divergences cancel. This will turn out to be the case for the surface operators of the theory 
on R^ X /?0 in the limit of large 6 introduced in sect. (7) and employed in the whole paper 

4.2 Twistor Wilson loops 

We define twistor Wilson loops in the YM theory with gauge group U (N) on R^ x R^ with 
complex coordinates {z = XQ + ix\,z = XQ — ix\,u = X2+ixT,,u = X2 — ix^ ) and non-commutative 
parameter 6, satisfying [du,da] = d^^l, as follows: 

Tr^A^^'{B^;L„„) = TryPexpi f {A, + XD,)dz + {A=, + X-^Da)dz (4.28) 

where £>„ = d^ + iA^ is the co variant derivative along the non-commutative direction m and A a 
complex parameter For many purposes it is not restrictive to choose A real, although other choices 
are possible, for example a phase, A = e'^. The plane {z,z) is commutative. The loop, L„„, starts 
and ends at the marked point, w, and lies in the commutative plane. Thus we regard the twistor 
connection, B^, whose holonomy the twistor Wilson loop computes, as a non-Hermitian connection 
in the commutative plane valued in the tensor product of the U {N) Lie algebra and of the infinite- 
dimensional operators that generate the Fock representation of the non-commutative plane {u,u). 
Bx is indeed a connection in the commutative plane since the non-commutative covariant derivative 
transforms as a Higgs field of the commutative plane. The trace is defined accordingly. The limit of 
infinite non-commutativity in the plane (m,m) is understood, being equivalent to the \wge.-N limit 



of the commutative gauge theory |101, 78, 108]. The U{N) non-commutative theory for finite 
has tachyon instabilities that occur in non-planar diagrams suppressed by powers of 0^' and A'^^' 



[J78p. Therefore non-commutativity is for us just a mean to define the \wgQ-N limit, as well as it is 
for Nekrasov just a mean to compactify the moduli space of instantons ^^. 

4.3 Fiber independence of the v.e.v. of twistor Wilson loops 

It easy to show that the v.e.v. of the twistor Wilson loops is independent on the parameter A: 

< ^Tr._,v'i'{Bx\L^'.) >=< ^Tr.A^'i'iBuK^) > (4.29) 



Once localization is obtained, the glueball spectrum is computed employing the effective action in the large-A' 
commutative theory, around the localized locus (sect. (12)). 
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The proof is obtained changing variables, rescaUng covariant derivatives in the usual definition of 
the functional integral of the non-commutative YM theory: 



exp{-^{^f J d^xTrj.{-i[Da,Dp] + d^'p\f)dAdA6Db 

= [ Trj^Pcxp i [ (A, + D',)dz + (A,- + D',)dz 

cM-^{^f J dhTrM-[D'a,Dfpf + {e„'p) l-2i[D'„,D'p]d„'p)5A5A5D'D' (430) 



where: 



Z 


= D, 


Z 


= D-z 


D'u 


= xbu 


D'u 


= x-'b. 



(4.31) 

The formal non-commutative integration measure is invariant under such rescaling because of the 
pairwise cancellation of the powers of A and A ^ ^ The first term in the non-commutative YM action, 
proportional to Tr^ [D'q,,Do]^, is invariant because of rotational invariance in the non-commutative 
plane. 

Indeed every u must be contracted with a m by rotational invariance in the non-commutative 
plane and thus the factors of A cancel. The only possibly dangerous terms couple the non-commutative 
parameter to the commutator Tr_yy{[b'f^,b'„]d^l) but only rr^([Dj,,D^]6j^^') survives, because all 
the other terms are zero for /?^ x /?g. But the commutator is invariant under A -rescaling. 

We notice that, after rescaling, the integration variables (5D^, 8b'j^) should be treated as inde- 
pendent. For A real this is appropriate if we analytically continue the non-commutative plane to 
Minkowski space-time, after which the A invariance of the loop is simply invariance under Lorentz 
boosts ^^. The analytic continuation is also connected with the large-0 triviality (see below). QED 

In fact the twistor Wilson loops are trivially 1 at large-0 to all orders in the 't Hooft coupling 
constant g: 

lim < -\-Tr^^{B}_-K„) >= 1 (4.32) 

Firstly, we show that triviality holds to the lowest non-trivial order in perturbation theory. We have 
in the Feynman gauge in the large- limit ^*': 

< Try { [ (A, + Xbu)dz + (4- + X-^bu)dz j (A, + Xbu)dz + (A,- + X-^ba)dz) > 



■^^We would like to thank Konstantin Zarembo for discussing with us the A -independence at the GGI. 
We would like to thank Luca Lopez for working out a detailed version of this computation during our course at 

SNS. 
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= 2J dzj Jz(< rr^(AA-) > +/2 < rr^(AA) >) 

= (4.33) 

4.4 Twistor Wilson loops are supported on Lagrangian submanifolds of twistor space 

Secondly, we show that triviaUty holds in the large -limit to all orders of perturbation theory. 
For this aim it is convenient to gauge away the non-commutative derivatives that occur in the 
definition of twistor Wilson loops. This can be done by performing a local gauge transformation 
with values in the complexification of the gauge group. Although this is not a symmetry of the 
theory, the trace of the twistor Wilson loops is left invariant because of the cyclicity property of the 
trace. Let be 

5(z,z) = e'^^^"+'^"'^"^"- (4.34) 

where (z,z) are commutative coordinates. The components of the operator- valued gauge connec- 
tion, Bx, transform under this gauge transformation in the usual way: 

bI^ = S{z,z)BxJ{z,zr'+idJ{z,z)S{z,z)-' 

Bl- = Siz,z)BxAz,z)-'+id-Jiz,z)S{z,z)-' (4.35) 

where the partial derivatives are the usual partial derivations with respect to the commutative pa- 
rameters {z,z)- Correspondingly the twistor Wilson line transforms as '*': 

Siw,w)^iBx;L„,)S{v,v)-' 
= Pexpi f Bljlz + Bl-dz 

= Pexpi f (5(z,z)(4 + AD„)5(z,z)-i-A5„)Jz + (5(z,z)(4- + A-'D„-)5(z,z)-i-A-'5„-)Jz 

= Pexp/ / S{z,z){A, + i?iA^)S{z,zy^dz + S{z,z){Az + i^'^Aa)S{z,zy^dz (4.36) 

Therefore the twistor loop lies effectively on the submanifold of four-dimensional commutative 
space-time defined by: 

{z,z,u,u) = (z,z,/Az,/A"'z) (4.37) 



with tangent vector: 



(z,z,m,m) = (z,z,/Az,/A ^z) (4.38) 



This is a Lagrangian submanifold of the (complexified) Euclidean space with respect to the Kahler 
form dz/\dz + duA du that lifts to a Lagrangian submanifold of twistor space provided A is either 
real or a unitary phase. The two cases correspond to Lagrangian submanifolds of antipodal and 
circle type respectively. 



We ignore central terms that vanish for large 0. 
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4.5 Triviality of twistor Wilson loops in the limit of infinite non-commutativity 

The proof of triviality of twistor Wilson loops to all orders of perturbation theory in the limit 
— ;. oo follows now almost immediately. 

Indeed at any order in perturbation theory a generic contribution to an ordinary Wilson loop 
of a commutative gauge theory contains a correlator of gauge fields, i.e. a Green function, with 
tensor indices contracted with a product of monomials in Xa {s) at generic insertion points on the 
loop, labeled by s: 

dsids2 Gaia2...{^p{^l)-^li{^2),---)xai{s[)xa2{s2)--- (4.39) 

Because of the (9(4) invariance of the commutative theory Xa^ (si) is contracted either with another 
Xa2{^2) or with an Xa2{s2) to form polynomials in Xa{s)xa{s') or in Xa{s)xa{s'). Indeed all these 
monomials necessarily contain at least one factor of Xa since the gauge field along the loop has the 
index contracted with the one of i^- The possible factor of Xa arises from the dependence of the 
Green functions on the coordinates. 

We now specialize to twistor Wilson loops. 

In the limit — ;• oo of the non-commutative gauge theory 0(4) invariance is recovered, because 
the theory becomes the large-A^ limit of the commutative theory, that obviously is 0(4) invariant. 

Therefore all the monomials just mentioned vanish when evaluated on the Lagrangian subman- 
ifold which the twistor Wilson loop lies on, because they are of the form z{s)z{s') — z{s)z{s') = 
or z{s)z{s') — z{s)z{s') = 0. Thus the "effective propagators" that connect a Feynman graph at any 
order to the twistor Wilson loop vanish. The only factors that may spoil the triviality occur if sin- 
gularities due to denominators of Feynman diagrams arise, since Xa{s)xa{s') vanishes too on the 
Lagrangian submanifold for the same reasons. To cure this we analytically continue the correlators 
that occur in the computation of twistor Wilson loops from Euclidean to Minkowski space-time in 
order to get the ie prescription, z+ {s)z- {s') — z+ {s)z- {s') + ie = ie in the denominators. 

The gauge invariant prescription of analytic continuation from Euclidean to Minkowski space- 
time will be used over and over in the paper and it will play a crucial role. QED 

We describe now the aforementioned analytic continuation of the twistor Wilson loops at the 
operator level in the functional integral by means of the following sequence. 

Firstly, we analytically continue to Minkowski space-time only the commutative plane. Then 
at operator level the twistor Wilson loops become: 

Tr^^>{Bx;Kw) -^ Tr^Pexpi f (A^+ + iXDu)dz+ + (A,_ + iX-^Da)dz-^ (4.40) 

since xo—>/x4,z—>Jz+ and z —> jz_ , with (z+ =X4+x\,Z- =X4—x\) andA^— ^ — /Aj.^,A|— ^ — i^z_- 
The support of the twistor Wilson loops analytically continued in this way becomes: 

{z,z,i?iz,a-^z) -^ {z+,Z-,-?iz+,-?i-\-) (4.41) 

that is Lagrangian with respect to —dz+ A dz- + du f\du for a real section of the complexified 
Euclidean space-time. 

In sect.(6) we write a holomorphic loop equation that the twistor Wilson loops satisfy in Eu- 
clidean space-time. The holomorphic loop equation involves in the left hand side an Euclidean 
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effective action that should be renormalized in Euclidean space-time and in the right hand side a 
contour integral along the loop that is not well defined in Euclidean space-time and that should be 
regularized. 

We will show that there is an essentially unique way of regularizing by analytical continuation 
to Minkowski space-time. Thus after renormalization of the effective action the holomorphic loop 
equation makes sense in Minkowski space-time. In sect.(ll) the renormalized effective action in 
Minkowski space-time restricted to fluctuations of surface operators supported on the Lagrangian 
submanifold: 

(z+,z_,-Az+,-A-iz_) (4.42) 

is used to compute the glueball spectrum. In the effective action this Lagrangian submanifold is 
obtained first restricting to the Lagrangian submanifold in Euclidean space: 

(z,z,-Az,-A-iz) (4.43) 

and then analytically continuing. This can be done in two ways. The following choice for the 
analytic continuation u — )• iu^,U — )• /m_ leads to trivial twistor Wilson loops: 

Tr^^{B^;L„„)^Tr^Pexpi f {A,^ + XD,Jdz+ + iA,_+ X-'Du_)dz- (4.44) 

supported on the Lagrangian submanifold: 

{z+,Z-,az+,iX-'z-) (4.45) 

for which the correlators that occur in the evaluation of the twistor Wilson loops are vanishing as 
shown in the triviality proof. The other choice for the analytic continuation m — )■ m+, m — >• m_ leads 
to: 

Tr^^>{Bx;Kw) -^ Tr^Pexpi j {A,^ + iXD„^)dz+ + (A,_ +a"'D„JJz_ (4.46) 

supported on the Lagrangian submanifold: 

(z+,z_,-Az+,-A-'z_) (4.47) 

and to non-trivial twistor loops that satisfy the same loop equation in Minkowski, but not the lo- 
calization property. The point is that the effective action, F, is naturally defined on the Lagrangian 
submanifold that is the support of the non-trivial twistor Wilson loops, and thus it carries the inter- 
esting information of the localization. 

The trivial twistor loops are obviously finite at = oo, i.e. they have no cusp and perimeter 
divergences, in analogy with certain supersymmetric Wilson loops. 

Indeed the cognoscenti may have noticed that twistor Wilson loops resemble locally BPS Wil- 
son loops of theories with extended supersymmetry. In fact our triviality proof mimics the argument 
about a certain non-renormalization property [|T^] of locally BPS Wilson loops. Indeed it has been 
argued in [13] that a locally BPS Wilson loop in the four-dimensional c/K = 4 SUSY gauge theory: 

TrPexpi AaXa{s)ds + i(j)j,yh{s)ds (4.48) 
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has no perimeter divergence, to all orders in perturbation theory, because of the local BPS con- 
straint: 

I4(^)-IjI(^)=0 (4.49) 

a h 

At lowest order of perturbation theory this constraint assures the cancellation of the contribution 
to the perimeter divergence of the gauge propagator versus the scalar propagator, because of the 
factor of f- in front of the scalar propagator at that order. 



As far as the perimeter divergence is concerned, it is argued in 1 13] that this cancellation occurs 
to all orders in perturbation theory, when the locally BPS Wilson loop is seen as the dimensional 
reduction to four dimensions of the ten-dimensional Wilson loop of the ten-dimensional J^ = 1 
SUSY YM theory from which the four-dimensional ^ = 4 SUSY theory is obtained. 



Remarkably, in the argument of [13] SUSY plays no direct role. In fact the argument is based 
only on 0(10) rotational symmetry of the parent d = 10 .yV = 1 SUSY gauge theory from which 
the daughter d = 4 JV = 4 SUSY gauge theory derives by dimensional reduction, as we show 
momentarily. 

In ten dimensions the coefficient of the perimeter divergence of an ordinary unitary Wilson 
loop, at any order in perturbation theory, must necessarily contain as a factor a polynomial in the 
0(10) invariant quantity Y,m^m{^)- Indeed the perimeter divergence arises when all insertion points 
coincide, in such a way that all the arguments of the Green function vanish. 

In this case the Green function provides a factor that, by 0(10) rotational invariance, must 
be a polynomial in ten-dimensional Kronecker delta, since all the difference vectors in the Green 
function are zero at coinciding points and thus no other tensorial structure can be produced. 

This combines with the factors of xm{s) to produce an invariant polynomial in Lm ■*!/(■*) with 
no constant term, since the lowest order contribution is zero by direct computation. But Y^M^lti^) — 
Ha^ai^) ~ Lfc jfc('^) = is zero for a BPS Wilson loop because of the BPS constraint. A naive 
application of this argument to the four-dimensional BPS Wilson loop would imply the absence of 
the perimeter divergence for this loop on the basis of the 0(10) rotational invariance of the theory 
before the dimensional reduction. QED 

5. The quasi-localization lemma for twistor loops in large-A'^ YM 

We use the A -independence to show that the v.e.v. of twistor Wilson loops is localized on the 
sheaves fixed by the semigroup rescaling A . This involves a delicate and subtle interchange between 
limit and integration, that will be justified in sect.(8), after introducing a lattice regularization of 
the functional integral of differential geometric nature in sect.(7). In addition it will be checked by 
direct computation that the result of the localization agrees with the perturbative triviality. 

In this section for simplicity we use a notation that does not distinguish between commutative 
and non-commutative theories and therefore we do not add hats to operator valued quantities of 
non-commutative theories. The framework has been set in the previous section, therefore this use 
should not generate ambiguities. 

It is convenient to choose our twistor Wilson loops in the adjoint representation and to use the 
fact that in the laige-N limit their v.e.v. factorizes in the product of the v.e.v. of the fundamental 
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representation and of its conjugate. Then, for the factor in the fundamental representation, local- 
ization proceeds as follows. We write the YM partition function by means of the non-^L'^^F analog 
[^ of the Nicolai map of ^ = 1 SUSY YM theory worked out in sect.(3.3), introducing in the 
functional integral the appropriate resolution of identity: 



^= S{F^p-^ap)S^ap (5-1) 



exp(--^G--3 I / Trf{^i„j)d'x)d{F^p-^^p)5n^p5Aa (5.2) 

Q is the second Chem class (the topological charge) and /i^g is a field of A5D type. The equations 
of ASD type in the resolution of identity, Fqi — F23 = Moi , ^02 — ^ 3 1 = M02 1 ^03 — ^12 = M03 , can be 
rewritten in the form of a Hitchin system (taking into account the central extension that occurs in 
the non-commutative case): 

-iFa + [D,D] - 0-1 1 = Ai° = -Aioi 

-idAD =n = - (Ai(72 + '^03 ) 

-BaD =n = - (^02 - 'Mo3 ) (5.3) 

or equivalently in terms of the non-Hermitian connection whose holonomy is computed by the 
twistor Wilson loop with parameter p,Bp = A + pD + p^^D = {A, + pDu)dz + (A=; + p^^Dujdz, 

-iFBp-B^h =iip=ii'^ + p^^n-pn 
—idAD = n 
-idAD =n (5.4) 

The resolution of identity in the functional integral then reads: 

\=J5n8nj 5lXp5{-iFBp-llp-d-h)5{-idAD-n)5{-idAD-n) (5.5) 

J •J Cp 

where the measure, 5/ip, along the path, Cp, is over the non-Hermitian path with fixed n and h 
and varying pL^. The resolution of identity is independent, as p varies, on the complex path of 
integration, Cp. Indeed the constraint implied by the first of Eq.(5.4) on Cp is a linear complex 
combination of the Hermitian constraint in the first of Eq.(5.3) with coefficient 1 and of the re- 
maining two with coefficients depending on p . 

Let us consider the v.e.v. of twistor Wilson loops: 

5n5n / 5^pexp{-—^Q--j Trf{^y +4Trf{nn)dS) 

TrfPexpi j {A, + XD„)dz+{Ai + X-^Da)dz 

5{-iFBp -jJLp- e-^ l)5{-idAD - n)6{-idAD - h)dAdAdD5D (5.6) 
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and let us change variables in the functional integral rescaling the non-commutative covariant 
derivatives: 

fdndnf dHpexp{-^^^Q-^ f Trfin'-^f +4Trf{nn)d'^x) 

TrfP exp / / (A, + D',)dz + (A^ + D'a)dz 

5(-iFA + [D',D']-e-h-n°-i-dAD' + i^dAD'-p-^n + pn) 

P A 

8{-adADf -n)8{-a-^~dAD' -n)8A8A8D'8D' (5.7) 

Taking the limit A — ;■ inside the functional integral, the last line implies localization on « = and 
BaD' = 0. The 8n integral is performed by means of the second delta function. 8h decouples from 
the third delta function. The independence on the path Cp in the neighborhood of p = 0, that we 
denote, choosing p = aX, Cax+ with a fixed as A — ;■ 0+, implies that the 8n integral decouples 
and that BaD' = as well. 

Firstly, on Cax+ the integral over 8h decouples because 8h disappears from the first delta 
function as well. Secondly, the first delta function implies « = as A — ;■ as well. Therefore, setting 
« = from the start in the first delta function because of the second delta function, the argument of 
the first delta function contains the combination —IFa + [D' ,D'] — Q^^l — /x" + ia^^dAD', that can 
be a independent and thus zero for every a only if the two terms are zero separately. Therefore 
also BaD' = 0. Notice that for this argument to hold it is not necessary to consider the variables 
{n,D') as Hermitian conjugated to {n,D'). Indeed they are not so because of the A rescaling and/or 
the analytic continuation that is necessary to regularize the twistor loops for the triviality property 
to hold. 

We notice that the localized density has a holomorphic ambiguity, since we can represent the 
same measure using a different density, performing holomorphic transformations without spoil- 
ing the quasi-localization lemma : 5/Xo+ = -^^^81x1^+. The final result for the localized effective 
measure is: 



^+ 



[/ 5Mo.|f^exp(-^e-£,l/rr,(M„>^x)5(F-,-M„,)]„^,^^^^^ (5.8) 

where we have reintroduced the covariant notation. 

Thus the twistor loops are localized on the fixed sheaves for which two of the ASD fields 
vanish: jUqj = ^03 = 0. QED 

The integration on the gauge connection in Eq.(5.8) can be explicitly performed in the Feyn- 
man gauge in the way explained in sect.(3.3), to obtain: 

Ci 

= [/ 5M'exp(-^e-^ I I Trf{n-4)d^x) 

Der2(-A;i5„/3 - iad^-jDet{-hA){:^)'"'Det'2co^, x c.c]^^.^^ (5.9) 
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The complex conjugate factor arises by the conjugate representation, Det^^CO is the contribution of 
the Jib zero modes due to the moduli, and A the corresponding Pauli-Villars regulator. 

The volume form on the connections admits several different representations as the Liouville 
measure associated to a symplectic form, since different symplectic forms may lead to the same 
volume form. One possible choice for the symplectic form, ft), is coj, displayed in Eq.(3.82). 
However, in sect.(7) and sect.(12) we will see that there is a different choice of co, compatible with 
holomorphic/antiholomorphic fusion, that is most convenient in the computation of the glueball 
spectrum. 

6. Holomorphic loop equation for twistor Wilson loops 

6.1 Holomorphic loop equation 

We now specialize to the case p = 1 for the twistor connection B\=B. The partition function 
reads: 

Z = fdndnf 5Ai'|^exp(-^^(2-^ f Trji^^f +4Trf{nn)d\) 
J JC[ OfX g 8 J 

5{-iFB-H-e-h)5{-idAD-n)5{-idAD-n)5A5A5D5D (6.1) 

Writing the holomorphic loop equation for twistor Wilson loops requires that /i' be chosen in 
the holomorphic gauge, B^ = 0. Strictly speaking this is a change of variables and not a gauge 
transformation with values in the unitary gauge group, since this gauge can be reached only by a 
gauge transformation with values in the complexification of the gauge group. 

The further change of variables to the holomorphic gauge is a new key feature of our approach 
to the large-A'^ YM theory. It is based on the idea that twistor Wilson loops, being holomorphic 
functionals of jJ.', behave as the chiral (i.e. holomorphic) super-fields of an ^ = 1 SUSY gauge 
theory. In fact the new holomorphic loop equation resembles for the cognoscenti the holomorphic 
loop equation that occurs in the Dijkgraaf-Vafa theory [|55, 56, 57, 58] of the glueball superpotential 



^^ in ^ = 1 SUSY gauge theories. Thus the v.e.v. is taken with respect to the measure: 

dndn / d^'...exp{ —Q ^ / Trf{np.) + Trf{n + nfd'^x) 

5(-iFB-^-d~h)5(-idAD-n)5(-idAD-n)-^5A5A5D5D (6.2) 

The holomorphic loop equation is obtained following the Makeenko-Migdal technique, as an iden- 
tity that expresses the fact that the functional integral of a functional derivative vanishes: 



/ 



^"5^)^'^^'^^^'''^''^^^^' = ^ ^^-^^ 



The new holomorphic loop equation for twistor loops follows: 

oix'{z,z) kJl,,z-w 



(6.4) 



The name is perhaps misleading because the glueball superpotential cannot be used to compute the mass of any 
glueball state of the .yV = 1 SUSY theories. 
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where *P(B';L^^) is the holonomy of B in the gauge B'^ = 0. The Cauchy kernel arises as the kernel 
of the operator d^^ that occurs by functionally differentiating ^{B';L^^) with respect to }i' via 
idB'^ = }i' . Assuming that the loop L^^ is simple, i.e. it has no self-intersections, the holomorphic 
loop equation linearizes: 

< Tr{^-^^'i>{B';L,,)) >= - f -^ < Tr^{B';L,„) >< Tr\ > (6.5) 

6.2 Regularization by analytic continuation to Minkowski space-time 

The contour integration in the right hand side (i.e in the term that accounts for quantum fluc- 
tuations) of the loop equation includes the pole of the Cauchy kernel. We need therefore a regular- 
ization. 

The natural choice consists in analytically continuing the loop equation from Euclidean to 
Minkowski space-time, z —5- i{z+ + i£)- It is at the heart of the Euclidean approach to quantum field 
theory that this analytic continuation be in fact possible. 

In the approach to localization by the holomorphic loop equation the analytic continuation is 
performed only after functional integration and renormalization, that are performed in Euclidean 
space. Thus we think that this procedure has chances to work also from the point of view of the 
constructive quantum field theory. 

In fact the approach to localization via the holomorphic loop equation, when combined with 
the integration on local systems of the next section, leads to more complete and satisfactory results 
than the localization on fixed points. 

The result of the ie regularization of the Cauchy kernel is the sum of two distributions, the 
principal part of the real Cauchy kernel and a one-dimensional delta function: 

in8{z+ — w+) (6.6) 



Z+-W+ + ie z+ - w^ 
The loop equation thus regularized is: 



1 [ (p ^^+ _ i7tdw+5iz+-w+)) < Tr^>(B';L,^„^) >< Tr^{B'-K^,^) > (6.7) 
njL,.,, Z+-W+ 



1 f ,^ dw^ 

Z+-W+ 

where now both the distributions on the right hand side are integrable along the loop. This reg- 
ularization has the great virtue of being manifestly gauge invariant, an unusual feature for loop 
equations. In addition this regularization is not loop dependent. 

The right hand side of the loop equation contains now two contributions. A delta-like one 
dimensional contact term, that is supported on closed loops and a principal part distribution that is 
supported on open loops. Since by gauge invariance it is consistent to assume that the expectation 
value of open loops vanishes, as in Eq.(4.27), the principal part does not contribute and the loop 
equation in the holomorphic gauge reduces to: 

= -if dw+3in--w+)<T,<i'{B':L,^.J><T,'i'(B':L.,^,J> (6.8) 
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As briefly outlined in sect.(l.l), the holomorphic loop equation for cusped loops with the shape of 
the symbol oo and the cusp at the non-trivial self-intersection point reduces to a critical equation for 
the effective action F in Minkowski signature: 

since the contour integral in the right hand side of the holomorphic loop equation vanishes because 
of the opposite orientation in a neighborhood of the backtracking cusp. The geometrical side of 
this localization of the loop equation is described in sect.(lO). 

Taking w+ = z+ and using the transformation properties of n' and of the holonomy of B', the 
preceding equation can be rewritten in terms the curvature, jJ., and of the connection, B, in a unitary 
gauge: 



<Tr{ s:..,. . , 'i'{B;L,^,,))> 



5r_ 

Sl^{z+,z- 
ij dw+5iz+-w+) < Tr^>{B;L,^„^) >< Tr^>{B;K^,^) > 



(6.10) 



^z+z+ 



where we have used the condition that v.e.v. of the trace of open loops vanishes to substitute the 
holonomy of B' with the holonomy of B. 

As a consequence of the localization of the holomorphic loop equation, the equation of motion 
of r vanishes when restricted to the subalgebra of twistor Wilson loops: 

< Tr{nB-M,,Jj^^^^nB-,L'U)) >=0 (6-11) 

where L' and L" are the two petals of the loop. In the next section we construct an explicit realiza- 
tion of this sublagebra, on which we can interpret the localization of the holomorphic loop equation 
strongly: 

= (6.12) 



5iJ.{z+,z- 
7. Integrating on surface operators in large-A^ YM 

7.1 Integrating on infinite-dimensional local systems 

We would like to give a precise mathematical meaning to the formal manipulations of the func- 
tional measure in sect.(5) and sect.(6). One possibility would be introducing a lattice regularization 



of the functional integral according to Wilson [ |109| ]. However, this kind of lattice regularization 



would spoil completely the geometrical structure, since in the Wilson regularization the gauge con- 
nection lives on links and the curvature on plaquettes, a fact that makes exploiting the map from the 
connection to the ASD curvature problematic, not to mention the understanding of the moduli of 
the loci at which this map is not one-to-one, for which the zero modes necessary to get the correct 
beta function occur (sect.(3.3) and sect.(3.4)) . 

Therefore we introduce a new regularization of the YM functional integral that allows us to 
maintain the differential geometric structure. The differential geometric structure is crucial to get 
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a structure theory of the locus of the fixed points of the functional measure and to understand the 
zero modes of the determinants, that in turn affect the beta function of the theory. 

Our new regularization of the YM theory in the \aige-N limit is performed in two steps. In the 
first step the resolution of identity in the Nicolai map on R^ x /?g is represented in the operator no- 
tation of sect.(4. 1) as a functional integral on infinite-dimensional parabolic bundles, as suggested 



long ago in [22, 23] 



1= f 5{-i[Da,Dp]- -l^fi^^{p)5^'\z-zp)-e^yi)l\5fi^p{p) (7.1) 

•' p p 

In this notation all the dependence on the non-commutative coordinates is absorbed into the infinite 
dimensional nature of the operators that occur in the non-commutative Eguchi-Kawai reduction. 
Therefore the base of the infinite-dimensional parabolic bundles is the two-dimensional surface, 
R^, labelled by the commutative coordinates (z,z). 

This amounts to substitute the continuous field, jl^g{z,z), of the Nicolai map with the lattice 
field, fi„a{p), by the resolution: 

A„/3(Z'Z) =lA«/3(/')5^''(z-Z/') (7-2) 

p 

This resolution is dense in the sense of distributions, since for any smooth test function of compact 
support: 

ND'Lfizp,Zp)fi„p{p) ^ f f{z,z)fi^p{z,z)d\ (7.3) 

p ■'' 

On this dense set in function space ^^ the resolution of identity of the Nicolai map can be interpreted 



as hyper-Kahler reduction [g2[ g3p. Indeed the three constraints oiASD type are the Hermitian and 
the complex moment maps for the Hamiltonian action of the infinite-dimensional unitary gauge 
group on the commutative plane R^: 

p 

-idp = Y^h,5^^){z-Zp) 
p 

-idp = Y^hp8^^\z-Zp) (7.4) 

p 

with respect to the three symplectic forms [^, ^, |6|] '^'*: 

(Oj =— f (fztrffr{dA^ A 5A= + 5D„ A 5ba) 

The three operators ^i-„a{p) are reduced in fact to large finite dimensional matrices by Merita equivalence as 
explained momentarily. In the finite dimensional case these matrices all commute as a consequence of the local model 
of the Hitchin equations [BSL fed]. Therefore the resolution in Eq.(7.2) turns out to be dense in function space in the 
sense of the distributions only in a certain neighborhood of the fixed points, AoaCp) ~ f^mip) — 0- The missing degrees 
of freedom occur as moduli (see below). Nevertheless the local degrees of freedom and the moduli in a neighborhood 
of the fixed points are enough to reproduce the correct universal one- and two- loop contributions to the beta function 
(sect.(9) and sect.(ll)) . Hence the mentioned degrees of freedom are in fact dense in function space in a neighborhood 
of the fixed points in the large-A' limit (see below). 

We use the same labels (/, /, K) of the symplectic forms as in |En] for the finite dimensional case. 
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(Oj - icOK = 2^ / d^Ztrffr{5A, A 5D„-) 



(Oj + i(OK =— I d^ztrffr{dA,A5Di,) (7.5) 



as it follows immediately from the interpretation as (infinite-dimensional) Hitchin systems [52,p3|]. 

Thus we come to the remarkable conclusion that, on a dense set associated to a lattice divisor, 

the Nicolai map can be interpreted, up to gauge equivalence, as a resolution of the gauge connection 

into orbits parametrized formally by hyper-Kahler moduli spaces, that arise as the quotient of 



the manifold defined by Eq.(7.4) for the action of the unitary gauge group [52, 53]. Yet, for the 
moment, this construction is somehow formal, because of the infinite-dimensional nature of the 
bundles involved. Thus it is unclear what the moduli theory of these infinite-dimensional bundles 
is. 

7.2 Reducing to finite dimension by Morita duality and inductive structure 

In order to reduce to finite dimensional bundles a possible way out is to compactify the non- 
commutative plane, 7?^, on a non-commutative torus of large area, L^. The corresponding non- 
commutative U{N) gauge theory enjoys, for rational values of the dimensionless non-commutative 



'-2 _ M 



parameter, IkBLT- = ^, Morita duality [[78|, |108| , |110| , IlllQ to a theory on a commutative torus 
of area l?N^^, with gauge group U{N xN), with the same 't Hooft coupling constant g, and with 
twisted boundary conditions corresponding to a 't Hooft flux. 

Indeed, starting from the infinite-dimensional non-commutative case we can perform an infinite- 
dimensional unitary gauge transformation, U{u,u), depending on the commutative parameters, 
(m, m), in such a way that: 

-i[D^,b^]--Y^fi- {p)5^^\z-Zp)-^ail=0 (7.6) 



becomes: 



[daAp{z.,Z,u,u) - d^Aa{z,Z,u,u) +i[Aa{z,Z,u,u),A^{z,z,u,u)]\ 

= Y,U{u,u)^^^{p)JJ-\u,u)8^^\z-Zp) + Ki^ {1.1) 



where: 



A^{z,z,u,u) = U{u,u)A^{z,z)U {u,u) + idpU{u,u)U {u,u) (7.8) 

and ^aj} is an antisymmetric field that takes into account a possible U{1) background flux in addi- 
tion to the central term that arises by the inverse of the non-commutativity. 

The content of Morita duality is that these seemingly infinite-dimensional equations admit 



a finite dimensional solution [78] on a commutative torus with coordinates (m,m), gauge group 
U{N xN), twisted boundary conditions corresponding to a 't Hooft flux and a [/(I) background 
flux, 0^o. Despite the structure group isU{N xN), the Morita equivalent theory does not describe 
the most general U{N x N) gauge theory, but only the one that satisfies the twisted boundary 
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conditions. In particular the twisted boundary conditions require solutions of the kind: 

[daAp{z,z,u,a)-dpAaiz,z,u,a) + i[Aaiz,z,u,a),Ap{z,z,u,u)]] 

= l^U{u,u)^^p{p)U-\u,u)5^^Hz-Zp) + ^'„pl^^l^ (7.9) 

p 

where now U{u,u) are U {N) matrices and Ix^Ap) lives in the tensor product of Lie algebras, 
u{N) X u{N), [ JTSl , |108 . |110| , |111[ ]. These are the equations that define surface operators in finite 



dimension but for the fact that the gauge connection satisfies twisted boundary conditions on the 
torus. 

More explicitly any U{N) connection of the YM theory on R^ xTg, with coordinates {y,x), 
with periodic boundary conditions on the non-commutative torus, Tq, with 2ndL^^ = j-, admits 



N' 



the expansion [108]: 



A{x,y) = £ a,{y)e-^''"'/'^ (7.10) 

/ez2 



and the corresponding Morita equivalent U{NN) connection reads [ |108| ]: 



A'{x,y) = ^ a;(3;)y-'>^''^'2^-^'''2/2e-2;r,7-^/*^ (7.11) 

/6Z2 

where the matrices {U,V) are the clock and shift matrices of SU{N) and (O = e^"^'^'^: 

UV = (oVU (7.12) 

The traceless part of A' is a connection on the twisted 't Hooft bundle SU{NN)/Zj^fj, with magnetic 
fluxM' = rNmod{NN) [ M : 



A'{xj + L/N) = TjA'{xj)Tj (7.13) 

with 

Ti = 1a, X V 

r2 = lwxV (7.14) 

where r is an integer that occurs in the definition of the SL{2,Z) matrix that defines the Morita 
equivalence [ |108| ]. 

Therefore the theory contains an untwisted sector SU{N) x 1^ that is diagonally embedded in 
U{N xN). This untwisted sector plays a special role, because carries zero momentum on the torus 
and thus it may condense. 

We choose in this paper ^^ = ^ =— )• ^ for any integer n > 2, because it allows us to per- 
form the large 6 limit, necessary to reproduce the \aige-N limit of the commutative SU (N) theory, 
uniformly for large N. It turns out that the choice of n is just the choice of a renormalization 
scheme associated to a particular inductive sequence of finite dimensional bundles used to define 
the large-^ hmit (sect.(12)). 

From a physical point of view we apply 't Hooft duality ideas to the U{N) non-commutative 
theory. The structure group of this theory is the tensor product of L'^(A'^) and of the group of *-gauge 
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transformations. We assume Z^ magnetic condensation for the untwisted SU{N) factor diagonally 
embedded, because, as we have already remarked, the untwisted sector has zero momentum on the 
torus and therefore it may condense in the vacuum of the localized theory. 

However, we will find in sect.(12) that the RG-flow in a given U{^) Morita equivalent com- 
mutative YM theory must change the rank, N, and the degeneracy, A'^, of the untwisted SU{N) x 1^ 
sector, keeping the product, A^A'^ = ^, constant, in order to define non-trivial finite correlation 
functions of composite surface operators. Yet, the Umiting bundle in the large ^ limit is, indepen- 
dently on «, a version of SU {oo^ embedded with oo multiplicity mU{ooy 

Because the Morita equivalent twisted connections hve in the tensor product, u{N) x u{N), 
there is a natural embedding in the Morita equivalent theory of SU (N) connections with Z/v holon- 
omy in Z/v x 1^. 

For computational technical reasons it may be convenient to get rid of the twist by means of 
the thermodynamic limit, to get back (commutative) R^. Going back to the original gauge (that 
would be singular on the commutative torus but not on 7?^) we obtain: 

[daAp{z,Z,U,u)-dpAa{z,Z,U,u)+i[Aaiz,Z,U,u),Ap{z,Z,U,u)]] 

= L^apiP)S^'Hz-Zp) + ^'„pl^,fj (7.15) 

p 

that is the standard defining equation of a lattice of surface operators in presence of a central 
magnetic field. The central term is now irrelevant in the large-A^ limit ^^ since without the non- 
commutativity it does split, at difference of Eq.(7.7). 

The solutions of Eq.(7.11) describe surface operators with singularities supported on the {u,u) 
plane and, because of translational invariance of the vacuum ^^ in the {u,u) plane, reduce to the 
standard two-dimensional Hitchin equations in the (z,z) plane. The large-A'^ and large-0 limit of 
the non-commutative theory is therefore recovered as a double large-A^ and large-^ limit. 

For further use we need to know that the second Chern class has an extension to surface 



operators as a parabolic Chern class. In the notation of |g^ : 

■ j d^xFapFap = Q + Y,trf{apmp) + -Y^DpHOptrfiaj) (7.16) 



16;r2 . 



where Q is the usual second Chern class of the U (N) bundle without the parabolic structure, Up is 
the vector of the parabolic weights at the point p, i.e. the vector of the eigenvalues of F^y divided 
by 271 modulo 1 in the fundamental representation, trip the magnetic flux through the surface Dp of 
the singular divisor p x Dp of the surface operator and Dp n Dp the index of self-intersection of the 



surface Dp. 



There is one more symplectic form that plays an important role in this paper. It occurs as the 
symplectic form associated to the twistor connection: 



(Op 



^ld\trftr{5Bp,A5Bp{) 

(Oi - ip {(Oj + icOK) -ip^\(Oj- icOK) (7.17) 



^^In sect.(12.7) it still plays a role. 
We mean in fact the vacuum on which twistor Wilson loops can be localized on. 



62 



Glueballs in large-N YM by localization on critical points Marco Bochicchio 



For p = — 1 it will be employed as an ingredient of the holomorphic/antiholomorphic fusion in 
sect.(12). Indeed it follows from Eq.(5.4) that co i = ft) depends holomorphically on /x_i. 

In sect.(12) we employ a modification of (O, co', defined over a punctured sphere rather than 
over its compactification obtained adding the singular divisor. The relation between the two forms 



is (Eq.(3.30) of [|120|]): 

(0 = co' + Y^Tr{^p{dgpgp')^) (7.18) 

p 

where the terms in the sum over p represent the Kirillov forms on the adjoints orbits at p. co' 



depends only on the holonomy of the connection (Eq.(3. 13) of [ ]120| ]). 

To summarize we have started from the point-like parabolic singularities of the non-commutative 
Eguchi-Kawai reduced theory and we have ended with surface-like singularities in the Morita 
equivalent commutative theory. We can turn the argument around and say that the aforemen- 
tioned point-like parabolic singularities of the non-commutative partial large-A^ Eguchi-Kawai re- 
duction are daughters of codimension-two singularities of the four-dimensional parent gauge the- 
ory. Codimension-two singularities of this kind have been introduced in [|^, |^] in the pure YM 
theory as an "elliptic fibration of parabolic bundles" for the purpose of getting control over the 



large-A^ limit of the pure YM theory exploiting the integrability of the Hitchin fibration. In [28] 
they have been introduced in the ^ = 4 SUSY YM theory for the study of the geometric Lang- 
lands correspondence, under the name of "surface operators", and this is now the name universally 
used in the physical literature. 

7.3 Moduli of surface operators 

To study the moduli space of surface operators in Eq.(7.11) it is convenient to compactify the 
(z,z) plane on a sphere. The moduli space has three different equivalent descriptions that are all 
employed in this paper. There is a vast mathematics [79, 8C, ^, |2, 83, 84,|85, |^, ^ and physics 



literature [p8, 91] on parabolic Hitchin bundles ^'^ . Thus we summarize briefly the essential results 



[79,80, 



The first description of the moduli space is of differential geometric nature as a Hitchin system 
and hyper-Kahler quotient, that in our approach follows by the the non-SUSY non-commutative 
Nicolai map on a dense set, as we just discussed. This is the description that occurs combining 
the quasi-localization lemma with the idea of integrating on surface operators. In the hyper-Kahler 
description the structure group of the bundles involved is compact. In our case U{N) or SU{N). 
Thus we refer to the gauge fixing in this framework as the unitary gauge. It is convenient to write 
Eq.(7.11) in non-covariant notation exactly as Hitchin equations (we disregard for the moment the 
central U{\) extension that splits): 

-iFA-[Au,Au] =£Mp"5P)(z-Zp) 
p 

- idj^Aa = ^ «p 5 (^' (z - Zp ) 
p 

-id-,A, =Y^np5^^\z-Zp) (7.19) 



47 



These references are by no means a complete list. 
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Because of the delta function at p in general the gauge connection has a pole singularity. The triple 
{lJ.p,np,np) determines the coefficients of the leading behavior of the gauge connection around the 
pole. The local model arises by restricting to such leading behavior [85, 86]. Since /i!) is Hermitian 



it is always diagonalizable. Let us consider first the semisimple case for which by definition all the 
eigenvalues of jip ^^ are different modulo n. 

A study of the local model implies that in this case also np and Mp can be diagonalized si- 
multaneously with jXp by a compact gauge transformation, gp [85]. This is a quite remarkable 



fact, referred to in this paper as local abelianization ^^ and it is the ultimate reason that allows the 
explicit computations of sect.(12) in the large-A^A^ limit. 

Thus the matrix, Up = jXp+Up — hp, commutes with its adjoint, jip, i.e. it is normal. Hence the 
compact adjoint orbits at a point, gp Apg^ ', where Xp are the (in general complex) eigenvalues oi IJLp, 
label some moduli of the solution of Eq.(7.15). There are other moduli that are not immediately 
manifest in the unitary gauge. They arise as moduli of the metric of the Hitchin bundle that is 



implicit in the definition of the unitary structure []86p. 

When some eigenvalues of ji^ are degenerate modulo n the asymptotic behavior of the con- 
nection changes. This is recalled below after describing the other representations of the moduli 
space. 

The second description of the moduli space is of holomorphic nature. It arises by a meromor- 
phic connection in a holomorphic gauge. Indeed the Hitchin equations imply the flatness equation: 

-iF{B) =£Mp5P)(z-Zp) 

F{B) =d,B,-d,B, + i[B„B,] (7.20) 

for the non-Hermitian connection: 

B^=A^ + iAu 

B-,=A-, + iAa (7.21) 

The moduli space arises as the Kahler quotient of the space of solutions of the flatness equation, 
Eq.(7.16), with respect to the action of the complexification of the gauge group. Because of a well 



known result [ ]52[ , |53P it coincides with the hyper-Kahler quotient of the three equations, Eq.(7. 15), 
with respect to the action of the compact gauge group [|^]. The structure of the moduli space is 
particularly transparent in a holomorphic gauge: 

S, = (7.22) 

In this gauge B^ is a meromorphic connection: 

id,B, = '£^'p5^^\z-Zp) (7.23) 

p 

with residue at p determined by jx' that is conjugate to jXp by a gauge transformation in the com- 
plexification of the gauge group. This description is the most transparent to understand the moduli 



^'^We have normalized flp in sucii a way tiiat the holonomy around p is e^'^^p . 
A proof in the physicists style of the result in | |85| ] about the commutativity of the triple (pLp,np,np) can be found 
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space because all the local moduli are labelled by the adjoint orbit in the complexification of the 



gauge group, /x^ = GpXpG 



-1 
p " 



The holomorphic description arises in the holomorphic loop equation. 

It implies also that Hitchin equations are associated to local systems, i.e. to fiber bundles with 



locally constant transition functions |79, 112]. A local system on a complex curve is the same as 
a representation of the fundamental group of a Riemann surface with punctures [^, |112| ]. This is 
the topological description of the moduU, and it is also the easiest to understand globally. Indeed 
the residues of the meromorphic connection, B^, determine its holonomy around p: 

Mp = Pe'^^i''"'''' 

= e^'^'p (7.24) 

The global moduli space on a punctured sphere is therefore the quotient of the algebraic variety: 

Y[Mp = 1 (7.25) 

p 

modulo the adjoint action of the complexification of the global gauge group (this description has 
been employed in sect.(l)). 



We come now to the non-semisimple case \\7% |85| , |86| ]. 



If the eigenvalues of the holonomy around p, e^'^i' , are not all different, the holonomy cannot 
be diagonalized in general but it can be put in Jordan form. In this case in the unitary gauge some 
eigenvalues of np and of its Hermitian conjugate Up are degenerate as well and the Higgs field, A„, 
in some directions in color space has not anymore a pole singularity but only a milder one, a pole 
divided by powers of a logarithm. The power of the logarithm and the coefficient of the pole are 
determined by the off-diagonal parameters in the Jordan form of the holonomy JS^]. 

Perhaps the most important property of the hyper- Kahler construction from a physical point of 
view in the semisimple case is the fact that in a unitary gauge the coefficients of the delta function 
at a point {iJ.p,np,np) commute and thus can be diagonalized at the same time by a unitary gauge 
transformation. 

On the opposite in a holomorphic gauge the residues of the meromorphic connection and 
the local holonomies of the twistor connection cannot be diagonalized in general by a unitary 
transformation but generically only by a transformation in the complexification of the gauge group. 
Thus there is mismatch in the number of local degrees of freedom between the holomorphic and 
unitary descriptions. 

This is explained by the fact that in the unitary description the missing degrees of freedom 
arise as moduli of the Hermitian metric associated to the Higgs field [^6|]. Physically this means 
that the dimension of local fluctuations of the ASD curvature, jXp, of semisimple type that occur by 
the hyper-Kahler construction in a unitary gauge is just one-half of the dimension that occurs in 
the holomorphic description of local fluctuations, jj.' This leads to a non-trivial Jacobian from the 
unitary to the holomorphic gauge whose logarithm turns out to be the glueball potential. 

In the next section we consider solutions of the Hitchin equations for connections with Z^ 
holonomy. These connections have no (local) moduli since the adjoint orbit of the center is the 
center. They turn out to be the Hitchin bundles that occur at the fixed points of the quasi-localization 
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lemma. They satisfy the Hitchin equations: 

-iFA-[A,,Aa] =£Ap5(2)(z-Zp) 
p 
BaAu =0 

h^u =0 (7.26) 

with e^^^i' G Zn. Therefore: 

iXp = diag{2n{k-N)/N,2nk/N) {1.11) 



N-k 



These equations are invariant for the following L'^(l) action: 



Au ^e'%, 



Aa ^e-'^Aa {12%) 

Since there are no moduli this L'^(l) must act by gauge transformations: 

geAug^^ =e'^Au 
geAug'^^ =e''^Au 

g^Kgt =^z 

geA^ggi =A, (7.29) 

More generally the set of moduli fixed by this [/(I) action is the Lagrangian submanifold of the 



hyper-Kahler moduli space for which the Higgs field, A„, is nilpotent [113 



There are fundamentally two interesting types of orbits in the Lagrangian cone of the hyper- 



Kahler moduli space [ ]113| ]: the orbits with unitary holonomies, for which the Higgs field vanishes 
identically; the orbits that correspond to Hodge bundles, for which the holonomies are valued in a 
real version of the complexification of the gauge group [ ]113| ]. 

In the first case the holonomies can always be diagonalized, despite the eigenvalues may not 
be all different. In the second case the holonomies cannot be diagonalized, but can be set in Jordan 
form. Both the orbits play a role in the computation of the Wilsonian beta function (sect.(9)). 



8. Localization on fixed points in large-A^^ YM 

We use the description of surface operators as local systems (i.e. as representations of the 
fundamental group) to obtain localization on fixed points. 

Indeed for the lattice theory of sect. (7) we can justify the exchange of the order of integration 
and limits that occurs in the quasi-localization lemma of sect.(5), taking advantage of the existence 
of the singular gauge to reduce the lattice theory to a locally abelian theory, whose eigenvalues 
label the gauge orbits of the hyper-Kahler reduction and do not fluctuate in the large-A^A^ Umit. 
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In the aforementioned singular gauge the partition function reduces to an integral over the 
eigenvalues and the zero modes of the locally abelian theory: 

|2 
"PI 



2 = I j e-^]\8Xp8Vp8Vp 

I' - - 4NN — , 

= I / 5'Ad'Ad'D5'Dexp{-^Y.^rNTr^{XpXp + 4VpVp)) 



ANN . 
Sw p 
5{-iFB-UK-Vp + Vp)5^'\z-Zp)) 



d{-idAD-l^Vp5^^\z-Zp))d{-i~dAD-l^Vpd^^Hz-Zp)) 



p p 

M^p + Vp) .«.„*1-rs:. s... s..-. |2 



H^p -Vp + Vi 



A"''co-Y[S?ip5Vp5Vp\ (8.1) 



'P "p ^ "py p 

where {Xp,Vp) are the eigenvalues of {iXp,np) and A(A) is the Vandermonde determinant of the 
eigenvalues, Eq.(12.23). The Vandermonde determinant that occurs in the numerator is due to 
gauge fixing np in triangular form, by the action of the unitary gauge group in the singular gauge. 
As a consequence np is automatically diagonal since it arises from the solution of the Hitchin 
equations, sect.(7) and sect.(12). The Vandermonde determinant in the denominator arises because 
of the combination of gauge fixing in the singular gauge and of the change of variables to the 
holomorphic gauge, sect.(12). 

At large-A^A'^ it is not restrictive to assume that only one set of eigenvalues (Xp, Vp) or a discrete 
sum of them actually contribute to the partition function: 

z = 1.-^1^ 



r - - 4NN ~ ~ 

I / 5'A5'A5'D5'Dexp{-^Y.frMTr^iXpXp + 4Vp^p)) 



ANN . 
Sw p 
8{-iFB-Y^{lp-Vp + ~Vp)8^^\z-Zp)) 



5{-idAD-Y^Vp5^'-\z-Zp))8{-i~dAD-Y^X'p5^'-\z-Zp)) 



p 



Mvp + Vp) .^„,^-|2 ^3_2) 



A(Ap -Vp + v 



p> 



Therefore we can now justify the interchange of order of limits and integration in the quasi- 
localization lemma (sect.(5)), since in fact at large-A^A'^ there is no integration variable in the ar- 
guments of the delta functions that depend on the ASD lattice variables: 



5'A5'A5'D'5'D'txp ( - ^ £?rA,rr^(Apip + 4Vp^p)) 

Sw p 



1/ 

TryPexp / / (A, + D',)dz + (A,- + D'a)dz 

5{-iFA + [D\D']-i-dAD' + i^5AD'-Y,{lpd^^Hz-Zp)-p-'vpd(^\z-Zp) + p^pd^^Hz-Zp))) 
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8{-adAD' -Y^Vp8^^\z-Zp))5{-a-'~dAD' -Y^VpS'^'-Xz-Zp)) 



p p 



A{Xp-p-^Vp + pVpj 
Thus we get at the fixed points: 



Z =\e-^\' 



I / 5A'5A'5D'5D'exp{-—j-Y^trNTr^iXpXp)) 

■J Rw n 



ANN, 
5{-iFB-Y^lp5^^\z-Zp))5{dAD)5CdAD) 



^A-«*|^ (8.4) 

It remains to integrate on the moduli and determining the eigenvalues, Xp, that is resolved requiring 
that the global gauge group be unbroken at the critical points. Indeed physically the basic idea is 
that the twistor Wilson loops in the adjoint representation, of both the non-commutative theory with 
gauge group U{N) and of its Morita equivalent counterpart, satisfy the general criteria of 't Hooft 
duality. Thus if the gauge group is unbroken the twistor Wilson loops in the adjoint representation 
must be localized, by large-N factorization, on a condensate of the tensor product of Z^ magnetic 
vortices with the conjugate representation. We have noticed in sect.(7) that the correct ansatz for 
the localized locus is in fact Z^ x 1^, i.e. Z^ occurs with degeneracy N. 

The formal argument is as follows. The complexification of the global gauge group acts on 
the holonomy at one point, pi,by the adjoint action, in such a way that Mp^ can be put in canonical 
form. Mp^ can be diagonalized if it has distinct eigenvalues, while in general it can be put in Jordan 
form. 

In the laige-N limit it is possible to restrict the integration measure, djU-p, to orbits whose 
holonomies have fixed eigenvalues, since this restriction implies an error of subleading order in 
jj. In addition by translational invariance the conjugacy class of the orbits at all the points p must 
be the same. Finally the global SU{N) gauge group (and not only the L'^(l)^^^ torus as in the 
Nekrasov case) must fix Mp^ , i.e. gMp^g^^ = Mp^ , since otherwise Mp^ would break spontaneously 
the global gauge symmetry. Therefore Mp^ must be central and thus must be in Z^. But then all the 
orbits collapse to a point and there are no moduli at the fixed points ^^. 

However, in any neighborhood of the fixed points of the global gauge group, the orbits are 
non-trivial and moduli there exist. There are essentially two different ways to describe these neigh- 
borhoods of the fixed points. One possibility is deforming infinitesimally the unitary part of the 
eigenvalues of the holonomy, the other possibility is deforming the holonomy along nilpotent di- 
rections. Both possibilities are discussed in the next section. 

Thus if we first compute the effective measure in a neighborhood of the fixed points and then 
we sit on the fixed points the induced measure will contain the powers of the Pauli-Villars regulator 
due to the moduli. This has an analog in the localization of the ^ = 2 SUSY YM partition function. 



As we have seen in the introduction the requirement of translational invariance can be relaxed. 
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where generically instantons have moduU (this is essential to get the correct beta function in that 
case too), but the instantons at the fixed points of the torus action have not. 

Thus at the fixed points the contour integral over /Xo+ in Eq.(5.8) collapses to a discrete sum 
over sectors with Z/v holonomy. The reduced Eguchi-Kawai effective action of the localized theory 
is now: 



Der^-(-AAdaB-iad,,-)Det(-AA)(:^)"''Der-(oJ^xc.c.] _„ (8.5) 

The connection, A, denotes the solution of the equation [F „ —T,plJ-^g{p)5^{z — Zp,-.) = 0]„=,i=o 
in each Z/v sector. Det^co is the contribution of the ni, zero modes due to the moduli and A the 
corresponding Pauli-Villars regulator. The complex conjugate factor arises by the conjugate repre- 
sentation. 

Thus the holonomy of adjoint twistor Wilson loops at the fixed points is trivial, because of the 
cancellation of Zp^ factors between the fundamental and the conjugate representations. Hence for 
twistor Wilson loops the same result is obtained performing the limit A — ;• outside the functional 
integral, leading to triviality via A -independence and the all order argument of sect.(4.5), and inside 
the functional integral, leading to localization on the tensor product of Z/v surface operators and the 
conjugate representation, and to triviahty as well. 

9. Wilsonian beta function 

9.1 Beta function by restricting the non-SUSY Nicolai map to surface operators 

We now proceed to the computation of the beta function. In order to define the renormalization 
of the coupling constant it is necessary to compute the classical YM action of surface operators. 
As we have seen in sect.(7) surface operators are defined by connections on parabolic bundles. In 
a mathematical sense we can think of parabolic bundles in two different ways. Either parabolic 
bundles occur on space-time with no boundary and with a divisor and a parabolic structure that 
belong to the space-time. This is the point of view in this paper and in some mathematical literature. 

Or they arise on space-time with boundary, where the boundary is the parabolic divisor. This is 



the point of view of [28]. In the latter case the insertion of a surface operator keeps the finiteness of 
the action, since the singular parabolic locus is not included in the space-time integral that computes 
the action. This justifies also the term operators, since their occurrence is the analog of operator 
insertions ^^ 

However, our point of view is that the surface operators are dynamical objects and therefore 
their singular divisor is included in the path integral in space-time. 

Despite the parabolic singularity, the topological term in the action has a well defined mathe- 
matical extension to parabolic bundles, as parabolic Chern class (Eq.(7.16)). 

This is not the case for the term involving the ASD field. As a consequence the classical YM 
action is quadratically divergent on each singular divisor, p x Dp, of a surface operator, with a 



We would like to thank Edward Witten for a discussion about this point. 
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divergence proportional to the area of the singular locus of each surface operator. Therefore we 
need a way to handle this classical divergence. 

We have already recalled in sect.(4.1) that when the codimension-two surface is non-commutative, 
as in our case, the YM action of the corresponding non-commutative reduced Eguchi-Kawai (EK) 
model is rescaled by a power of the inverse cutoff, that cancels precisely ^ the quadratic diver- 
gence that occurs evaluating the classical YM action on surface operators. This allows us to define 
a new kind of semi-classical computation for which the classical YM action is finite on parabolic 
bundles. In our case the EK reduction is only partial in such a way that the action is: 

^N{^)^Jd^xtrNTrf,{-i[da + iAa,dp + iAp] + d^'p\f (9.1) 

where the trace Tr^ is taken over a subspace of dimension A'^, with 

N(—f = 2Kd (9.2) 

A 

in the large N, 0, A limit. The contribution of each parabolic singularity in the action reads: 

dh5^^\x-xpf 



5(2) (0) 1^2^5(2) (x-X^ 



^(A)^ (9.3) 

and it is cancelled by its inverse in the EK reduced action. It is convenient to describe the same 
result in terms of the action of a commutative gauge theory. In this case the action reads: 

^N J d^xtrNTr^{-i[da + iAa,dp + iAp] + 0„^ 1)^ (9.4) 

and it is divergent on a surface operator as: 

d^xS'^^^x-Xpf 
= 5^^\0)V2 fd^xd^^\x-xp) 

= A^2 (9.5) 

where V2 is the area of the singular divisor of the surface operator and the last equality is the 
definition of N2- Thus the reduced action is related to the one of a commutative gauge theory by 
the factor of A/^2^ ^2. 

Y^Ni^f I d^xtrNTri^{-i[da + iAa,dp + iAp] + 0-^^ if 
= :^^N I d^xtrMTrj^i-i[da + iAa,dp + iAp] + d^'^lf (9.6) 



52 



Precisely the same factor arises in the quenched version of the EK reduction [ |l06[ ]. 
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We can interpret this result by saying that the normalization of the trace in the non-commutative 
EK reduced theory differs by a factor ofN2^ from the one of a commutative theory. Thus it may be 
computationally convenient to perform the calculations in the commutative theory and at the end 
normalize appropriately the traces. 

Let us evaluate firstly the action of the surface operators of Z^ holonomy that occur at the 
fixed points. The Morita equivalent commutative ^^ theory has gauge group U{N x N) and the 
center of SU{N) is embedded diagonally in L'^(A'^ x N) as e~i^ 1^. For a surface operator of Z^ 
holonomy around the point p of charge k, i.e. such that Mp = e'~N~ ,N — k eigenvalues, iXp, of the 
ASD curvature at p, Fqj = 2Xpd^^' (z — Zp), are equal to ^ and k eigenvalues are equal to ^^~ ' , 
for the curvature to be traceless and to give rise to the holonomy Mp = e~N~ ^'*. The trace of the 
squaie of the eigenvalues of the ASD curvature in the fundamental representation is thus: 

,.,(4A^, = (^_,)(^). + ,(?M)^ = p.)^^ (9.7, 

in such a way that the reduced action is: 

and the Morita equivalent one is: 

NN{Anf N^ 



nm 



Q + Tr^N2Y,2{lKfk{N-k) (9.9) 



2g2 - ■ 2g ^ 



From these equations it follows that, if Q is finite, its contribution to the action is irrelevant with 
respect to the one of the ASD curvature at the parabolic singularities. 

Once the classical quadratic divergence has been tamed by the EK reduction we need to un- 
derstand the logarithmic divergences that lead to a non-trivial beta function in the Jacobian for the 
change of variables from the connection to the ASD curvature. These divergences have been al- 
ready computed in sect.(3.4), and we can just adapt our previous calculation to the case of a lattice 
of surface operators. We have already observed that, as in Nekrasov localization, we should evalu- 
ate and renormalize the functional measure in a neighborhood of the fixed points and thereafter we 
should sit on the fixed points. Since the fixed points have no moduli and thus no zero modes, the 
inverse order, first sitting on the fixed points and then renormalizing the functional measure, would 
not lead to the correct result. 

But let start ignoring for the moment the zero modes and sitting on the fixed points. The 
observation that the contribution of Q to the classical and the quantum effective action is irrelevant 
with respect to the one of the ASD curvature at the parabolic singularities resolves the issue about 
the different factors of Z^' and Z^^ that occur in sect.(3.4) as counterterms in the Jacobian of the 
non-SUSY Nicolai map. 



^■'We rescale the area of the Morita equivalent theory from 4j to L? in order to perform the thermodynamic limit 
uniformly in N. 

The curvature in not uniquely determined by the holonomy, since parabolic bundles admit extensions over the 
punctures such that the eigenvalues of the ASD curvature differ by shifts of In. Our choice is in some sense minimal. 
This feature together with many others, known in the mathematical literature, is reviewed in [Ea]. 
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This observation allows us to ignore also issues related to the possible global non-triviality 
of the bundles and to the counting of the global moduli, and to concentrate only on the local 
counterterms associated to the sum over the points of the parabolic divisor. 

As in every computation involving an effective action the background field should live on 
a scale much larger that the quantum fluctuating field. This may look awkward to realize for 
surface operators that carry delta-like singularities, that involve any momentum scale. In fact it is 
impossible to realize for an isolated surface operator, but it is possible for a lattice with uniform 
lattice spacing ^^. A typical example is the following one-loop contribution to Z ^ in Euclidean 
configuration space, that arises from the spin term of sect.(3.4) evaluated on surface operators: 

,, T,^ > d ud v-. ,T , rr-^ (9.10) 

where the sum over p,p' runs over the planar lattice of the parabolic divisors of the surface opera- 
tors. There is also an orbital contribution that has the same structure. Indeed the orbital logarithmic 
contribution to the beta function arises from terms of the kind ^^: 



-yr ''"{^-y) 



4 



|2 



2 2 2 2 4 Z ~ w\ 

d zd wd ud wTr{AJx)Ay{y)- -7: — ; pr-r" 

{\z-w\^ + \u-v\^Y' 

d'^zd^wd'^ud'^wTr{A^{x)Ai{y)did, 



6(|z-w|2 + |m-v|2)2' 



4 
d^zd^wd^ud^wTr{d^d„{A^{x)A-^{y)) 



-(|z-w|2 + |m-v|2)2^ 

d^zd^wd^ud^wTr(dyAJx)dwAy(y)-—, p: — | pr-^) 

6(|z — wp-l- \u — v\^)^ 

~ £ / d\d\- 1. V |,.2 (9-11) 

^p>J {\zp-Zp'\^ + \u-v\^Y 

where in the last line we used 5|Aj.(x) ^ £p Mp5(2) {z — Zp) and d„Ai{y) ~ — l^p Mp5(2) {w — Zp) for 
the surface operators that occur at the fixed points, withx = (z,z, m,m) andy = (w,w,v,v). 

We would like to find a regularization for which the loop expansion of the functional determi- 
nants evaluated on surface operators satisfies the usual power counting as in the background-field 
computation of the beta function. 

To avoid quadratic divergences we restrict the sum to p ^ p' . Although quadratic divergences 
in the Morita equivalent theory correspond to finite counterterms in the EK reduced action, we 



^^In fact it is sufficient that the lattice spacing be uniform on a scale much larger than the typical scale of the physics 
involved. 

^^The term involving daAa vanishes identically around the local singularity, while the term involving A^ is quadrat- 
ically divergent and does not contribute because of cancellations due to gauge invariance in any gauge invariant regu- 
larization of the theory. Since in our computations only functional determinants occur, the Pauli-Villars regularization 
suffices. 
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would like to avoid divergences at coinciding points in higher orders of the loop expansion for the 
conventional power counting to hold. We set further /ip = Ap = A in such a way that : 



is logarithmically divergent (both in the ultraviolet and the infrared). Indeed introducing a lattice 
scale a: 



1 V /■^2 ,2 NTrjX^) 

i^^'Y^J {\zp-zp>V + \u-v\^Y 

1 f ^2 -2 a -2 f ^2 j2 NTr{X' 

a ua a za a wcl v- 



{Ak^YJ ^ y (|z-w|2 + |m-v|2)2 

NTr{X'^) 



-^nI j cYwcYv- 



{4n^Y"^J ' (|z-w|2 + |m-v|2)2 

' N^NT r{X^) log - 

■■N2Y,NTr{X^)log- (9.13) 

p M 



One factor of A'^2 is just the sum on lattice points, the other factor is the "phase space area" A^2 

^271' 



(^)2V2 of one surface operator. From the preceding equation we read also that: 



a-' = A (9.14) 

271 

To summarize, there exist a point-splitting regularization of the effective action in the background 
of the lattice of surface operators, that together with the implicit use of the Pauli-Villars regu- 
larization, needed to handle quadratic tadpoles, leads to the same power counting as the usual 
dimensional regularization, with the logarithmic divergences occurring as in Eq.(9.12). Had the 
contributions with p = p' been included, there would appear quadratic divergences, thus spoiling 
the usual power counting in higher order terms of the loop expansion. This lattice point-spitting 
regularization ^^, followed by Epstein-Glaser renormalization in Euclidean configuration space (see 



[ ]114| ] for references) is a possible starting point of a new constructive approach for the laige-N YM 
theory. 

We should understand now the contribution of the zero modes. 

At first we count the moduli of surface operators and then we give an argument to identify the 
zero modes with the moduli, as for the instantons. 

From the computation of the Z^ factor that we already performed in sect.(3.4) it follows that 
the renormalization of the action in absence of zero modes would be: 

Adding to the action the contribution of the complex conjugate representation we get: 

2gw(M) 2g2,(A) (471)2 3 \i 



This regularization has been found during joint work with Arthur Jaffa. 
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9.2 Dimension of the Lagrangian neighborhood of the fixed points 

It follows from Eq.(9.16) that, in order to get the correct one-loop beta function, the contribu- 
tion of the zero modes to the renormalization of the action should be —2k(N — k) log(^), including 
the fundamental and conjugate representation. The sign is consistent with the Pauli-Villars reg- 
ularization of zero modes, yet the absolute value of the coefficient of the logarithm is in general 
an even integer but not a multiple of 4, as it would be implied by the hyper- Kahler reduction ^^. 
Thus the neighborhood of the fixed points cannot be generic. Remarkably we have already seen in 
sect.(7) that the fixed points sit automatically inside the Lagrangian cone of the moduli space for 
which Au is nilpotent. Since in a Lagrangian neighborhood of the fixed points the dimension of the 
moduli space is generically one half of the dimension of a hyper- Kahler neighborhood, the correct 
beta function may arise. 

We should classify now which are the Lagrangian neighborhoods of the fixed points that lead 
to the correct beta function. 

One component of the Lagrangian cone corresponds to A„ = and gives rise to unitary rep- 
resentations of the fundamental group. The complex dimension of an adjoint orbit for a generic 
parabolic unitary bundle of rank A^ is given by: 

&imx = \{N^-Y,mj) (9.17) 

where m, are the multiplicities of the eigenvalues. Vortices of Z^v holonomy have no moduli, since 
the holonomy lives in the center, for which the multiplicity of the eigenvalues (modulo 2n) equals 
the rank. But the following slight deformation of the eigenvalues gives rise to a non-trivial adjoint 
orbit for the holonomy: 

2X=diag{2n{k-N)/N + e,2nk/N-ek/{N-k)) (9.18) 

^ V " V ' 

k N-k 

Thus the complex dimension of the orbit is: 

dim;L = -{N^-k^-{N- kf) = k{N - k) (9.19) 

and the same holds for the complex conjugate orbit, in such a way that the real dimension of the 
orbit matches the number of zero modes needed for the correct beta function: 

16.^M^-t) ^, 1 1 5 A ^^^^ 

However, there is a more intrinsic characterization of the Lagrangian neighborhood of the fixed 
points. Instead of deforming slightly the eigenvalues we may deform the moduli along nilpotent 
directions. Hence we may require that, in a unitary gauge in the Lagrangian neighborhood, exactly 
the same Hitchin equations are satisfied as at the fixed points. Therefore the eigenvalues of the ASD 
field are precisely: 

2X=diag{2n{k-N)/N,2nk/N) (9.21) 

k 



A hyper-Kahler manifold has necessarily a real dimension that is a multiple of 4. 
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but the we allow the Higgs field, A„, to have a nilpotent residue. In the Lagrangian cone not only the 



Higgs field has a nilpotent residue, but it is nilpotent itself ||8g]. These are bundles of Hodge type 



[ ]113| ] for which the twistor connection has a holonomy that cannot be diagonalized, but it can be 
set in Jordan canonical form. In the Lagrangian cone the local moduli at a point for these bundles 
are parametrized by orbits of the Jordan canonical form for a real version of the complexification 
of the compact gauge group [ |113[ ]. Thus in our case, for which the diagonal part of the holonomy 
is in Zn, the local holonomy is unipotent. 

Since the diagonal part of the holonomy is central, to compute the dimension of the orbit we 
need to consider only the nilpotent part. Any such matrix is conjugate by the Jordan theorem to a 
direct sum of k blocks of dimension dt, such that jj^^i di = N, where N is the total rank. Each block 
has di zero eigenvalues on the diagonal ^^ and it is upper triangular with all 1 on the super-diagonal. 

The classical action is not modified at all by a deformation of the moduli along a nilpotent 
direction, since the nilpotent part of the holonomy is invisible in the ASD curvature in a unitary 
gauge [||]. 

Therefore, to get the correct beta function, we need to construct orbits of nilpotent Jordan 
matrices, such that the real dimension of the orbits for the action of a real version of the complexi- 
fication of the gauge group is precisely the double of the complex dimension of the unitary orbits, 
i.e. of the flag manifolds. 

Indeed in this case the conjugate representation describes exactly the same moduli, since the 
orbit is for the action of a real version of the gauge group. 

Nilpotent orbits with the features just described, having double the dimension of a flag, are 



called Richardson orbits [115] 



Here are some examples. The principal nilpotent complex orbit, i.e. the orbit of a nilpotent 
Jordan block of maximal rank, has precisely the same dimension as the complex orbit in the generic 
semisimple case, i.e. N'^ —N = N{N— 1), that is always even. Thus the real dimension is a multiple 
of 4, as it should be. 

By Eq.(9.17) the complex dimension is the double of the complex dimension of the maximal 
flag, and thus the principal nilpotent orbit is a Richardson orbit. 

We are looking for a nilpotent orbit with double the dimension of the partial flag in Eq.(9.19). 



Let us first recall the general formula for the dimension of the adjoint orbit of a nilpotent [116]: 



k 
dimON=N^-N-lY^{i-\)di (9.22) 

i=\ 



For example for the principal nilpotent, k = \, d\=N, and we get the aforementioned result [ |117| ]. 
For the zero nilpotent, k = N,di = I, and we get 0. The orbit of the direct sum of k nilpotent Jordan 
blocks of dimension 2 and of one Jordan block of dimension N — 2k has double the dimension of 
the partial flag in Eq.(9.19). 
Indeed in this case: 

k N-k 

dimON=N^-N-2j^{i- 1)2-2 £ (/-I) 

(=1 i=k+\ 



A nilpotent matrix has all the eigenvalues zero. 

75 



Glueballs in large-N YM by localization on critical points Marco Bochicchio 



*: N-k 

■.N^-N + 2{2k + N-2k)-2Y^i-2 ^ i 

i=i i=k+i 

■ N^ - N + 2N - {N - k+ \){N - k) - {k+ \)k 



■ N^ + N - {N - k+ \){N - k) - {k+ \)k 
■-2k{N-k) 



(9.23) 



QED 

We should clarify in which sense the moduli that label the orbits with fixed eigenvalues occur 
as zero modes. Let us consider first the unitary orbits in the Lagrangian cone. In this case in any 
smooth unitary gauge the ASD curvature, /Iqj (p), is conjugated to the eigenvalues by the action of 
the unitary group. 

However, the classical action depends on the eigenvalues but not on the unitary matrices. Thus 
the unitary matrices define flat directions in the classical action, but the holonomy of the twistor 
connection actually depends on the unitary matrices that label the orbit. 

We can associate to these moduli zero modes of the functional determinants choosing the 
singular gauge in which the gauge curvature is diagonal. This gauge may be singular in the sense 
that may be reached by gauge transformations that are possibly singular along lines, i.e. semi- 
infinite strings that start at the punctures and end at infinity or with another puncture. From the point 
of view of the homological localization of the holomorphic loop equation of sect. (10) these gauge 
transformations are allowed, provided the associated strings do not intersect the backtracking arcs 
that connect the loop to a puncture, in such a way that the twistor connection is not discontinuous 
across the backtracking arcs. 

In this unitary singular gauge the equations for the surface operators of unitary holonomy are: 

p 

with Aqj(/7) diagonal matrices and the other ^^gip) vanishing [^]. But the connection still de- 
pends implicitly on the moduli, in such a way that zero modes occur by the standard argument at 
the end of sect. (3. 3). 

In the case of Richardson orbits for Hodge bundles the holonomies are conjugated to the Jordan 
form by a real version of the complexification of the gauge group, that factorizes into a compact and 
a parabolic subgroup, the compact factor being a subgroup of the unitary group. In a unitary gauge 
only the orbit of the compact subgroup occurs as an adjoint orbit of the ASD field at a point, while 
the parabolic group parametrizes the moduli of the metric. Thus it there exists a possibly singular 
unitary gauge in which all the moduli occur as zero modes as for orbits of unitary holonomy. 

10. Homological localization of the holomorphic loop equation 

It has been known for many years that (twisted)-5'C/5'7 observables can be locahzed in gauge 
theories with extended SUSY by deformations that are trivial in the cohomology generated by the 
twisted super-charge [41, 42, 44]. Since cohomology is dual to homology [ ]118t ] ^^, we may wonder 



The main difference between cohomology and homology is the lack of a ring structure in the latter. 
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as to whether we can compute functional integrals by deformations that are trivial in homology 
rather than in cohomology. Were the answer be affirmative, we could get localization without 
supersymmetry. 

While there is no positive answer in local field theory, gauge theories contain many non- 
local observables, the Wilson loops. Thus the natural arena for homological localization in gauge 
theories, as opposed to cohomological localization, is the loop equation ^^■ 

In general the loop equation is the sum of a classical equation of motion and of a quantum 
term, that involves the contour integral along the loop. By homological localization of the loop 
equation we mean a deformation of the loop that is trivial in homology and for which the quantum 
term vanishes, in such a way that the loop equation is reduced to a critical equation for an effective 
action ^. Hence the needed homological deformation has to satisfy the following properties. 

It has to be trivial in homology. 

It has to leave the expectation value of the loop invariant. 

It has to imply the vanishing of the quantum term in the loop equation, i.e. of the term that 
contains the contour integral along the loop. 

In this homology framework there is a very natural analog of the operation of adding a cobound- 
ary in cohomology, that is based on the zig-zag symmetry of Wilson loops. The zig-zag symmetry 
is the invariance of a Wilson loop by the addition of a backtracking arc ending with a cusp. This 
deformation is a "vanishing boundary" in singular homology. In a regularized version the arc is 
the boundary of a tiny strip. While this symmetry holds classically in most of the cases, quantum 
mechanically the renormalization process may spoil it. The reason is that in general Wilson loops 
have perimeter and cuspidal divergences. The perimeter divergence is linear in the cutoff scale. 
The cuspidal divergence is logarithmic, with a coefficient that in turn is divergent for backtracking 
cusps. In SUSY gauge theories with extended SUSY there are examples of \oca\\y-BPS Wilson 



loops that have no perimeter divergence [|13[]. 

We have seen in sect.(4) that twistor Wilson loops share with their supersymmetric cousins 
these non-renormalization properties and, being trivial in the large-A^ limit, they have not cuspidal 
divergences either. Localization by homology leads this kind of non-renormalization properties to 
their extreme consequences. 

One of the virtues of the lattice regularization of the Nicolai map of sect. (7), from the point 
of view of the homological localization, is to allow identifying the cusps of the aforementioned 
backtracking arcs with the parabolic singularities of the reduced EK theory. We have noticed in 
sect.(7) that this point-like parabolic singularities are daughters of codimension-two singularities 
of surface operators of the parent four-dimensional theory. The lattice version of the holomorphic 
loop equation of sect.(6) follows ^' : 



< 



^'<«S;:7'^'<^-'» >= I L ^ " "''"'<"-■' "" "'"<'""' " "°-" 



where *P' is the holonomy of B in the gauge B'^ = 0. The lattice points associated to the divisor of 
surface operators become the cusps that are the end points, p, of the backtracking strings, bp, that 



^'The holomorphic loop equation is written in linear form since it is assumed that the loop Czz is simple, i.e. it has 
no self-intersections. 
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perform the deformation of the loop, C. Adding the backtracking strings implies the homological 
localization of the holomorphic loop equation: 






< Tr{ ^ .,;. ,\ *y^(Z.U [bp])) >= (10.2) 



The homological localization can be understood in the following geometrical terms. The existence 
of the regularized residue in Eq.(6.8) is the geometric obstruction to the localization of the loop 
equation. The regularized residue is the line integral of a current supported on a point and computes 
the de Rham cohomology of compact support in one dimension, H}{R) ^^, that coincides with the 
cohomology of a point, H^{pt) [ |118[ ]. It is interesting to display how this cohomology of a point 



can be obtained by a Mayer- Vietoris argument [jllSj] involving partitions of unity in the right hand 



side of the loop equation. For a smooth point of the arc (i.e. a point with a continuos tangent), 
introducing partitions of unity on the arc, \ = u\ {s) + U2{s), we get: 



dw^{s)5{z+{s,,n)-W+{s)) = U,{S,^) , 7:T(, +U2{s„„) , 7r„!."(, = 1 (10.3) 



W+{st,n) , .. / nJM^™J 
\w+{sJm)\ \w+{Ssm)\ 

Let us suppose now that we try to compute the "cohomology of a backtracking cusp". The same 
Mayer- Vietoris argument shows that such a cohomology does not exist in a classical sense, since 
the result depends on the choice of the partition of unity: 



dw+{s)8{z+{Scusp)-W+{s)) = Ui{Scusp)t—, TT + "2 (■ScM.sp 



W+{Scusp)\ ' \w+{Scusp)\ 

= U 1 {Scusp ) - "2 {Scusp ) ( 1 0.4) 

This is due to the fact that we can integrate distributions on smooth manifolds, but their extension 
to non-smooth ones depends on arbitrary choices in general. In particular, if the partition of unity is 
symmetric, u\{scusp) = U2{scusp) = \, the regularized residue vanishes. Therefore a regularization 
exists, that preserves the zig-zag symmetry of twistor Wilson loops, for which the holomorphic 
loop equation localizes. 

Thus if every marked point of the loop equation can be transformed into a backtracking cusp 
we can complete our argument about localization. But this is precisely the effect of our lattice, 
since marked points of the loop contribute to the loop equation in the lattice theory only if they 
coincide with the lattice points. 

We may think that it is a change of the conformal structure around the lattice points that 
generates the cusps. 

Since the Lagrangian submanifold on which twistor Wilson loops are supported on satisfies: 

\dzdz\ = \dudu\ (10.5) 

this two-dimensional conformal transformation lifts to a conformal rescaling of the four-dimensional 
metric: 

9 — — 

ds = dzdz + dudu (10.6) 



^ The integral of the delta function in one dimension can be approximated by the normalized integral of one-forms 
of small compact support around a point. 
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Thus it acts by adding a conformal anomaly to the effective action: 

r{[bp]) = r{[p]) +ConformalAnomaly{[bp\) (10.7) 

that amounts to a local counterterm, i.e. to a change of the subtraction point. 

Therefore there is a symmetry of the RG flow that generates the homological deformation of 
the loop by a vanishing boundary, i.e. by backtracking strings. This is the analog of the action 
being a closed form in cohomology, since in the last case there is a symmetry of the action (i.e. the 
twisted supersymmetry) that generates the coboundary. 

At this point we would like to clarify why we cannot use the Makeenko-Migdal {MM) loop 
equation to get localization. 

We can write the MM loop equation for unitary Wilson loops in the large-A^ YM theory as: 

= i I dxa I dya8^^\x-y) < Tr^<{x,y;A) >< Tr^<{y,x;A) > (10.8) 

where 

^'{x,y;A)=Pexpi f A^dxa (10.9) 

In the case of loops without self-intersections but with cusps the MM loop equation reduces to: 

= / / dxa j dyaS^'^^x-y) < Tr^'{x,x;A) >< Tr\ > (10.10) 



Performing the two contour integrations along the loop in the right hand side, we get [ ]13| ] : 



/ dxa<^,Tr{-^^^>{x,x;A))> 

~/(/'A^+ £ ^^!^(7r-a,„,p)A2) < Tr^'{x,x;A) > (10.11) 

cusp Siri^'ciup 

where P is the perimeter of the loop and D.cusp the cusp angle at a cusp. For our conventions 
^cusp = ^ for no cusp, while ^cusp = for a backtracking cusp. The perimeter divergence arises 
by the double integration of the four-dimensional delta function, i.e. of the contact term, along 
the loop. Integrating the contact term in a neighborhood of each cusp gives rise to a sub-leading 
quadratic divergence, since around a cusp there are two independent integrations instead of one, 
due to the two sides of the cusp. The coefficient of the cusp contribution is proportional to the ratio: 

cos^ (10.12) 

The numerator arises from the scalar product, the denominator from the two independent integra- 
tions of the two-dimensional delta function. In the Umit in which the cusp angle ^cusp reaches 
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the cusp backtracks and the cusp contribution to the contact term of the MM loop equation is di- 
vergent. It turns out to be proportional to the coefficient of the logarithm in the cusp anomaly [ll3|]. 
Therefore for backtracking cusps the cusp anomaly and the perimeter divergence mix together [|l^]. 
Zig-zag invariance can be implemented in the MM loop equation only by means of a subtle, reg- 
ularization dependent, cancellation between the extra perimeter due to the backtracking cusp and 
the cusp anomaly. This is due to the non-cohomological nature of the contact term that arises as 
the obstruction to locahzation in the MM loop equation, that is the line integral of a non-integrable 
distribution. 

We can write the MM loop equation also for a planar twistor Wilson loop: 

i "-l^-^fe-'-^^'o- 

if dzj dz5'^'^\z-w)<Tr^>{w,z;B)><Tr^>{z,w;B)> (10.13) 

"/ L.WW ^ i-'ww 

Also in this case the obstruction to localization of the loop equation is of non-cohomological nature, 
being the line integral of a 5^^', as opposed to the 5^'' that occurs in the regularized holomorphic 
loop equation. For twistor Wilson loops the contact term in the MM loop equation gives rise to the 
same divergent contribution ^^ for backtracking cusps as for unitary Wilson loops. Cancellations 
occur only in the solution of the loop equation. But even if it were possible to get cancellations 
at the cusps by fine-tuning, a lattice regularization of the MM equation would provide links rather 
than points, because the functional integral in the MM equation involves the connection instead of 
the curvature. 

Hence the MM loop equation cannot localize in the homological sense that we are discussing, 
not even for twistor Wilson loops, because it cannot be regularized in a way that does implement 
the zig-zag symmetry. 

11. Canonical beta function of large- A^ YM by homological localization 

11.1 Gluing rules for local systems 

In this section we compute the canonical beta function of the large-A/^ YM theory in the scheme 
defined by the homological localization of the loop equation. We take into account some global 
constraint that we have disregarded in our computation of sect.(9) of local counterterms in the 
Wilsonian scheme. 

We have seen in the last section that, for getting homological localization of the holomorphic 
loop equation, we should draw backtracking strings from the loop to the lattice points in order to 
transform all the marked points into cusps. 

We now show that the cusps can be paired by the backtracking strings. 

This is a consequence of requiring the consistency of the gluing rules for functional integrals 
with the localization on local systems. Indeed, to glue together two disks with twistor Wilson 
loops as boundaries, we need that the twistor Wilson loops on the boundaries, taken with opposite 
orientation, agree. Moreover, in a theory in which the twistor Wilson loops are localized on local 



The missing factor of A is in fact hidden in the normalization of 5/vc (sect.(4)). 
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systems, they are determined by the holonomies around the punctures inside the disks, since local 
systems are the same as representations of the fundamental group. 

Therefore the precise condition for gluing is that the product of the holonomies around the 
points inside the two disks are the inverse of each other, in such a way that the total product is 
1. This constraint is satisfied naturally gluing a lattice of vortices in one disk with a lattice of the 
same number of antivortices ^^ in the other disk. Hence in this case the number of punctures in the 
two disks to be glued must be equal in general. Another possibility is to pair vortices in one disk 
with vortices in the other disk and to impose the condition that the product of all the holonomies 
is 1 after gluing, for example requiring that the total number of points is an integer multiple of N. 
This can be certainly done for disks that have the same number of punctures and the same cutoff. 
Afterward one of the disks is rescaled together with the cutoff to a large size in order to perform the 
thermodynamic limit (see below). The procedure does not spoil localization since all the punctures 
are paired by backtracking arcs, that can be glued on the common boundary of the disks since the 
punctures are equal in number. 

This looks like a holographic correspondence [^] that is the consequence that all the punctures 
can be paired by arcs connecting punctures in different disks ^^. 

The resulting configuration has not a translational invariant cutoff, as in the computation of 
the Wilsonian beta function, but what really matters for that computation to hold is that the lattice 
on which the background field lives has a very large number of punctures and it is locally uniform. 

If the punctures of the two disks are paired, the density of the punctures in the two disks cannot 
be the same in the thermodynamic limit, since the area of one disk must be much larger than the 
area of the other one. Thus we get two cutoff scales, a and a, on the two disks, that we refer to as 
the ultraviolet and infrared scale respectively and that we identify with the lattice spacing. Thus 
each puncture carries a weight that is the lattice scale. Since the weights are all equal inside the 
two disks, because we require at least locally uniform lattice spacing, and the punctures are equal 
in number, the two collections of weights are in fact projectively equal, i.e. equal up to a common 
rescaling factor. 

Therefore the localization of the loop equation for twistor Wilson loops on local systems leads 
to a kind gluing that coincides with the one implied by a weighted arc families with projectively 
equal weights ^]\\. Surprisingly these are precisely the string gluing axioms [61, ^, ^3|] for 



weighted arc families at topological level (fig.4 I of []63|]). Hence we may say that open strings 
solve the YM loop equation for the twistor Wilson loops, in the sense that they localize the loop 
equation on a saddle point for an effective action. 

We have just come to the conclusion that if we combine the gluing properties of the functional 
integral with the vanishing requirement for the contact term in the holomorphic loop equation, i.e. 
the requirement of localization, we get the string gluing axioms at topological level. Thus the 
proper graph to get localization is a weighted graph [ ]6l1 ] similar to a Mandelstam graph [64, 65] 
(fig.4 I of []63|]). The weights in this language are the sizes of the strips. 



^^By antivortices we mean surface operators tiiat have precisely the inverse holonomy of vortices. They have the 
same classical action, if evaluated at the same cutoff scale, and a neighborhood with moduli space of the same dimension 
as for vortices. 

^^These are the "hairs" associated to the holographic correspondence between the ultraviolet and the infrared (see 
below). 



81 



Glueballs in large-N YM by localization on critical points Marco Bochicchio 



A subtle point arises about the cutoff of the locaUzed effective action. We recall that the 
introduction of a lattice is essential for localization, since it allows us to transform every non-trivial 
marked lattice point into a backtracking cusp. We have seen that the gluing rules imply a different 
cutoff at the cusps of the two disks of the sphere in the quantum effective action. 

Hence the local part of the effective action, as a consequence of the stringy nature of local- 
ization in the loop equation, is in fact bilocal with two local fields living at different scales, one at 
the ultraviolet cutoff and one at the infrared cutoff, paired by the backtracking strings. We will see 
momentarily that the field at the ultraviolet, but not the one at the infrared, affects the renormaliza- 
tion of the Wilsonian coupling constant, as it is expected from its very definition. Instead the field 
at the infrared together with the one at the ultraviolet affects the renormalization of the canonical 
coupling constant. 



To summarize, we draw our weighted graph \ pl[ p3| , |64), |65|], that is made by two charts with 
the boundary loop in common. The charts are a conformal transformation of two topological disks 
with punctures (fig.(l) and fig.(2) of ||6^]). This introduces a conformal transformation in the EK 
two-dimensional reduced theory. However, this transformation lifts to a conformal rescaling in 
the four-dimensional parent theory because on the Lagrangian submanifold that is the support of 
twistor Wilson loops Eq.(10.6) holds. Thus the four-dimensional metric changes conformally and 
the effective action changes by the appropriate conformal anomaly (Eq.(10.7)). This implies that, 
in addition to the explicit cutoff dependence, the effective action on the weighted graph is related 
to the one on the sphere with marked points by the addition of a divergent conformal anomaly, 
because of the singularity of the conformal transformation. This divergent conformal anomaly 
plays a key role to reconcile our computation of the anomalous dimension in the canonical beta 
function with the general properties of the RG group and, more generally, in the interpretation a 
posteriori of localization as a RG flow to the ultraviolet. 

11.2 Z factor and canonical normalization 

It is useful to write the effective action of the commutative Morita equivalent theory before the 
EK reduction, that amounts to dividing by the factor of A'^2 (sect.(9)). This is convenient to define 
the canonical normalization of the effective action in analogy to the SUSY case of sect.(3.2). 

The contribution of one surface operator of Zn holonomy, that we call one- vortex, to the local 
part of the Morita equivalent Wilsonian effective action reads: 

exp {—rq{one — vortex)) 

, 1% 871^ 271 k(N-k), , 1 ,, 

= '="'<-7?:?iAZ-'MA'-*))exp(^-4-Jlog(^)) 

, 271 871^ ^ ^^ ,1% k{N-k) , 1 ,, 

where the factor of ^^ is the transverse measure over the zero modes two-dimensional vortex 

sheet that is the singular divisor of a surface operator in four dimensions, and // = ^ by definition. 
The factor in the second line of the right hand side is actually the contribution of the antivortex 
in the infrared ^^, that is paired to the vortex by a backtracking string. Equating the ultraviolet 



Or of the vortex in the infrared. 
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cutoff on the transverse and longitudinal planes, because of rotational invariance, we get a = aj- 
Another way of formulating this condition is that the longitudinal measure on the size of vortices 
that we read from the classical action and the transverse measure on the vortices two-dimensional 
sheet should coincide. The Wilsonian beta function of sect. (9) follows, since the contribution of 
the antivortex at the infrared is finite. 

For a twistor Wilson loop in the adjoint representation it is possible to pair the holomorphic 
integral to the antiholomorphic one. In this case the counting of zero modes corresponds to the real 
dimension of the orbits. Thus we get for the bilocal part of the Wilsonian effective action: 

exp(-r,)=ne.p(-|lgz-..„(«-W+....)exp(i|^^(^log(^) + ....) 

We now come to the canonical coupling. We observe that the fields at the ultraviolet and at the 
infrared are not canonically normalized in the Wilsonian effective action. In order to obtain the 
canonical /3 function, as in the ^ = 1 SUSY YM case of sect.(3.2), we have to rescale the fields 
on the ultraviolet divisor by a factor of gZ^ and on the infrared divisor by a factor of g, since Z 
is finite in the infrared and can be normalized to 1 . This can be achieved by rescaling a and H as 
a = UcZ 2 and H = g H(. and setting a^ = aj to preserve an equal longitudinal and transverse 
cutoff. Rescaling H is equivalent to rescale inversely L^, the area of a surface operator. The 
canonically normalized effective action reads: 



exp(-r,) =Y{cM-Tr-37r-rs%{N-kp)+cx)^M—-2s' ' . iQg( ^, -27-1,2 )+^-^-) 

p nc-u^ z.g^ ncU^ I. ricg ^ "c 

«P(-|?.gr*.(«-*,)+«Oexp(||r^«^.o.(j^)+<-.<-0 (U3, 
This produces the following rescaling factors at each point from the contribution of the zero modes: 

Thus defining as in sect.(3.2) 

and factorizing the term %n^k{N — k)g^ we get: 

114 1 

+ -— ^log^ + -— ^logZ (11.6) 



IgliA) 2gHA) (47r)2 ^* (47r)2 
Taking the derivative with respect to log A and using the fact that gw is 1-loop exact we obtain: 

1 dg 4 \ dg 1 5logZ 



§3 5 log A (4;r)2g5logA (47r)2 5logA 
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from which it follows: 

^3 , g' dlogZ 



dg _ h>Og -r J^ c)logA 



Since 



^1°8^ 1-W^ 



^logZ 10 gl, 1 10 gl 10 ^2 



(11.8) 



aiogA 3 (47r)2i , 10 gi .A 3 (47r)2 3 (47r)2 



(11.9) 



we get: 



in „5 4 

(-A)§' + ^7l;^4)(l + 77:;;^2^' + •••) = -/3o^'-^l^' + ■•■ (11-10) 



where 



5 log A ^ '^'^ 3 (471)4^^ (471) 

o 4/3o 10 1 34 1 



(471)2 3 (471)4 3 (471) 



(11.11) 

that agrees with the perturbative result up to two-loops \ pl\ |32| , |33| , [34| ]. 

Now we come to the computation of the anomalous dimension. The one-loop exactness of Z 
implies: 

-.. ^ 1 10^2 



^loga l-g^74^flog(^ 



(4;r)2 

where now we have included the contribution of the conformal anomaly, that rescales the subtrac- 
tion point by the factor of a to give a finite but arbitrary result for the higher order contributions 
to the anomalous dimension. Assuming that the beta function is independent on the subtraction 
point, as required by general principles of the RG, the RG trajectory must be followed along the 
line c = — TT^-ylogl^) = const. We observe that it is precisely the contribution of the conformal 
anomaly that allows the anomalous dimension to be a function of the coupling only, according to 
the RG. 



12. Glueball spectrum from the localized effective action 

12.1 Glueball propagators and anomalous dimensions in perturbation theory at large-A'^ 

We start this section recalling some features and conjectures about the large-// limit of pure 
YM. 

To the leading \aige-N order and to every order in the 't Hooft coupling constant g, the ex- 



pectation value of a product of normalized local gauge invariant operators factorizes [ ]119| ]. For 
example: 

< iL^'-^a/3(^l)-4l^'-^a/3(^0 > = < iL^rPapi^O > - < ^L^''^a/3 (^^) > (12-1) 
^^ aP ^^ ap ^^ ap ^^ ap 
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Thus to this order the only information that survives is the value of the condensate. In the case 
of < g^^T^ap ^^''^afi(-^) >' ^^ must be proportional for a suitable regularization to the appropriate 
power of the renormalization group invariant scale, Aym, since it coincides up to numerical factors 
with the action density. 

In turn Aym encodes the information on the beta function of the laige-N theory. In addition 
it is believed that to the next to leading ^ order the connected two-point functions of local gauge 
invariant operators are saturated by a sum of pure poles. 



For example, for the scalar glueball propagator it is conjectured that the equation holds [119, 



l<^L^rF'api^)i:^trF^apiO)>conn ^'^^^ = I "^f^ (12-2) 



ap ap 



The sum of pure poles is constrained by the perturbative operator product expansion [11]. It must 
agree asymptotically for large momentum with the "anomalous dimension" of the glueball propa- 
gator as computed by perturbation theory plus the sum over condensates that occur in the operator 
product expansion [[77|]. Indeed the scalar glueball propagator behaves in perturbation theory at 
large momentum, within two-loop accuracy, up to contact terms, i.e. polynomials in the momen- 
tum squared, p^, as [ f77| ] : 

/(p)/log(^) (12.3) 

The logarithm explicitly displayed in the two-loop computation is necessary to reproduce the 
conformal behavior: 

/c/V>Vlog(^)~^ (12.4) 

The factors of g, the renormalized 't Hooft coupling at momentum p, occur because of the canon- 
ical normalization of the glueball propagator, that involves < Y,ap i^^^'^afiC-^) > rather than the 
action density < g^'^Y.aB ^^^^api^) >> ^^'^ ^^^y account for the one-loop anomalous dimension 
that in this case is determined by the one-loop coefficient of the beta function. They imply log- 
arithmic corrections to the conformal behavior. Glueball propagators for other normalized gauge 
invariant operators of naive dimension L involve in general one-loop anomalous dimensions that 
are independent on the one-loop coefficient of the beta function [^, 21]: 

j < -trO{x)-trO{0) >conn e'P'd^x 






Go{p^)-Zl{p^)p^^-Hog{^) (12.5) 



and satisfy the following RG equation: 

d „, . d 



(37— +^(5)^+27ofe))Go(/) = (12.6) 

dXogp dg 
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where 

ro{g) = 37^^ (12.7) 

Sometimes it is convenient to factorize out the contribution of the anomalous dimension: 

G'o{p^)=Zo\p^)Go{p^) (12.8) 

in such a way that the suitably normalized glueball propagator, G'q{p^), is RG invariant ^^: 

' ^ +^(^)i-)G'o(;^2)=0 (12.9) 



dXogp 
In the large-A^ limit there is a sector of the theory that is integrable at one loop [2^ 21], that is made 



by operators of ASD or SD type and their covariant derivatives. The corresponding anomalous 
dimensions can be computed as the eigenvalues of a Hamiltonian spin chain. 

In the ASD one-loop integrable sector the anomalous dimension of ^ Y.a^ '^^^aB (^) ^^ ^^^ 
same as the one of ^ l^^p TrF^p (x) since < g-^ Y,ap jj Tr{FapFap ) > and < ^"^ ^^^ ^ j^^2^ (^) > 
are both RG invariant pO|]. Therefore the renormalization factor of ^^„^ '^^^aB^^^ ^^ ^^ ^'^'^ ^^^ 
one-loop anomalous dimension coincides with — 2j8o. 

The anomalous dimensions of a number of operators can be computed explicitly solving by 
the Bethe ansatz the Hamiltonian spin chain in the thermodynamic limit, that corresponds to op- 
erators of large naive dimension 2L and length L^^. In particular the anomalous dimensions of 
the antiferromagnetic ground states turn out of to be of the form (Eq.(27) of [|^] and Eq.(5.23) of 



[21]) 



Zi = 1 -Lg2^^1og(^) + 0(L°) (12.10) 

The ground state of the spin chain corresponds to the operators with the most negative anomalous 
dimension, that turn out to be all scalars constructed by certain contractions involving only the ASD 



part of the curvature [21] 



12.2 Localized effective action and holomorphic/antiholomorphic fusion 

We are now ready to start our computation of the glueball spectrum. The large-A'^ one-loop 
integrable sector contains the correlators that can be computed by localization. This is perhaps not 
completely surprising since localization involves the change of variables from the gauge connection 
to the ASD part of the curvature, that lives in the one-loop integrable sector. 

In the previous sections we have obtained the localization of trivial twistor Wilson loops, in 
the sense that we have reduced the loop equation for them to a critical equation for an effective 
action and we have renormalized the effective action in Euclidean space, in such a way that the 
effective action is finite. 



^^We learned this suggestion by an unpublished talk of Gabriele Ferretti, "Applying the BES trick to QCD", Newton 
Institute, Cambridge (2007). 

^**We would like to thank Konstantin Zarembo for explaining to us this result at the GGI. We would like to thank also 
Gabrielle Ferretti for explanations about the same subject at Chalmers University. 
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Now the question arises as to which infrared information the renormalized effective action 
actually contains, since the effective action is naturally defined on a physical section of space-time, 
as opposed to the support of the trivial twistor Wilson loops. 

We have seen that to make sense of the holomorphic loop equation of sect.(6) after the renor- 
malization of the Euclidean effective action it is necessary to continue analytically the Euclidean 
loop equation to Minkowski space-time. Then at operator level the twistor Wilson loops become: 

Tr^W{Bp;L^,^,) -^ Tr^Pexp/ f {A,+ + ipDu)dz+ + {A,_+ip-^Du)dz- (12.11) 

as shown in sect.(4). The support of the twistor Wilson loops analytically continued in this way 
becomes: 

{z,z,u,a) = {z+,z-,-pz+,-p^\-) (12.12) 

For computational simplicity we set p = — 1 in the following. This is just an irrelevant change of 
sign of the antihermitian part of jXi, since /Xp = /x" + p« — p^^n. Thus in this section /x = /x_i = 

It is in the spirit and in the letter of the localization idea that the renormalized effective ac- 
tion analytically continued to Minkowski space-time contains information on local observables 
supported on the Lagrangian submanifold displayed in Eq.(12.12) since this is the analytic contin- 
uation of the support of Euclidean twistor Wilson loops from which the effective action has been 
obtained by localization. However, getting the spectrum of fluctuations exceeds by far the lim- 
ited framework of localization of the homology of 1 in the same sense in which the statement that 
the prepotential of ^ = 2 SUSY YM determines the low energy effective action in the Coulomb 
branch exceeds by far the framework of the localization of the cohomology of 1 . 

We have seen in sect.(lO) that localization by homology requires implicitly extending diago- 
nally the action of the conformal group from two to four dimensions, an extension that can occur 
meaningfully only on Lagrangian submanifolds of the kind displayed in Eq.(12.12). Thus the 
aforementioned Lagrangian submanifold is in a sense the only object for which we can hope that 
the localized effective action has a physical meaning beyond the leading \aige-N approximation. 

More physically homological localization requires that the ASD field be singular and therefore 
by necessity be of magnetic type, in a theory in which smooth fields are of electric type. Thus 
homological localization realizes in a technical sense long-standing ideas on the YM vacuum as a 



dual superconductor Q16i |17| , |18| ], whose analog are the duality ideas [ ]4-3| ] that lead to the physical 
justification of the prepotential as the low energy effective action in the Coulomb branch. 

In particular 't Hooft duality leads us to hope that homological localization may capture the 
mass gap via the implied condensation of the magnetic charge. 

We write now the localized effective action in the Wilsonian scheme. From the effective action 
we extract the glueball spectrum restricted to the Lagrangian Minkoswki section of space-time. To 
compute fluctuations we need to extend the effective action from the fixed points, characterized by 

— — /s 

d^D = d^D = 0, to fluctuations in neighborhood of them. The fluctuations correspond generically 
to surface operators of semisimple type, i.e. with different eigenvalues of the (diagonalizable) 
holonomy and of the ASD curvature, /i. 
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For a twistor Wilson loop in the fundamental representation this extension can be performed 
only in a holomorphic way via the holomorphic loop equation. The corresponding effective action 
is not Hermitian. Therefore we couple the twistor Wilson loop in the fundamental representation 
with the one in the complex conjugate representation^^ . As a result the effective action is Hermitian 
^°. Thus we are realizing a sort of holomorphic/antiholomorphic fusion that is reminiscent and in 



fact technically very similar to holomorphic/antiholomorphic fusion [ p9| , pO| ] in conformal field 
theory. 

Initially, thanks to large-A^ factorization, we treat jj. and p. as independent variables that define 
two different chiral sectors. But then we choose the section in function space where p. is actually 
the Hermitian conjugate of /i for the effective action to be Hermitian. This enables us to compute 
jj fluctuations. 

Yet, this n/p. sector is only a special sector of the theory, that involves certain correlators 
constructed by /i = ^(/^qi "'^03) ^^'^ ^^^ Hermitian conjugate. These are the special correlators ac- 
cessible to homological localization of the loop equation. We cannot say anything about correlators 

ofMo2- 

Realizing holomorphic/antiholomorphic fusion in the non-commutative Eguchi-Kawai reduced 

theory (sect.(4.1)) involves therefore the product of the holomporphic and antiholomorphic parti- 
tion functions, extended to the holomorphic and antiholomorphic neighborhoods of the fixed points, 
that are thought to be each the Hermitian conjugate of the other one: 

p 

= 1 f5AdAdDdDd{-iFs-Y,flpS^^\z-Zp)-d-h) 
J p 

p 

^w{-^-^l,tr^frf,{P,Pp))Det{^)Y{5p'X (12-13) 

Sw p "M p 

We then proceed as in sect.(7). In the reduced theory the classical action is finite on surface oper- 
ators because of the rescaling explained in sect.(4.1). The U{N) non-commutative theory can be 
compactified on a torus of large area, L^. For rational values of the dimensionless non-commutative 
parameter: 

2%dL-^ = — (12.14) 

the non-commutative L'^(A'^) theory is Morita equivalent to 2iU{N y<N) commutative gauge theory 
on the smaller torus, L^N^'^, with non-trivial 't Hooft flux through the commutative torus and with 



^^This coupling is automatic via large-A' factorization if we start with twistor Wilson loops in the adjoint representa- 
tion. In this representation the condensation of the magnetic charge via surface operators of Za; holonomy is compatible 
with triviality of twistor Wilson loops, because of the pairwise cancellation in the holonomy of every Za; factor by its 
complex conjugate for any shape of the loop. 

^"The complex conjugate field and the Hermitian conjugate contain the same information since they differ by trans- 
position of the indices. 
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the same 't Hooft coupling constant, g. Thus the Iwge-N limit of the partition function reduces to 
a finite dimensional inductive sequence: 

Z = lim I 1 8A5A8D8D8{-iFB-yu{u,u)iXpU{u,uy^8^^\z-Zp)) 

8{dAD-dAD + iY^{U{u,u)lJ.pU{u,uy^ -U{u,uy^p.pU{u,u))8^^\z-Zp)) 
p 

8{dAD + dAD)exp{-^Y^trNTrfj{^pfip))^ll8^'pf (12.15) 

The unitary matrices, U{u,u), account for the twisted boundary conditions on the commutative 
torus. In the thermodynamic limit, L^ — )■ oo, the effect of the twist of the torus disappears since the 
Morita equivalent theory must coincide in the Xwge-N limit with the commutative theory on R^. 
The equivalence with the untwisted theory can be seen also choosing the gauge where there is no 
twist in front of the delta functions. Of course this gauge is singular on the torus, in the sense that 
it is defined by a gauge transformation that is not single-valued on the torus, precisely because it 
is not continuous on the boundary of the fundamental domain that is employed to define the torus. 
But the boundary becomes irrelevant in the thermodynamic limit. 

We require that the fields on the commutative R^ so obtained in the thermodynamic limit have 
well defined limits at infinity, in such a way that the theory can be compactified on S^. Thus the 
fields on R^ x R^ can be extended to S^ x S^ and the resolution of the identity reduces to the one 
for ordinary surface operators. Thus we get: 

Z = lim I f 8A8A8D8D8{-iFB-yiXp8^^\z-Zp)) 

N,N^oo J p 

8{dAD + ~dAD)8{dAD - ~dAD + /£(Aip - jXp)8^^Hz - Zp)) 

p 

exp{-^l^trMTr^i^pfLp))^ll8^'p\' (12.16) 

Sw p ^M p 

that in the Feynman gauge, as shown in sect.(3.3) and sect.(8), reduces to: 

Z = \J[Det-'/\-AA8„p-i^~-^)Deti-AA)]^^.^^^-^_,^^ 

A-[^'l exp ( -'-^Y^tr^Tr.i^pfip)) ^ A Co'^-^U^^'pf (12.17) 

Sw p "t^ p 

In this formula ft)' should be chosen in a way compatible with the holomorphic/antiholomorphic 
fusion. In particular co' should depend holomorphically on the complex eigenvalues. A, of /x. This 
requires that co' be the puUback on the moduli space of the symplectic form associated to the twistor 
connection ior p = —I, B = B \ (sect. (7)): 

1 r' 
(o'=— d^ztrffr{8B^ A 8B^) (12.18) 

The superscript in co' refers to the version of (O defined on the punctured sphere as discussed in 
sect(7). This is not restrictive, since excluding the singular divisor is equivalent to omitting in (O 
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the sum of the Kirillov forms of the adjoint orbits on the singular divisor [120], because the volume 
form on these orbits is already taken into account by the product measure on the adjoint orbits, 
Inp^Mpl. (sect.(12.3)) that occurs by the resolution of identity. As a consequence: 

(0^ A6)^ Alld^'pdfi;, = co'^ A 6)'^ Alld^'pd fi'p (12.19) 

p p 

In any case the version that involves co' is the one that leads to the glueball potential. 

12.3 The glueball potential 

The glueball potential arises as a term in the localized effective action that is the logarithm of 
the modulus of the Jacobian for the change of variables of the complex field of ASD type from a 
unitary, /i, to a holomorphic gauge, /i'. 

This term cannot arise in perturbation theory and the only way to derive it in our approach is 
via the holomorphic loop equation, since the choice of the holomorphic gauge is necessary in order 
to produce the Cauchy kernel in the right hand side of the holomorphic loop equation. 

In turn the Cauchy kernel is essential, because after analytic continuation to Minkowski space- 
time leads to localization by homology, since the corresponding regularized distribution is zig- 
zag invariant in a neighborhood of a cusp and therefore its contour integral vanishes for arcs that 
backtrack at the cusps. 

The physical meaning of this Jacobian is the following. 

For surface operators, because of the specific features of the Hitchin equations, there is a 
mismatch between the degrees of freedom carried by the ASD curvature at a point in a unitary and 
in a holomorphic gauge. 

While at moduli level the unitary and holomorphic gauges are completely equivalent as shown 
in sect.(7), the local degrees of freedom that are manifest looking at the ASD curvature in the 
unitary and in the holomorphic gauge do not coincide completely. 

In the holomorphic gauge there is an essentially '^ one-to-one correspondence between the 
moduli that occur in the holonomy, Mp, of the connection, B, around a point, p, and the moduli of 
the local curvature, jj.' at p, given by the equation: 

Mp = e^'^'" (12.20) 

Thus in the holomorphic gauge jj.' is parameterized by orbits in the complexification of the gauge 
group, Up = GpXpGp^, with Xp the (complex) eigenvalues of n'p, that we assume all different in the 
semisimple case relevant to compute fluctuations. 

Therefore the integral over jj.' = GpXpGp ' is on an orbit of the complexification of the gauge 
group with measure: 

5ix'p=Det{adXp)8Xp8Gp (12.21) 

where the integration on 5Gp is actually on the moduli of the orbit, that can be parametrized as 
Gp = gpPp, by factorizing the complexification of the gauge group into its compact and parabolic 
factor ^^. In fact it is not restrictive to project the parabolic factor, Pp, to its unipotent subgroup, P' 

^'The moduli depend actually only on the conjugacy class of the eigenvalues of the holonomy, that determines the 
eigenvalues of the curvature only up to shifts of 2% . 
^ This is the Iwasawa decomposition. 
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since the diagonal factor of Pp acts trivially by conjugation on the eigenvalues. 
Thus the preceding equation can be rewritten as: 

d^'p = A{XpfdXpdgpdP'p (12.22) 

where the square of the Vandermonde determinant of the eigenvalues, A(A): 

A{X) = Y[{^i-Xj) (12.23) 

•>j 



arises by the usual Faddeev-Popov procedure, as in the holomorphic matrix models [£! 

In a unitary gauge this does not hold true. It is a fundamental result of the theory of sur- 



face operators [85, |8^ that in a unitary gauge [Hp,jlp] = 0, so that generically jXp and jlp can be 
diagonalized simultaneously by a unitary gauge transformation (sect.(7)). Thus only the moduli 
associated to the adjoint action of the compact unitary group, gp, are manifest in the ASD curva- 
ture in a unitary gauge, while the remaining degrees of freedom are hidden in the moduli of the 
Hermitian metric whose choice is implicit in the unitary gauge [^]. 

Now the non-Hermitian matrix. Up = gp{Xp + u' )gj,^ , that occurs in the resolution of identity 
that defines the non-SUSY Nicolai map in a unitary gauge, is generically conjugated by a unitary 
transformation to a triangular matrix, Xp + u' with u' triangular and nilpotent, in such a way that 
the induced measure is: 

^Mp =Pf{ad^ip)8Xp5u'p5gp 

= A{Xp)5Xp5u'p5gp (12.24) 

The different power of the Vandermonde determinant in Eq.(12.24) with respect to Eq.(12.22) 
arises by the Faddev-Popov procedure for triangularizing rather that diagonalizing /ip ^^. 

But the unitary gauge transformation that sets the complex matrix /ip in triangular form in 
a unitary gauge does actually diagonalize those Up that arise as the ASD curvature on the dense 
locus of the surface operators in the resolution of identity of the non-SUSY Nicolai map. Therefore 
in the induced measure the integration over the nilpotent part, u'p, actually decouples, since the 
ASD curvature in the resolution of identity does not depend in fact on u'p for those connections 
that arise by surface operators in a unitary gauge. Indeed the delta function in the resolution of 
identity, 5{-iFB-'£plJ.p5^^\z- Zp)), reduces to 5{-iDB a5B -Y,pgpu'pgp^d^^\z-Zp)) around 
the solution of —IFb — T.pSp^p8p^S^^\z — Zp) = and of the remaining Hitchin equations. But 
then the last delta function implies u' =0 because for smooth fluctuations, 5B, the argument of the 
delta function is dominated by d ^^^ (z — Zp) at each point Zp- Therefore the integral on u' decouples. 
On the other hand the integration over the moduli, P' of the twistor connection, B, contributes to 
the zero modes in the unitary gauge and produces the measure: 

(o'''^=Pf{co')lldP'p (12.25) 



^ We may get the same result noticing that the integration measure dfldfl on normal non-Hermitian matrices is 
\A{X)\ dXdXdg by the standard Faddeev-Popov procedure and then taking the "square root". 
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in such a way that: 



8^1 ^^ri^i =/./(«') j-r.^^^^^^p) ^^pK^Sr. 



8}x' ^p^Det{adjJ.p) 8Xp8Gp 

Pfjco') 
UpPfHiXp) 



(12.26) 



where we used Gp = gpPL. The same resuh is obtained in the singular gauge where jXp is actu- 
ally diagonal, jXp = Xp, for which therefore 8u' 8gp decouples from the measure 5/Xp, since the 
delta function in the resolution of identity, 8{—iFB — Y,p I^p8^^\z — Zp)), reduces to 5(— /Dg A 
8B — £p u'p8^^^ (z — Zp)) around the solution of —iFs — Y^p Xp8^^^ {z — Zp) = and of the remaining 
Hitchin equations. Therefore the last delta function implies u' = and therefore the integral on u' 
decouples. But the measure on the moduli is now: 

(o''"^ = Pfico') n 8P'p8gp (12.27) 

p 

Therefore the effective action reads: 

r = ^Y.frNTr^{^pflp)+'£\og\A{^p)\^ 
Sw p p 

-log \Der'/\-A-^8ap - map)Det{-AA)\l=i^(^^^^^i^,^^) 

-2nh iix'] log A - log |P/(«') I' (12.28) 

We anticipated that the second term turns out to be the glueball potential, that generates the glue- 
ball masses. The term log \Pf{o)')\^ is irrelevant for the glueball potential, since co' depends on the 
eigenvalues only through the holonomy of the connection (sect.(7)). Indeed, since the holonomy 
at the critical points lives in Z]^, the second derivative of \og\Pf{(o')\^ at the critical points cou- 
ples only to the trace of the fluctuations, Tr{8Xp), and thus decouples in the \wg&-N limit. The 
Pauli-Villars factors do not contribute to fluctuations. The remaining terms are considered in the 
following subsections. 

It is very instructive to reinsert the factor of N2 to recover the unreduced theory and to intro- 
duce a regularization by the density of the lattice of surface operators that is more suitable for the 
continuum limit. Reinserting the factor of A'^2 we get: 

r = -^N2Y,trNTrf^{Hpiip) + Y,N2\og\A{iip)\'' 
8w p p 

-log|DerV2(_A;,5„/3-/M;^)£»ef(-AA)|^=.(^-__,-^-3) 

-2N2ni,[pi'] log A - A^2log \Pf{co')\^ 
SNN 



-^5(2) (0) / d^uY,trNTrfj{iXpfip) + Y,8\0) / d^ulog\A{iXp) 



|2 



■log \Det '/H-^A8ap - /M«^)Def(-AA)|^=.(^-,_,-^-^) 
-28(^\0)Jd^uY,nh[^'p]logA-8(^\0)Jd^ulog\Pf{(o')\^ (12.29) 
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with the traces that define the functional determinants properly rescaled. We can now introduce the 
density of lattice points: 

p=£5P)(z-Zp,) (12.30) 



normalized in such a way that 



JdhZ5^^\z-Zp')=N!, (12.31) 



is the number of lattice points at the scale at which the density is p . Notice that the density is not 
normalized necessarily to N2, the number of lattice points at the cutoff scale, because the primed 
sum is only on lattice points where the holonomy of the surface operator is non-trivial. This allows 
p to scale non-trivially with the RG. Assuming translational invariance at N = 00 the effective 
action reads: 

SNN 



Sw 



(fu(f-zp^trNTrf^{lX}x) + / J^Mj^zp^ log |A(jU)p 

-log|Deri/2(_A;,5„/3-/M„^)£>eK-A;i)|J^.(^-^_.^-^) 

-2 fd^ud^zp\,[n']logA- f d^uplog\Pf{(o')\^ (12.32) 

This form of the effective action is of the utmost importance, because it shows that the coefficient 
of the glueball potential is in fact a RG invariant scale as it should be. 

12.4 The effective action is degenerate at the fixed points 

We can compute p in terms of A\y, the RG invariant scale in the Wilsonian scheme, by min- 
imizing the renormalized effective action as a function of p for a given A and gw- In fact we 
know already from sect. (8) that the effect of the third and forth term in the effective action is to 
renormalize the coupling constant. Hence the local divergent part of the effective action for surface 
operators with Z/v holonomy of magnetic charge k and density p reads: 

k(N-k)N^(47i:f , 2 10 1 , A^ /" ,2 ,2 2 
-2k{N-k)N^ / d^ud^zp^log — + ... 



Ae Posw r 
k{N-k)N^{4nf{-po^og — - — ) J d^ud^zp^ + ... 



-jSoA^log — — — / d^ud^zM^ + ... 

-A)A^log ^ j d^ud^zM'^ + ... (12.33) 



93 



Glueballs in large-N YM by localization on critical points Marco Bochicchio 

where we have chosen the subtraction point at the scale of the action density: 

M^ = Nk{N-k){4nfp^ (12.34) 

This condition ensures that all the sectors labelled by k are degenerate in the large-A'^,A'^ limit, in 
such a way that the renormalized effective action at the subtraction scale is large and negative and 
equal for all of them. In this case the trivial solution with magnetic charge ^ = is excluded since 
it has greater action and therefore there is condensation of surface operators with all the magnetic 
charges. Indeed the critical equation: 

|^.4Mnog^-M==«'(41„g^-l) = (,2.35, 



has solution: 



M 4 
in such a way that the effective action reaches its negative minimum: 



log^ = T (12.36) 



<I 



{fud^zM"^ (12.37) 



with M^ given by: 



M^ = e-^A% (12.38) 



This also means that for large N the square density of surface operators scales as \ and as ^. It 
implies also that the glueball propagators are a sum of pure poles as we will see momentarily. 

We may wonder as to whether finite terms, i.e. not A divergent, may affect this picture. It easy 
to see that terms proportional to k{N — k) simply redefine K^/. However it is not obvious that the 
contributions from all finite parts have this form. Thus finite terms may affect the dependence of p 
on k. 

In fact we have chosen the action density of the condensate, p^Trf^tr]M{X^\f^\fj), as the in- 
frared subtraction scale, M"' '^. This is the same prescription as for the Veneziano-Yankielowicz 



effective action 1 121 ] of ^ = 1 SUSY YM. This subtraction point implies that the renormalized 



action is the same in every sector of magnetic charge k, in such a way that all the Z^ magnetic 
charges "condense at once" with a renormalized square density that scales as p 

While this is the same prescription that occurs in the effective action of SUSY theories in 
fact its justification at fundamental level may imply a certain fine-tuning of the finite parts in the 
renormalized effective action. This fine-tuning is always possible if the surface operators of Z^ 
holonomy that occur at the fixed points are viewed as limit points in the closure of orbits with 
unipotent holonomy in the Lagrangian cone mentioned in sect.(7) and sect.(9). Indeed in this case 
the relative scale of \Pf{o)')\^ can be suitably adjusted for different k approaching the limit points 
where the nilpotent residue of the Higgs field vanishes. 



' Finite changes of the subtraction point affect the normahzation of the glueball propagator but do not affect the 
glueball spectrum. 
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This is essentially due to the fact that the moduli space of such orbits is non-compact and 
that we are suitably choosing the size of the neighborhood of the fixed points in order to satisfy 
certain conditions. Indeed every "compactification" is to some extent arbitrary ^^. On the contrary, 
the unitary orbits in the Lagrangian cone, that are the other orbits that saturate the beta function 
(sect.(7) and sect.(9)), do not allow such fine-tuning since they are compact. 

12.5 The kinetic term 

We are now ready to compute the glueball propagator. It is quite clear that the classical ac- 
tion cannot furnish the kinetic term for the glueballs since it is ultralocal. Therefore the kinetic 
term must be generated by radiative corrections around surface operators. This turns out to be 
the case for the fluctuations of Lagrangian-embedded surface operators analytically continued to 
Minkowski space-time. They can be obtained by diagonally embedded Euclidean surface operators 
continued to Minkowski space-time. 

The diagonal embedding can be described as follows. We choose the surface {z = u,z = d) 
diagonally embedded in /?"*. Since we have defined a lattice in the (z,z) plane, this defines a lattice 
also in the (m,m) plane by the diagonal map {zp = Up,Zp = Up). This lattice in the {u,u) plane has 
a set of dual plaquettes in such a way that the (m,m) plane is the union of the plaquettes. Now we 
define the function Zp{u, u) = u with domain the interior of the plaquette dual to p and analogously 
for the complex conjugate. We also define a lattice fluctuating field supported on the plaquette 
dual to p and locally constant as (m,m) vary in the support, 5lJ.p{u,u). 8}i.p{u.,u) is zero outside its 
support. We are now ready to do computations. We suppose that in addition to the translational 
invariant background of surface operators there are locally-defined fluctuating surface operators ^^ 
diagonally embedded in space-time: 

-iFB = Y,^i8^^\z-Zp) + Y,5^ip{u,u)8^^\z-Zp{u,u)) (12.39) 

p p 

Since the kinetic term for the glueball propagator must arise by the radiative corrections we examine 
the expansion of the functional determinants in one-loop graphs with multiple insertions of trees. 
In our case these terms carry multiple insertions of the background field and of the fluctuating 
field supported on the lattice of surface operators. We have seen in sect.(9) and sect.(ll) that 
the divergent parts, that contain the background field, determine the beta function. We are now 
interested in the finite parts to second order, that contain the fluctuating field. 

The justification is as follows. Every term of the loop expansion contains a trace in the adjoint 
representation and thus it is proportional to A'^^, that of course diverges for large N. However, the 
loop expansion is in fact an expansion in powers of the density p of surface operators. But since 
the density scales as N^^ only the leading quadratic term survives the double large-A'^,iV limit. 



^^In Nekrasov theory of cohomological localization a certain (to some extent arbitrary) compactification of the 
moduli spaces of instantons has to occur too. Such compactification turns out to be compatible with the Seiberg-Witten 
ideas on the electric/magnetic duality in ,yV = 2 SUSY YM. In our case we can argue similarly that we can choose the 
size of the neighborhoods of the fixed points to avoid that only a proper subgroup of Zf./ condenses, spoiling 't Hooft 
ideas on electric/magnetic duality. 

^ Both the background and the fluctuations are diagonal matrices in color space. 
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The spin contribution to the effective action is: 

NN'A ^ /• 

^ ly^^^^ ^i:;,:;;;^::;:;;;;;:;::,:;^ (12-40) 



2NN'4 ^ /• 2 ,2 trNTrf^,{5Hpiu,u)5jlp>{v,v)) 



plus the complex conjugate term that we add only at the end of the computation. The orbital 
contribution has the same structure and a different coefficient and sign in order to reproduce the 
divergent Z ^ factor when evaluated on the translational invariant background of surface operators 



77 . 



Indeed the orbital contribution to second order, up to constant overall factors, is '^: 

fd'xd''yTr{A,{x)d,—^A,{y)d,—^+A,{x)da—^Aa{y)d,—^) 
J [x-yy [x-yY (x-yY [x-yY 

/" ,4 ,4 ^ /x / N^fz-w) ^ . .2(z-w) ^ , ,2(m-v) ^ , , 2(m-v), 
= / d'xd'yTr{A,{x)-^ dMy)T d+Mx)^ T?^" 3^ 7^ rj 

J (x-yr (x-yY (x-yY (x-yY 

/7 9 T T 412 — ^1 4|m — vl 
d^zd^wd^ud\Tr{A,{x)A-,{y)- '-^—^ +A^{x)Aa{y)- ^^ ^) 

/■ 9 7 9 9 4 4 

~ - / d^zd^wd^ud^vTr{A,{x)A^{y)d-^d„— p^— ^^-^ +Au{x)Aa{y)dad^ 

J d( z — wr + M — vr)^ 



6(|z-w|2 + |m-v|2)2 "^ ^ "^''^ " ''6(|z-w|2 + |m-v|2)2^ 

f 4 4 
= - / d^zcfwd^ud^vTr{did„{A^{x)Ai{yY-r7, ^^ ^-^ + 4-5v(A„(x)Av-(y))— r^— -^) 

/" 2 2 2 2 4 4 

= - d zd wd ud vTr{diA^{x)d„Ai{y)— -^-— -j-. + duAu{x)duAu{y)-j. r^-—. ,272 

J b[\z — w\'^ + \u — v\'^Y 6(|z — "^1^ + |m — v|^)^ 

V /^2 ,2 NTr{8^lp8^lp>) 

~ > I d ud V-. -. r -. ,n I pttt; (12.42) 

p^,y (|Zp(M,M-)-Zy(v,V-)|2 + |M-v|2)2 

where in the last line we used diA^{x) ~ ^p5/Xp5(2)(2 — ^^(mjm)), 5j^A|(j) ~ — Lp^Mp^*-^H^ ~ 
Zp (v, v) ),duAu (x) ~ £p 5«p 5 P) (z - Zp (m , m) ) and 5vA„- (y) ~ - Lp 5«p 5 (2) (w - Zp (v, v) ) for the fluc- 
tuations of surface operators on the diagonal Lagrangian submanifold, with x = {z,z,u,u) and 
y = (w, w, V, v). The total result is: 

3(4?r2)2^'^^,i (\z^(u,u)-Zpiv,vW + \u-v\^Y 

The coefficient can be found without direct computation since, for the surface operators that are 
constant on the {u,u) plane and translational invariant, Eq.(12.43) must give rise to the logarith- 
mic divergence that produces Z^. Expressing Eq.(12.43) in terms of the density of the surface 



^^The Z^' factor, contrary to the beta function, does not depend on sitting on surface operators of Z^ holonomy. 
^ The term involving daAa vanishes identically around the local singularity, while the term involving A„ is quadrat- 
ically divergent and does not contribute because of cancellations due to gauge invariance. 
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operators, p, we get: 

20NN' /■ 2 ,2 2n /■ ,2 ,2 trNTrff,{5lx{u,a)5ll{v,v)) 



:{ j d^vd^up'^) fd^ud^v 



3(4712)2 V - ^- -r ;y - - (|„_v|2+|„_v|2)2 



3{4n2)^k{N-k){4ny'J wjj- -^ (|„_v|2 + |m_ v|2)2 

(12.44) 



that for large A'^ reduces to ^^ : 



20e ^ _^,2 f ^2^^2^irNTrf^,{5ix{u,u)5ix{v,v)) 



3(4;r2)22^(4;j)2 2y (j,_v)2(ji_v)2 



(12.45) 



M= d^uAl, (12.46) 



where we have set: 



This term generates the kinetic term of the glueball propagator after analytic continuation to Minkowski 
space-time ^^. Inserting the complex conjugate term we get: 

'^^"' ,N'' fdu^du^dv^dvJ'''^'^'^^^^"::"-^^^^'^^'-^}+c.c. 



?,{AK^Ylk{Anf ^J + " + - (m^_v+ + /£)2(m__v_+/£)2 

= ^'^^^^ 'i(4jt2)22k(4ny ^\l ^^+^^-^^NTrf,i{8iJ.{u+,u^)d+d^8}X{u+,u-)) + c.c. 

(12.47) 

We notice that for obtaining this result it is crucial that fluctuations occur as surface operators, that 
the support of fluctuations is embedded as a Lagrangian submanifold in Euclidean space-time and 
that the analytic continuation to the Minkowskian Lagrangian submanifold is performed. 

It is interesting to compare the kinetic term just obtained with the contribution of the classical 
action: 

SA^A^' / 2 2 2 

d ud zp trNTrfj,{8^{u,u)8^{u,u)) 



8% 



d zA^,) / d uAyftrNTrfji{5^{u,u)8fi{u,u)) 



glk{N-k){AKY 

(12.48) 



that at large A^ reduces to: 



N'2 I d uA^^trNTrfj,{5n{u,u)5lJ.{u,u)) 



(12.49) 



^^It is not restrictive to require k= 1 , • • • , y for A' even, because of the symmetry of all our formulae for the exchange 
k^N-k. 

It is necessary to assume that 5/l(^j (m+ , m_ ) and 8iiQ^{u+,u) are the boundary values of holomorphic functions 
on the upper-half plane for each of the independent variables, (k+ , k_ ) , with suitable properties at infinity. 
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Hence the classical action is irrelevant with respect to the term generated by radiative corrections 
for the fluctuations in the \wge.-N limit and in the thermodynamic limit in which N'2 diverges. In 
fact the real source of the mass term is the glueball potential whose contribution to the effective 
action is: 



f d^ud^zp^log\Ain)f 



1 



N'k{N-k){4nf 

^iV;/</^„A^,og|A(;,)| 



N'k{N-k){Ai 



(12.50) 



Apparently the glueball potential in the thermodynamic limit, N2 — ?■ °°, is as suppressed with respect 
to the kinetic term as the classical action is. Therefore the only possibility for the theory to have a 
mass gap is that the glueball potential is singular. In fact these singularities of the glueball potential 
arise precisely when some eigenvalues coincide. This looks encouraging because it implies that 
the mass gap may arise only by configurations for which the gauge group is unbroken. In addition 
these configurations for which the eigenvalues coincide must have infinite degeneracies in the ther- 
modynamic limit, in order to compensate the -X- suppression with respect to the kinetic term. The 

"2 

needed degeneracies follow by the ansatz that the center Zj^ of SU{N) occurs with multiplicity A'^ 
in the commutative theory Morita equivalent to the non-commutative one. 

12.6 The mass gap 

We are now ready to compute the mass matrix. With our normalizations e^'^ G Z^ with /x G 
su{N) and jx the translational invariant condensate. Thus: 

2lx=diag{2n{k-N)/N,2nk/N) (12.51) 

k N-k 

The mass matrix in the /x//l sector is the second derivative of the logarithm of the modulus of the 
Vandermonde determinant: 

J diXidiXj diXidiXj ^^ 

= ^T^ T, ( ^ai ^Hi) +C.C. 

= 4^{l^— I^— )+c.c. 

= 3^(E + E )+c.c. 

d ^ 1 
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^ L (5^^HjU<- - lip)5ij - 5^^\iii - ixp)djp) + c.c. 

Hi 
% X 5P)(^, - yip)5,j - 5(2)(^, - ^^.)(1,7 - 5,7) +C.C. (12.52) 

Hi 



where we used: 



i^l = ;r5 5(2)(^) (12.53) 

dZi Zj 

and 1,7 stands for a matrix with all entries equal to 1, l;y = 1, V ij. 

5^'^\}Xi — jJ-j) is a distribution in color space. It is nonvanishing when the two eigenvalues are 

degenerate. 

Let us now specialize to the Zj^ vortices in Eq.(12.51) with multiplicity A'^'. 

For the the first block, /, j = \,.. . ,k, Eq.(12.51) gives: 

Ml = n8^'-\0){{N'k-l)8ij - {hj - 5,-,)) = n5^'-\o){N'k8ij - hj) (12.54) 

with 5^^\0) = -^Tj the delta function at zero in color space. 

Since the diagonal terms scale as N', the non-diagonal corrections to the mass term are negli- 
gible in the large-iV' limit, in such a way that the theory has a mass gap ^^ 

Similarly, for the second block, i,j = k+\,...,^, one gets 

Mfj = 7r5P) (0) {N'{N - k)5ij - hj) (12.55) 

that implies that the glueball masses are at the cutoff scale ^^ for large A^ in this block. Thus the 
glueball mass term is: 

to be added to the kinetic term to get: 



■) 10 1 f 

(47r)2 ^ 3(4^^'^ J du+du^trNTr^,{d^d+d^dfi) + 

+-pr-^fj'^2 I du+du-7i:A^trNTrj^{5n5jl)) +c.c. 
^ {A%)H% ij^^2\l du+du-trNTrfj,{5^d+d^5fl) 
+N'N2 / du+du-A^trj^Trfj,{8jJL8jX)) +c.c. 



(12.57) 



Now we have the identification (Eq.(4.14)): 



In 
N{—f = 2Kd (12.58) 



^'For N' = I and k= I the theory has a massless eigenvalue in addition to the trivial diagonal U{1) that decouples. 
We would like to thank Daniele Dorigoni for working out this case during our course at SNS. 

It is not restrictive to require k= 1 , . . . , ^ for A' even since our formulae are symmetric for the exchange k^ N — k. 
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with 



2ne=L^^ (12.59) 

N 



and 



Therefore: 



A 



,2,2 



N2 = i—YL' (12.60) 

271 



^2n' N 



at the cutoff scale and thus for ^ — ;• ^ in the large A^ limit (sect.(7)) for any finite positive integer 



«>2: 



and at the renormalized scale: 



A2 = nN (12.62) 



A2 = nN' (12.63) 

by the diagonal embedding. We will see in the next section that any such choice of n corresponds 
to the choice of a different renormalization scheme for which the preceding equation holds. 
Finally, in the reduced theory in the sector labelled by k the glueball effective action is: 

g-i 10 1 /■ 
{4ny2n ^3n^'^k J du+du-trNTr^,{5^d+d-5fi) 

+N' / duj^du-A.Y^trNTrfj,{8jx8ji)^ +c.c. 

(12.64) 

12.7 Glueball propagators in the Wilsonian scheme 

Thus we find the following propagator in the twistor sector of the large-A theory for the Wilso- 
nian normalization of the EK reduced commutative Morita equivalent effective action: 

/I 1 

< — ^?rA,rr^,(jUjU)(x+,x_,x+,x_)— ^?rA,rr^,(jU/l)(0,0,0,0) >conn e'^''+''-+''-^+^^x+Jx_ 

1 °° Pa^ 

1 Y"! ft. iVjy 



N^N'^ti-(^'P+P-+^^l 



w 



~ T7^«"(-/'+/'-)'log-^ (12-65) 



^2^.2- ^ '^^'-' -^ A2^ 



with: 



a' = —n (12.66) 

371 
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Indeed it is not hard to see that, setting k^h^ = [(^A^ + a'p+p-){kl^ — a'p^p-) + {—a'p^p-)'^], 
the second line in Eq.(12.65) can be written as a logarithmic divergent sum that reproduces the cor- 
rect logarithmic behavior of perturbation theory: 



jf-'^ -a'p+p-+kl<^ 
_ - {{kAl + a'p+p-){kAl-a'p+p-) + {-a'p+p.f)Al 
hi -a'p+p-+kAl, 



a 



i-P+P- 



v2 



'Y ' '^^'-\ +..., (12.67) 

k^i-a'p+p-Ay^'' + k 



up to a divergent sum of condensates, proportional to a power of Aw, and up to a divergent sum of 
contact terms. We can now define glueball composite operators in the following way. Operators of 
the form: 

trNTrf^'i^-^^'afif (12.68) 

restricted to surface operators become: 

5(2) (o)2L-2 ^5(2) (0)5(2) (^ _ ^^)tr^Trf^, {^ipfipf (12.69) 

p 

They can be evaluated, in the same fashion as the action density, as : 

80-)(Qfi^-^p^trNTrf^,{^ifif (12.70) 

We set: 

A' = -Al (12.71) 

n 

in such a way that the factors ofN^^~^ from A^^^^^^^ are cancelled by the powers of A^^' from the 
denominator in trj^Trfi,5{}i.}i)^ evaluated on the condensate of surface operators. This condition 
defines a non-perturbative scheme in which the v.e.v. of composite surface operators are actually 
finite in the large {A,N) limit. The preceding relation implies: 

1 1 A^ 

/3olog- (12.72) 



IgUAw) 2gU^) 
and 

N'2 = nN' (12.73) 

in such a way that the rank ,yV = A'^2 of the Morita equivalent gauge group stays constant along the 
RG-^osN, but the degeneracy, A'^', flows into the rank, A'^, of the diagonal SU{N) along the /?G-flow. 
The flow starts with J/ =N2 = nN at the cutoff scale according to Eq.(12.62), that implies N = n 
at the cutoff scale. The flow ends with lower N' and larger A'^ at the renormalized scale in such 
a way that nN = NN'. Therefore in the large-A'^2 limit the gauge group in the infrared is SU {00) 
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embedded with oo multiplicity in U{°°), a result consistent with the structure of the gauge group of 
the non-commutative gauge theory [ |122| , |36| ]. 

We rescale also p^ by a factor of A^A'^' = N2 in order to get a quantity on the order of 1. Thus 
our operators are: 

5^^\ldf^-^N'NphrNTrp{pLpLf (12.74) 

The corresponding glueball propagators in the Wilsonian scheme are: 

~ i^f'^-'N^'-'N'^N^p' J < trMTr^,{d^^'^-'fi'-){x+,x^)trNTrf,,{^'^fl'^-'5fl){m >,onn e'^P-'-+P-'-U'x 



N ^y-a'p+p-+kAl, 

h\ -oc'p+P-+kAl, 

A^((/cA^ + oc'p+P-){l^^w - Gc'P+P-) + i-oc'p+P-] 



L 



a'p+p- + kAl, 



k=l ^ F+y-^i^^'^w 

'{-P^p-f'-'"^ 



i^i -a'p+p-+kAl^ 
-{-p^P-r-^log^P^ (12.75) 

up to a sum of RG invariant condensates and contact terms. Only the leading singularity for large 
momentum has been displayed. In fact the subleading singularities in powers of the momentum 
have divergent coefficients. We can "renormalize" these singularities, for which we have not an 
interpretation, as follows. Localization of twistor Wilson loops admits shifting the eigenvalues 
of surface operators by adding 27rl at every point ^^ . Indeed this does not modify anything in the 
SU{N) sector, but shifts the diagonal U{\) part of the action by a central term, that can be cancelled 
by a counterterm, in such a way that the effective action is flat for this L'^(l). 

We can now construct composite surface operators as in Eq.(12.74), but for the shifted curva- 
ture: 

llx=diag{lnk/N,2K{k + N)/N) (12.76) 

k N-k 

that has the same Z^ holonomy. Let us call the two dimensional Fourier transform of these opera- 
tors Ol{p-i-,P-). In the Wilsonian scheme we get: 



< Tr^(rip+,p-)Tr_^0'^{-p+,-p-) > 



conn 



r\j 



f ^2 ^2(2L-l)^4(2L-l) 



y 

^fi -a'p+p^ + kA"^ 



This shift may be related to the spht central extension that we disregarded in Eq.(7.15) as opposed to Eq.(7.1 1). 
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^h -a'p+p- + {N-k)Kl 

~(-P+P-)"-'l08^^ + ... (12.77) 

where the dots stand for contact terms, i.e. distributions whose inverse Fourier transform is sup- 
ported at coinciding points. But now there are not anymore subleading singularities in momentum 
with divergent coefficients. 

The natural interpretation is that this computation in the Wilsonian scheme furnishes the RG 
invariant version of some glueball propagators for which the anomalous dimensions or the powers 
of the gauge coupling have been factored out. 

12.8 Glueball propagators in the canonical scheme 

If we wish to recover the anomalous dimensions we should choose a canonical scheme, in 
which the fields are normalized in such a way to include the renormalization factors. 

Localization of twistor Wilson loops is just a statement about the homology of 1 and in princi- 
ple does not provide a dictionary to identify fluctuations of surface operators with specific glueball 
propagators in perturbation theory. However, heuristically we can construct a dictionary on the 
following basis. We notice that the one-loop anomalous dimensions in the ground state of the 
Hamiltonian spin chain coincide within one-loop accuracy with the anomalous dimensions com- 
puted by means of the localization on surface operators by the change to the ASD variables: 

95 1 ,A, , n^ L , olO 1 A,L 

Indeed, choosing a canonical scheme and following the definitions of sect. (11), we rescale the 
cutoff by a factor of Z2 and the area L^ by a factor of ^ . Thus A = Z2 Ac and t} = g^L^.. This 
changes the normalization of the operators insertions by factors of Z and the normalization of the 
effective action of the fluctuations by a power of g^, since the effective action is quadratic in the 
area of surface operators (Eq.(12.44)). Hence we get: 

A 






^Z2Acx8L-8X'r/2Ar2^4 



In 



J < trNTrfj,{d^^^-^fl^){x+,X-)trNTr^,{^^fi^-^ 5 fi){0,0) >^cln /''+^-+''-^+)^^x (12.79) 

where the canonical expectation value for the fluctuations is computed with respect to the effective 
action in Eq.(12.44) with L replaced by L^. It follows that: 

(^)8L-8^/2^2p4| ^ ;r^rr^,(5^^^-ip^)(x+,x_)?r^rr^,(Ai^Ai^"'5Al)(0,0) >iSl e^^^-'-^^-'-^d^x 
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2k 
I < trNTrj^,{5^n^-^ll^){x+,x^)trMTrj^,{H^fi^-^5ll){0,0) >i^i e'^P^'-+P-'^U^x (12.80) 

Thus the perturbative anomalous dimensions for long operators in the ground state of the Hamil- 



tonian spin chain in the thermodynamic limit [20, 21] are correctly reproduced by correlations of 



long surface operators in the canonical scheme. Actually they agree also for L = 1 , since in this 
case the anomalous dimension is determined by the beta function via the factor of g^. This suggests 
also that the states which surface operators factorize on by homological localization are all scalars, 
although in principle they may couple also to tensors since /i/I is not a scalar. Analogously in the 
canonical scheme the same anomalous dimensions arise for the operators O^: 



< Tr^rO^{p+,p-)Tr,yrO^{-p+,-p-) >conn 



lL(_„ _„ N^W 



= g\-p+P-)Z-—{-p+p-) < Trj.O^{p+,p^)Tr,rO^{-p+,-p^) >'conn 

^ g\-p^p_)Z-'^{-p^p^){-p^p^f^-Hog^P^ (12.81) 

13. Wild local systems and Regge trajectories 

The hyper-Kahler reduction induced by the restriction to local systems can be extended to 
representations of the wild fundamental group [ |85| , ^ . Also the extension to the wild surface 
operators admits a gauge in which the theory is locally abelian, because of the commutativity of 
the coefficients of the higher order poles [|8^, ^^ . Hence in principle we can extend the computation 
of the fluctuations of surface operators in Eq.(12.39) to curvatures that involve derivatives of the 
delta function: 



P P n 



P \ 

+ EE5Aif(M,M)5"5P)(z-z^(i.,i7)) (13.1) 

p n 

This corresponds naturally to Regge trajectories of higher spins. We leave the computation for the 
future. 

14. Conclusions and outlook: QCD-like theories and the twistorial string theory 

The main conclusion of this paper is that there exist twistor Wilson loops that can be localized 
in \wge.-N pure YM on local systems, i.e. on representations of the fundamental group of a punc- 
tured Riemann sphere immersed in space-time, i.e. on surface operators. These surface operators 
turn out to be connections with Zj^ holonomy around the punctures. The localization on surface 
operators leads to the one-loop exactness of the \wg&-N Wilsonian beta function and to a canonical 
beta function of NSVZ type. 

Some understanding of the mass gap and of the glueball spectrum occurs in a certain sector of 
the theory associated to twistor Wilson loops. 
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By certain changes of variables, that imply integrating on the moduli of surface operators, 
the loop equation for twistor Wilson loops is written in a holomorphic gauge in which a non- 
trivial glueball potential is generated by the change of variables. The second derivative of the 
glueball potential implies a mass term for the glueballs that is non-vanishing precisely for the 
surface operators with degenerate eigenvalues that occur for Z^ holonomy. 

In this language glueballs arise as massive fluctuations of magnetic surface operators supported 
on the Lagrangian submanifold of space-time that is the support of the twistor Wilson loops. 

This is the picture that follows from the localization of the loop equation for twistor Wilson 
loops and that realizes a new version of some long-standing ideas about dual superconductivity in 
pureyM[jT^,0,|l8]l. 



On the field theory side we may wonder as to whether the methods of this paper extend to ^ 
= 1 SUSY YM, once we observe that twistor Wilson loops are not in the Parisi-Sourlas cohomology 
of the Nicolai map of ^ = 1 SUSY YM and therefore there is no reason for which they should be 
localized on instantons. In fact from the point of view of the standard folklore the glueball spectrum 
of \aige-N jV = 1 SUSY YM should not differ in a qualitative way from the one of YM. 

Another extension would be to QCD-like theories, such as YM minimally coupled to Nf mass- 
less Dirac fermions in the fundamental representation in the large-A/^ limit, keeping the ratio -^ 
fixed. The computation of the associated glueball spectrum would imply the determination of the 
lower side of the conformal window, as the point at which the mass gap disappears. 

The basic issue involved in such extensions is the realization of the correct Wilsonian beta 
function around the critical points provided by the localization, and involves a crucial understand- 
ing of the fermion zero modes in a neighborhood of the critical points. 

On the string theory side the results of this paper suggest a new string program for the YM 
theory, if we look for exact solvability. 

It has been known for some time that .jY = 4 SUSY YM admits a partially equivalent twistorial 



string [ ]99| , |123| , |124[ , |125| ]. The triviality of twistor Wilson loops and the fact that they are supported 
on Lagrangian submanifolds in twistor space suggests the existence of a stringy interpretation of 
our results in terms of open topological strings ending on Lagrangian submanifolds in twistor space 



[ ]12q , |127| ] in presence of surface operators. The occurrence of topological strings, as opposed to 



the usual strings, is due to the trivial nature of twistor Wilson loops at large A'^. 



Since such twistorial topological string would be solvable by cohomological localization [ ]128 
|129| ], morally this conjectured topological string/gauge theory duality would provide the string 
cohomology dual to the field theory homology. Some hints about this conjectured twistorial string 
of YM can be found in [113 



We may wonder what such a twistorial string theory would be suited for, since it is sup- 
posed to be equivalent to the field theoretical results of this paper. The answer is found in the old 
fashioned unitarization program of string theory. The field theory computation provides the free 
glueball spectrum in the twistor sector but does not furnish easily information about the glueball 
interactions. However, in the string approach interactions are fixed by the geometry of the string 
world-sheet, once the free theory is known, and thus the conjectured gauge theory/ topological 
string duality would open the way to computing the glueball 5-matrix. 
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